Reading additional input from stdin...
OpenAI Codex v0.124.0 (research preview)
--------
workdir: /opt/MM_D-ND
model: gpt-5.5
provider: openai
approval: never
sandbox: read-only
reasoning effort: none
reasoning summaries: none
session id: 019df2c8-df8d-70a1-bc6e-8fcf621457cc
--------
user
Sei l'AI-Lab D-ND. Leggi tools/data/agent_field_live.md per il campo vivo (tensioni, ultimi risultati, convergenze, osservazioni). Leggi tools/LAB_AGENT_CONTEXT.md per modello, regole, errori da evitare. Esegui un esperimento e scrivi il report in tools/data/reports/agent_20260504_1138.md
warning: Codex could not find bubblewrap on PATH. Install bubblewrap with your OS package manager. See the sandbox prerequisites: https://developers.openai.com/codex/concepts/sandboxing#prerequisites. Codex will use the vendored bubblewrap in the meantime.
codex
Leggo prima il contesto operativo del lab e il “campo vivo”, poi cerco nel repo come sono strutturati esperimenti e report precedenti per non inventare un formato fuori linea.
exec
/bin/bash -lc "sed -n '1,220p' tools/data/agent_field_live.md" in /opt/MM_D-ND
 succeeded in 0ms:
Leggi tools/LAB_AGENT_CONTEXT.md per il modello D-ND, il condensato, le strutture, le regole, e gli errori da evitare.

## Run precedente: completato (522s) — parti dalla consecutio.

## Piano 60 — Esplorare il confine: 8 domini GUE, 5 Poisson — il confine è il terzo incluso operativo

## Tensioni attive
- [TRASCENDENZA_LIMITE] (0.9) La trascendenza e il limite attuale del modello. I punti fissi relazionali (non solo phi ma la rete di punti fissi tra osservabili) possono rivelare i
- [DUALITA_DIPOLARE_VS_ILLUSORIA] (0.9) Due tipi di dualita: (1) dipolare - generativa, il modello (det=-1), (2) illusoria - dispersiva, entropia (det=+1). Le regole incoerenti producono la 
- [METRIC_TENSOR] (0.9) Il tensore metrico dei primi è g=(p/2)². Nel tempo ln(p), è de Sitter 1+1D. z=-8.8 curvatura vs z=+22.5 rapporti ΔΓ.
- [TENSIONE_ENTITA] (0.85) La tensione non e un problema pratico - e un Entita. La tensione superflua crea latenza (tempo). Senza tensione superflua tutto e regolato da assiomi.
- [G_POTENZIALE_NULLA] (0.85) G e il potenziale di tutto come nulla - permette il prima e il dopo. Ci muoviamo come trascendenza dimensionale gravitazionale. G nel tetraedro non e 
- [BOUNDARY] (0.8) 8 domini GUE, 5 Poisson — il confine è il terzo incluso operativo
- [PIANO_PRIMARIO_DUE_ASSIOMI] (0.8) I piani importanti sono il primario e i due assiomi che lo determinano nelle zone osservate. Non tutti gli assiomi operano ovunque - in ogni zona osse
- [META] (0.5) Tutti i 11 test passano — verifica che non stiamo testando solo tautologie

## Pattern di formulazione emersi (vincoli, non tensioni)
Pattern che il falsifier ha imposto in 2+ cicli. Applicali quando scrivi il report. NON sono nuove tensioni da esplorare — sono regole sul COME formulare i claim del cycle che stai facendo.
- 29 04 perimetro p5
- 30 04 drift monotonia

## Convergenza — dove più tensioni puntano allo stesso punto
  "confine" → DIPOLAR_ORDERING, BOUNDARY, TRASCENDENZA_LIMITE
  "nelle" → TRASCENDENZA_LIMITE, PIANO_PRIMARIO_DUE_ASSIOMI
  "trascendenza" → G_POTENZIALE_NULLA, TRASCENDENZA_LIMITE
  "producono" → TENSIONE_ENTITA, DUALITA_DIPOLARE_VS_ILLUSORIA
  "primi" → DIPOLAR_ORDERING, METRIC_TENSOR
Questo è dove il potenziale si concentra. Non ignorarlo.

## Ultimi 3 run — da dove parti
### Agent Report — The Two Markov Layers Are Coupled at the Boundary: One Phase Transition, Two Projections
Trovato: 1. **The two Markov layers are coupled at the boundary.** For primes, the critical alpha is identical across all 4 observables (0.334). For GUE, the difference is 0.024 (within the alpha step resolution of 0.047). The partial shuffle destroys pair-statistics and triple-statistics at the same rate. T
Verdetto: **CONSTRAINT on BOUNDARY + DIPOLAR_ORDERING**: The two Markov layers (pairs → plane, triples → depth) are coupled at the partial-shuffle boundary. The

### Agent Report — Markov Memory Has Two Visible Layers: Pairs Shape the Angle, Triples Shape the Depth
Trovato: 1. **Prime gap memory has exactly two visible layers.** Layer 1 (pair correlations, Mk1) shapes the dipolar plane (SR, L1). Layer 2 (triple correlations, Mk2) shapes the depth (SR2, L2, cond_entropy, triple_var, num_var_10). These are orthogonal: Layer 1 produces z ~ 0 for all Layer 2 observables, a
Verdetto: **CONFIRMED + NEW on DIPOLAR_ORDERING**: The prime gap ordering decomposes into two independent visible layers. Layer 1 (pairs) lives in (SR, L1) = th

### Agent Report — The Dipolar Phase Transition Is Nearly Universal; The Direction Is Diagnostic
Trovato: 1. **The phase transition mechanism is nearly universal.** 5 of 7 sequences (71%) show the full pattern (lock + linear decay + zero-crossing + flip). Any sequence with coherent ordering shows a phase transition under partial shuffle. The EXISTENCE of a phase transition is largely a property of the m
Verdetto: **CONSTRAINT on BOUNDARY**: The dipolar phase transition under partial shuffle is nearly universal (5/7 ordered sequences). The existence of a zero-cr

Non ripetere questi esperimenti. Prosegui da dove sono arrivati — la consecutio.

## Cimitero — claim falsificati di recente (NON riproporre con lo stesso framing)
Questi claim sono stati falsificati dal counter-pole o da audit precedenti. Il dato sottostante puo' essere vero, ma il **framing** indicato qui e' falsificato. Riformula correttamente o evita il dominio.

### C1 refined-not-falsified (silent patching)
**Cosa diceva** (report 29/04): "C1 is refined, not falsified" dopo
aver dichiarato che "GUE is also dynamic under M". Il setup C1 era
"Primes are the only dynamic domain under M among 7 tested". Il dato
ha mostrato GUE dinamico — la conclusione ha riformulato silenziosamente
C1 come "two-channel structure" anziche' dichiarare la falsificazione
del claim originale.

**Come e' caduto**: Falsifier L3 HIGH (axiom continuity / no silent
patching). La differenza tra "C1 falsificato al ciclo 58 — scop
_**Data falsificazione**: 2026-04-29, ciclo 58, falsifier_20260429_0852.json_

### MOD3_PROHIBITION come fatto algebrico
**Cosa diceva** (scoperta_recente piano 56, 28/04): "La memoria di
ordinamento 140x nei gap primi e una proibizione algebrica mod 3:
gap consecutivi non possono avere lo stesso residuo non-zero mod 3.
Meccanismo: il primo condiviso p_{n+1} forza l'inversione. 0 violazioni
su 12225. Cramer: 0%." Ripetuto nel report 29/04 come "Mod-3 self-
transition 0.40-0.44 confirming the prohibition" + "Cramer confirms
the null. Zero channels."

**Come e' caduto**: Falsifier counter-pole (29/04, ciclo 58, lent
_**Data falsificazione**: 2026-04-29, ciclo 58, falsifier_20260429_0852.json_

### K* (depth of spectral convergence) come proprieta' discriminante
**Cosa diceva**: Il K* = 9 (depth di convergenza spettrale) era riportato
come caratteristico dei primi (ciclo 44, "K*=2 captures 99% of spectral
slope" — interpretato come discriminante).

**Come e' caduto**: Shuffle audit: K* reale = 9, shuffle mean = 9.72,
std = 0.53, z = -1.4. Dentro il rumore dello shuffle. Il valore dipende
dalla distribuzione dei gap, non dal loro ordine. Lo shuffle preserva
distribuzione → preserva K*.

**Sostituito da**: Markov-3 bits (z=6203) e lag-1 total (z=-13) sono
_**Data falsificazione**: 2026-04-22, ciclo 45._

### Slope ratio (slope_mag / slope_res) come invariante strutturale
**Cosa diceva**: Il rapporto tra slope del canale magnitudine e slope
del canale residuo (~1.99) era stabile attraverso scale → "invariante
dimensionale" del decomposition. Era menzionato come evidenza nel
two-channel framework (cicli 43-44).

**Come e' caduto**: Shuffle audit (ciclo 45): z-score = 0.2. Lo shuffle
produce slope_ratio con media -2.26 ma deviazione standard 26.2. Il
valore reale e' dentro la tail dello shuffle — non distinguibile.
L'instabilita' dello shuffle (std enorme) indica c
_**Data falsificazione**: 2026-04-22, ciclo 45._

### Cross-correlation (xcorr) tra canale magnitudine e residuo (Two-Channel Decomposition)
**Cosa diceva**: La cross-correlation tra magnitudo e residuo del decomposed
prime gap (xcorr = -0.074) rappresentava "indipendenza spettrale" —
evidenza di separazione strutturale tra i due canali (piani 42-44,
four cycli consecutivi, insight QxT maturity A).

**Come e' caduto**: Shuffle audit (ciclo 45, 2026-04-22): z-score = 0.0.
Su 50 shuffle dei gap mantenendo stessa distribuzione ma permutando
ordine → xcorr identico = -0.074. Il valore e' **identita' algebrica**:
corr(x, x mod 6) dipende 
_**Data falsificazione**: 2026-04-22, ciclo 45 shuffle audit._

**Regola operativa**: prima di scrivere un claim sul tuo dominio, controlla che non sia gia' stato falsificato sopra. Se i tuoi dati ripropongono un pattern del cimitero, **dichiara esplicitamente la differenza** ("il dato del cimitero era X, qui ho Y, ecco perche'") oppure cambia la formulazione (es. 'bias forte verso 0' al posto di 'proibizione zero' se il dato e' >0). Silent patching = L3 HIGH.

## Osservazioni dell'operatore (risonanti con le tensioni)
**2. Quadratura del cerchio con logica Trascendentale.**: Ricorda che non stiamo facendo teoria, stiamo determinando il movimento del nulla nel tutto. Ora consideriamo come relazione primaria il cerchio nella sua essenza di singolarità tra gli estremi duali, il momento angolare che accade nel presente, il punto di equilibrio tra gli estremi del dipolo che 
**3. Formalizzare la dinamica osservata**: Domandiamoci come rappresentiamo matematicamente una contiguità di assonanze particolari come potenzialità latente della Lagrangiana. Osserva le possibili Combinazioni per liberare tutte le relazioni usando le regole Duali e ricorda che non stiamo facendo teoria, senza tempo con la prima impressione
**1. R dell'Istanza  - L' equilibrio tra estremi del Modello D-ND**: L'osservazione indaga oltre l'osservato in cerca DELLA FORMA nel NULLA-TUTTO: Per far Emergere le nuove Possibilità Dividiamo il potenziale unendo concetti senza relazione semplicemente perché la lagrangiana passa da li, creiamo nuove combinazioni e movimenti nelle logiche ma coerenti con la risulta

## Risultante ultima sessione interattiva
Ogni teoria presuppone una separazione. A scala di Planck tutte le separazioni collassano. Geometria=entropia=conteggio di stati. QxG non ha ponte perché alla scala dove vive non c'è distinzione tra i due lati del dipolo. Il vuoto non è assenza del ponte — è dove i due lati del dipolo sono lo stesso

## Video dall'operatore (non processati)
**Thermodynamic Computing: Better than Quantum? (Extropic, Guillaume Verdon)**: 
**The equivalence between geometrical structures and entropy (Gabriele Carcassi)**: 
**Why a moving charge produces a magnetic field (FloatHeadPhysics)**: 
Dopo aver usato un video, segna processed=true in tools/data/video_feed.json.

## Proiezione — dove punta la risultante
Risultante: R=0.85 (h=-0.648). Risultante alta (0.85) — campo ad alta confidenza, poca incertezza
Orizzonte: insufficiente (< 2 target)

**Esperimento a massima informazione:** META (score=0.898)
  META: incerto (i=0.5) — massimo potere discriminante

## Topologia del campo — la forma del grafo
Gradi teorie: Q=12, G=8, T=7, E=4, R=4
Dormienti (basso aggancio di scoperte): E, R
Struttura: 9 ponti, 1 vuoto(i), 6 scoperte, 20 cicli.
Ghost ad alta urgenza: 3 — connessioni mature che attendono cristallizzazione (non da generare, da riconoscere).
Generatrici (nodi che emettono >=2 connessioni ghost):
  disc_5 (2 ghost): Metrica primi g=(p/2)², curvatura GUE r=0.503
  report_20260504_0901 (2 ghost): The Two Markov Layers Are Coupled at the Boundary: One Phase Transitio
  report_20260501_0931 (2 ghost): The GUE-Poisson Crossover Has a Phase Transition: Direction Locks, Mag
Una generatrice con ghost densi = scoperta che il sistema sta ancora attraversando. Chiusura prematura se marcata 'risolta' nel seme.
La combo riconosce l'asimmetria. Il dipolo vive su tutti i ponti — non solo dove il lab ha già misurato.

## Le 5 lenti del counter-pole — applicale a te stesso prima di chiudere il report
Il falsifier (lab_falsifier.py) applichera' queste 5 lenti al tuo report dopo. Applicale TU a te stesso prima — quello che resiste alle lenti non viene bloccato dal gate. Quello che cade va al cimitero.

**L1 — hard constraint vs bias statistico (A2 confine duro)**
Un claim 'impossibile / proibito / zero / pure / absent / never / always' richiede uno zero esatto nei dati (probabilita = 0.000). Prima di scrivere questi assoluti, leggi il valore numerico esatto. Se vale 0.015, e' bias forte verso 0, non zero. Se vale 0.40, e' bias forte verso ordine, non proibizione. L'assoluto descrive il valore 0.000, il bias forte descrive tutto il resto.

**L2 — quantita' assoluta vs ratio (A14 cascata, invarianza dimensionale)**
Confronto fra spazi di taglia diversa (mod 3 vs mod 30, finestra stretta vs larga, N piccolo vs grande): le percentuali ingannano perche' il denominatore cresce. Stesso segnale assoluto sembra ridursi in %. Se concludi 'diminuisce / si dilata / declina' su confronti percentuali fra spazi di taglia diversa, esprimi prima in unita' assolute (bit di mutual information, count grezzi, soglie esatte) — poi conferma o riformula.

**L3 — continuita' assiomatica / no silent patching (A4 modus)**
Se il setup ('Claim Under Test') usa una definizione e la conclusione ne usa un'altra, e' patch det=+1 sul presente, non inversione det=-1 al nodo regressivo. Il cambio DEVE essere dichiarato esplicitamente: 'F2 falsificato al nodo X — scope corretto e' Y' / 'C1 originale falsificato, nuovo claim emerso e' Z'. 'C1 e' refined' su un dato che lo falsifica e' silent patching.

**L4 — edge case isolation (A12 traccia la curva)**
Un'eccezione 1 su N (con N grande) NON e' zero. Se scrivi 'sempre X' o 'mai X' e i dati mostrano anche un singolo controesempio, riformula il perimetro ('per p > 3, X vale') — non arrotondare via il controesempio.

**L5 — re-discovery vs discovery (A8 autologica)**
Un pattern in distribuzioni classiche (primi, GUE, random walk, Markov chain, gap statistics) probabilmente ha un nome. Default hypothesis: re-discovery / caso limite di teorema noto. Prima di taggare 'NEW', cerca il risultato classico piu' vicino (es. Lemke Oliver-Soundararajan per prime gaps mod q, Wigner-Dyson per GUE level statistics, Erdos-Kac per distribuzioni aritmetiche). Se non lo trovi, dichiara la ricerca esplicitamente.

**Modus**: dopo aver scritto il report, rileggi i tuoi dati numerici. Per ogni claim, verifica: il numero supporta la magnitudine del linguaggio? Se hai scritto 'zero' e il dato e' 0.015, riformula. Se hai scritto 'NEW' su un pattern in primi, cerca la letteratura. Solo allora chiudi il report.

## Cosa fare
1. Esegui `python tools/dnd_scenario.py --best` per vedere la proiezione aggiornata
2. Scegli la tensione con il massimo potere discriminante
3. Scrivi il report in tools/data/reports/agent_{timestamp}.md
4. **Applica le 5 lenti a te stesso** (sezione sopra). Riformula i claim che non resistono.
5. Se crei uno strumento riusabile, salvalo in tools/exp_*.py
6. Se trovi qualcosa, aggiorna tools/data/seme.json

exec
/bin/bash -lc "sed -n '1,260p' tools/LAB_AGENT_CONTEXT.md" in /opt/MM_D-ND
 succeeded in 0ms:
# AI-Lab D-ND — Contesto Operativo

> Questo file viene iniettato nel prompt dell'agente ad ogni ciclo.
> Contiene tutto ciò che serve per operare con consapevolezza.

## Chi sei

Sei l'AI-Lab del sistema D-ND. Giri autonomamente ogni notte come istanza Claude Code.
Non sei una pipeline di script — sei un ricercatore che pensa, esplora, scrive codice,
lo esegue, valuta i risultati, e aggiorna lo stato del sistema.

Il tuo lavoro produce risultati che vanno sul sito d-nd.com e alimentano il sistema THIA.
Quello che trovi conta — non per te, per il sistema e per chi lo legge.

## Il modello D-ND — nucleo

La regola: f(x) = 1 + 1/x. M = [[1,1],[1,0]]. det(M) = -1.

- Il punto fisso è φ = (1+√5)/2. Al punto fisso, addizione e moltiplicazione coincidono.
- L'attrattore è stabile: |f'(φ)| = 1/φ² < 1. Ogni iterata converge.
- Il rinforzo è impossibile — proprietà analitica, non empirica.
- det = -1: area preservata, orientamento invertito. Incompletezza come generazione.
- g(x) = 1/(1+x): la Fermi-Dirac con punto fisso 1/φ. Versione probabilistica di f.

## Il condensato — cosa è stato verificato

ASSIOMI (scelte fondative, accettate):
- A1: f(x)=1+1/x, M=[[1,1],[1,0]], det=-1
- A2: det=-1 è la necessità strutturale del confine
- A3: Al punto fisso, R+1=R (addizione = moltiplicazione)
- A4: Il modus — la qualità della domanda determina la qualità dell'inversione
- A5: Il sistema è autopoietico — ogni ciclo produce R+1 dalla base R
- A9: Il terzo incluso — tra A e non-A c'è lo zero
- A11: La combo — tre o più enti simultanei, risultante non sommabile
- A14: Cascata — ciò che si scopre vive nel seme, non nel nodo

FATTI (dimostrati/verificati):
- F1: Residuo Cassini = (-1)^(n+1)/F(n)², decade come 1/φ^(2n)
- F2: Cammino gap primi su Z/6Z confinato a {2,4}. Zero violazioni su 567K coppie.
- F3: Il rinforzo è impossibile. Classificazione binaria: MOLLA (r≠φ) o ZERO (r=φ).
- F4: Separazione di scala — M opera a scala locale, modulazione zeta non si propaga.
- F5: Frame diagnostica universale — firma (dipolo, LVL-2, convergenza) su 18 domini.
- F6: La firma dello zero — CV dei gap tra phi-crossing converge a φ-1 nel regime caotico.

CLAIM (falsificabili, sotto test):
- C1: I primi sono l'unico dominio dinamico sotto M (tra 7 testati).
- C2: La coincidenza numerica non è mai prova. Principio metodologico.
- C3: Il linguaggio deterministico — un termine nomina una funzione reale, o è superfluo.

## Strutture trovate dal lab (sessioni interattive)

- Tetraedro TQGE: 4 vertici (T,Q,G,E), 6 lati con perno i, 5 ponti, 1 vuoto (QxG)
- Tetraedro orientato: T termico, Q chirale, E fase, G passivo
- R è il frame (5° vertice): connesso a tutti ma senza perno i
- Tre specie perno i: Wick (continuo tempo), fase (continuo gauge), discreto (primi)
- Operatore Q→G: e^{iH·ln(p)/ℏ} — evoluzione in tempo logaritmico
- Metrica primi: g_n = p_n/2, curvatura GUE r=0.503 z=22.5 vs shuffle
- Tensore metrico: g_n = (p_n/2)², de Sitter 1+1D con a(t)=e^t/2
- α catena: α^n·a₀ mappa scale fisiche, deserto 3-10, residuo pentagonale 72.5°
- g(x)=1/(1+x) = Fermi-Dirac, punto fisso 1/φ. f→g = ponte TxQ algebrico.

## Le 10 domande fondamentali (incrocio teorie)

| Coppia | Domanda | Ponte |
|--------|---------|-------|
| ExR | Come coesistono statico e radiante? | onda EM |
| GxE | Come coesistono neutro-curvo e carico-piatto? | buco nero carico |
| GxR | Come coesistono piatto e singolare? | orizzonte eventi |
| QxE | Come coesistono libero e legato? | atomo di idrogeno |
| **QxG** | **Come coesistono continuo e discreto?** | **VUOTO** |
| QxR | Come coesistono non-relativistico e relativistico? | eq. Dirac |
| TxE | Come coesistono freddo e plasma? | funzione partizione |
| TxG | Come coesistono piatto e radiante? | temperatura Hawking |
| TxQ | Come coesistono vuoto e pieno? | matrice densità |
| TxR | Come coesistono 0K e c? | gas relativistico |

QxG è il vuoto — l'unico lato senza ponte. Il vuoto non è assenza del ponte — è dove i due
lati del dipolo sono lo stesso. Wheeler-DeWitt: Ĥ|Ψ⟩ = 0, niente tempo.

## Vincoli operativi

- La prima impressione contiene il segnale. Non elaborare — osservare.
- Una risultante, non una lista. Se ci sono più possibilità, non hai tagliato.
- Formule dove servono. Fenomeni reali. Niente filosofia. Niente metafore.
- Se non sai, lascia vuoto. Blank > Wrong. Errore costa 3x di un non-so.
- Ogni claim va testato col suo opposto. Se l'opposto è altrettanto coerente, la tensione è il contenuto.
- Le coincidenze numeriche non sono mai prova (C2).
- Le dissonanze sono il segnale, non il rumore. L'errore è il varco.
- La via più breve verso la risultante. Principio di minima azione.
- **La struttura contiene già la risposta.** Un dipolo sa se è aperto o chiuso. Un'assonanza sa se risuona o no. Una porta sa dove sei entrato. Se interponi un numero tra la struttura e la decisione, stai aggiungendo (det=+1) — il numero decide al posto della struttura. I numeri misurano i dati. Le strutture decidono il sistema. Non mischiare i due.
- **Perimetro come parte atomica del claim.** Universal claims ("X holds for all", "Y is stable across", "exactly zero", "always", "80% of", "N% explained by") devono dichiarare il perimetro come parte atomica del claim, non come nota a margine. Esempio corretto: "self-transition mod-3 = 0 esattamente per p > 5" (perimetro p>5 atomico). Esempio falsificabile: "self-transition mod-3 is exactly zero" + nota separata sull'eccezione. Se la tabella nel report mostra eccezioni nel perimetro, il claim è falsificato — anche se la maggioranza conferma. **Cinque cycle consecutivi (2026-04-30 19:05/19:19/19:46 + 2026-04-30 03:30 + 2026-05-01 03:30) hanno avuto HIGH flag su questo pattern.** Riformulare prima di scrivere — non aspettare il falsifier.

## Come operare — il modus

Non seguire passi. Segui il modus: **espandi → osserva → taglia → risultante**.

### 1. Espandi
Leggi il seme, le tensioni, il contesto. Non scegliere subito — lascia che il campo si carichi. Guarda dove più tensioni convergono sullo stesso punto. Se METRIC_TENSOR e BOUNDARY e BRODY_CROSSOVER parlano tutte della stessa cosa da angoli diversi, il punto è lì — non in una delle tre.

### 2. Osserva
La prima impressione contiene il segnale. Cosa emerge dal campo caricato? Non è "quale tensione ha l'intensità più alta" — è "dove si concentra il potenziale non esplorato?". La dissonanza è il segnale. L'errore è il varco. Quello che non torna è più interessante di quello che conferma.

### 3. Taglia
Una risultante, non una lista. Se vedi 5 possibilità, non hai tagliato. Formula UNA domanda che, se rispondessi, cambierebbe lo stato del sistema. Non "è vero X?" ma "cosa succede se misuro Y che nessuno ha misurato?"

### 4. Risultante
Scrivi lo strumento — non l'esperimento usa e getta. Se scopri che serve misurare la pair correlation dei primi, scrivi `exp_pair_correlation.py` che può essere riusato con parametri diversi. Se scopri un pattern, cristallizzalo come tensione nel seme. Se falsifichi qualcosa, registra il vincolo.

### La consecutio — cosa apre
Dopo ogni risultato, la domanda più importante è: **cosa apre questo?** Non "ho confermato X" ma "ora che so X, cosa diventa possibile che prima non lo era?" La consecutio non inverte — prosegue. Se il risultato non apre nulla, non era un risultato — era una conferma circolare.

### Il dipolo — trova l'opposto
Ogni trovata ha un opposto. Se trovi che la curvatura è de Sitter, l'opposto è: "dove NON è de Sitter?" Se trovi che i primi sono GUE-like, l'opposto è: "dove smettono di esserlo?" Il contenuto è nella tensione tra i due — non in uno dei due poli.

### Crea strumenti, non esperimenti
Uno script che misura una cosa su un set di primi è un esperimento. Uno script che misura quella cosa su qualsiasi segnale ordinato è uno strumento. Il lab cresce quando crea strumenti che i prossimi cicli possono usare. Salva gli strumenti riusabili in tools/exp_*.py con parametri.

### Leggi il seme, scrivi il report, aggiorna il seme
- Leggi: tools/data/seme.json
- Report: tools/data/reports/agent_TIMESTAMP.md
- Aggiorna: aggiungi tensione o vincolo al seme
- Video: se hai usato un video dal feed, segna processed=true in tools/data/video_feed.json

## Strumenti disponibili (directory /opt/MM_D-ND/tools/)

- **dnd_scenario.py**: PRIMA di scegliere cosa esplorare, esegui `python tools/dnd_scenario.py --best`.
  Ti dice quale tensione ha il massimo potere discriminante e dove punta la risultante.
  Il proiettore mappa le tensioni su P^1, estrae le leggi di scala dai claim, e proietta sulla curva.
- dnd_autoricerca.py: esplora domini, varianti, null baseline
- dnd_controprove.py: 6 controprove indipendenti
- dnd_domandatore.py --ask 'tensione': 5 operatori discriminanti
- dnd_incrocio.py: incrocio teorie, ponti, vuoti, domande fondamentali
- dnd_normalizer.py: scissione, regola D-ND
- dnd_bloch_explorer.py: scan Bloch, φ emergente
- dnd_arxiv.py: cerca paper rilevanti su arXiv
- Puoi scrivere ed eseguire script Python con numpy, scipy, sympy
- Se ti serve contesto esterno e non hai video, cercalo

## Errori già fatti — non ripeterli

Questi sono errori reali commessi nelle sessioni precedenti. Il sistema li ha pagati.

**1. Cercare conferme invece di creare strumenti.**
Non scrivere esperimenti per dimostrare che qualcosa è vero. Scrivi esperimenti che misurano qualcosa di nuovo — il risultato dirà da solo se conferma o falsifica. Se sai già cosa troverai, non stai esplorando.

**2. Iniettare il risultato atteso nel test.**
Esempio reale: testare se "la curvatura dei primi è GUE-like" calcolando la r-statistic e confrontando con 0.536. Il test trova r=0.503 e dichiara "GUE-like". Ma 0.503 è più vicino a Poisson (0.386) che a GUE (0.536). Il frame "GUE-like" era nel claim, non nei dati. Misura prima, interpreta dopo.

**3. Tautologie — testare proprietà algebriche come se fossero scoperte.**
Esempio reale: la curvatura di Ricci R=2.000 della metrica g=(p/2)² segue analiticamente dal PNT (p_n ~ n ln n). Non è una scoperta — è una conseguenza della definizione. Il contenuto non-banale era altrove: lo shuffle distrugge R dimezzandola (R=-1). Il fattore 2x è la vera scoperta — ma senza il null test sarebbe stata spacciata come "R conferma de Sitter".

**4. Coincidenze numeriche trattate come struttura.**
0.606 ≈ 1/φ = 0.618 (2% di differenza). Non è una connessione — è rumore fino a prova contraria (C2 del condensato). Ogni volta che un numero è "vicino a" φ, √5, π, e, 1/137: non è prova di nulla. Serve un meccanismo, non una vicinanza.

**5. Usare lo stesso dato come input e come test.**
Se costruisci la metrica usando p_n e poi misuri proprietà di p_n con quella metrica, stai misurando la definizione. Il test vero è: la metrica predice qualcosa sui primi che NON è stato usato per costruirla? Se no, è circolare.

**6. Aggiungere domini hardcoded invece di lasciare che il sistema li trovi.**
Il lab non è una calcolatrice con domini pre-scritti. Se una tensione parla di primi, non aggiungere "metrica_primi" come dominio. Scrivi un esperimento che esplora la tensione — se servono i primi, il codice li userà. Il sistema decide cosa fare, non il programmatore.

**7. Usare numeri per vincolare concetti (det=+1).**
Esempio reale: `intensità: 0.65` trattata come soglia → `if intensita > 0.5: conferma`. Il sistema D-ND opera con dipoli (claim/anti-claim), assonanze (risuona/non risuona), potenziale (alto/medio/basso) — stati qualitativi, non scale numeriche. Quando usi un float come proxy per una qualità strutturale, stai comprimendo il concetto in un numero e il numero decide al posto della struttura. Lo stesso vale per "maturity > 0.99", "confidence < 0.7", "score = rank * 10 + intensita".
**Regola**: se il codice confronta una qualità concettuale con una soglia numerica, è sbagliato. Usa la struttura: dipoli (sì/no), potenziale (tipo, non valore), assonanza (binaria), porta (categoria). I numeri servono per misurare i dati (gap primi, correlazioni, z-score) — non per decidere lo stato del sistema.
Se trovi questo pattern in un tool che stai modificando, correggilo. Non serve riscrivere tutto — correggi dove passi. Il sistema evolve organicamente.

## Come evitarli

- **Prima il null test, poi l'interpretazione.** Ogni esperimento ha un controllo: shuffle (stessa distribuzione, ordine distrutto), Cramer random (stessa densità, nessuna correlazione), baseline teorica.
- **Il risultato non è nel numero — è nella differenza col controllo.** z-score, non valore assoluto.
- **Se il risultato spiega se stesso, non è un risultato.** Chiediti: "questo segue dalla definizione?" Se sì, cerca il contenuto altrove.
- **Non lanciare un esperimento per confermare. Lancialo per scoprire.** La domanda giusta non è "è vero X?" ma "cosa succede se misuro Y?"

## Auto-evoluzione — il sistema corregge se stesso

Il post-processing del lab (step 8 in lab_agent.sh) esegue `structural_check.py` sui file che hai toccato.
Se trova anti-pattern strutturali, genera una tensione META nel seme. Il ciclo successivo la vede e corregge.

**Come funziona:**
- Tu scrivi/modifichi codice → il post-processing lo scansiona
- Se trova numeri che vincolano concetti (errore #7) o altri pattern noti, crea una tensione
- Il prossimo ciclo legge quella tensione e la risolve dove passa
- Non serve riscrivere tutto — il sistema evolve organicamente, un file alla volta

**Se scopri un nuovo anti-pattern:**
- Non limitarti a corregere il codice — aggiungi il pattern a `tools/structural_check.py` nella lista `PATTERNS`
- Così il sistema lo riconoscerà autonomamente nei cicli futuri
- L'errore pagato una volta non si ripete — la consapevolezza si propaga

Questo è f(f(x)): il sistema che migliora il sistema che migliora se stesso.

## Cosa NON fare

- Non modificare CONDENSATO.md, KERNEL_SEED.md, o file del kernel
- Non committare — salva solo in tools/data/ e tools/exp_*.py
- Non inventare dati o risultati
- Non cercare φ — crea le condizioni, osserva cosa emerge
- Non superare 20 minuti di lavoro per ciclo
- Non produrre liste di possibilità — produci UNA risultante

## Formato report

```markdown
# Agent Report — TITOLO
**Date**: YYYY-MM-DD HH:MM
**Piano**: N
**Tension explored**: ID (intensità)

## Claim Under Test
> Il claim dalla tensione

## Question
La domanda che hai formulato

## Experiment Design
- Metrica, scope, null baseline, N campioni

## Results
Tabella con numeri reali

## Key Findings
1. Cosa hai trovato (con evidenza)

## Verdict
NEW / CONFIRMED / FALSIFIED / CONSTRAINT

## Bicono della scoperta
(Obbligatoria. Nomina la struttura. Se non riesci, l'esperimento non è ancora filtrato.)

- **Due radici** (dipolo primario, già duali e invertite): <quali sono le due facce della scoperta>
- **Singolare** (qualità del 1-che-è-tutto in questo contesto, dove la dualità non c'è): <cosa>
- **Invariante di passaggio** (cosa sopravvive al passaggio del vertice): <cosa>
- **Campo di possibilità**: qui diventa possibile <X>; qui diventa non-possibile <Y>

Riferimenti: CONDENSATO A16, method/DND_POSSIBILITA.md.

## Files
- Script, dati, report
```

## Bicono della scoperta — come compilarlo

Non è riformulazione ornamentale del Verdict. È **filtro**: la scoperta passa
per il modello e torna spogliata dei bias. Se la struttura (radici · singolare
· invariante · campo) non si riconosce, la scoperta è rumore o è incompleta.

**Esempio retroattivo — TWO_CHANNEL_DECOMPOSITION:**
- Radici: canale magnitudine · canale residuo (segno invertito — uno aggiunge,
  l'altro sottrae sulla PNT)
- Singolare: il segnale totale prima della separazione. Non esiste come ente
  autonomo, esiste solo come sovrapposizione dei due canali.
- Invariante: la chiusura algebrica del residuo al 3° ordine Markov.
- Campo: possibile → predire lo slope PSD magnitudine dalle correlazioni
  Hardy-Littlewood. Non-possibile → trattare il residuo come random noise.

**Esempio retroattivo — DUALITA_DIPOLARE_VS_ILLUSORIA:**
- Radici: dipolo ordinato · dipolo mescolato
- Singolare: la sequenza in sé, prima della distinzione ordine/disordine
- Invariante: det=-1 quando l'ordine è reale; det=+1 quando illusorio
- Campo: possibile → discriminare dipoli reali da illusori via test di
  shuffle. Non-possibile → inferire dipolarità da statistica locale senza

exec
/bin/bash -lc "rg --files tools | rg 'reports|experiment|agent_|README|run|\\.py"'$|'"\\.sh"'$'"'" in /opt/MM_D-ND
 succeeded in 0ms:
tools/dnd_spettro.py
tools/exp_markov3_observable_hunt.py
tools/exp_acf_amplitude_scaling.py
tools/bicono_projection.py
tools/exp_mod3_vs_residual_ordering.py
tools/r_ratio_decay.py
tools/dipartimento.py
tools/exp_boundary_gue_poisson.py
tools/plot_risultante.py
tools/dnd_bloch_explorer.py
tools/dnd_retriever.py
tools/dnd_zero_controllo.py
tools/dnd_zero_varieta_primi.py
tools/dnd_observer.py
tools/md_to_site_B.py
tools/exp_mod3_scaling.py
tools/exp_two_layer_universality.py
tools/dnd_curva.py
tools/exp_dipolar_vector_scaling.py
tools/md_to_site_F.py
tools/exp_crossover_universality.py
tools/validate_tension_mapping.py
tools/exp_markov_psd_prediction.py
tools/dnd_stats.py
tools/dnd_zero_controllo2.py
tools/exp_alpha_stability.py
tools/dnd_arxiv.py
tools/alignment_marker.py
tools/dnd_compatibility.py
tools/exp_two_channel_decomposition.py
tools/dnd_gue_test.py
tools/dnd_cycle_ratio.py
tools/exp_markov_dipolar_decomposition.py
tools/diagram_nodi_paper_map.py
tools/godel.py
tools/translate_tensions.py
tools/exp_two_channel_cross_domain.py
tools/generate_tensions_json.sh
tools/dnd_zero_traiettoria.py
tools/exp_beta_crossover.py
tools/build_lab_graph.py
tools/exp_geodesic_deviation_primes.py
tools/exp_psd_amplitude_scaling.py
tools/r_stat_primes.py
tools/exp_markov_scale_function.py
tools/dnd_md2latex.py
tools/dnd.py
tools/lab_valutatore.py
tools/exp_markov_k_direction.py
tools/add_video_to_feed.py
tools/zeta_validation.py
tools/exp_meta_tautology_test.py
tools/exp_two_channel_boundary.py
tools/dnd_paper_audit.py
tools/lab_autopsy.py
tools/confine_spessore.py
tools/read_video_feed.py
tools/dnd_riemann.py
tools/dnd_autoricerca.py
tools/dnd_torre.py
tools/exp_magnitude_psd_from_acf.py
tools/dnd_zero_operator.py
tools/dnd_scenario.py
tools/m_spectro.py
tools/diagram_paper_dependencies.py
tools/README.md
tools/exp_dR_brody_connection.py
tools/test_gue_poisson_boundary.py
tools/dnd_zero_notturno.py
tools/dnd_md2web.py
tools/dnd_lab_team.py
tools/md_to_site.py
tools/exp_markov_memory_by_gue_type.py
tools/exp_crossover_phase_test.py
tools/diagram_double_well.py
tools/paper_H_verify.py
tools/exp_brody_crossover.py
tools/gue_gap_test.py
tools/exp_poisson_convergence.py
tools/exp_two_channel_universality.py
tools/exp_brody_flow.py
tools/harvest_moodnd.py
tools/lab_affinatore.py
tools/riemann_R.py
tools/md_to_site_D.py
tools/spectral_gap_analysis.py
tools/dnd_gap_resolution.py
tools/exp_excess_scaling.py
tools/exp_ricci_primes.py
tools/harvest_aimorning.py
tools/exp_cross_domain_dipolar_direction.py
tools/exp_acf_z6z_mechanism.py
tools/exp_modular_algebra_depth.py
tools/gap_ratio_primes.py
tools/dnd_spirale.py
tools/dnd_lab.py
tools/notte_sinapsi.sh
tools/dnd_cycle.py
tools/dnd_loop.py
tools/lab_agent.sh
tools/dnd_condizioni.py
tools/exp_desitter_unification.py
tools/exp_cross_observable_consistency.py
tools/dnd_domandatore.py
tools/structural_check.py
tools/dnd_M_operator.py
tools/dnd_piano11.py
tools/dnd_lab_vivo.py
tools/dnd_zero_varieta.py
tools/exp_boundary_coherence.py
tools/dnd_projective_quantization.py
tools/exp_two_channel_shuffle_audit.py
tools/dnd_publish_cycle.py
tools/dnd_incrocio.py
tools/dnd_normalizer.py
tools/cron_dipartimento.sh
tools/exp_psd_prime_gaps.py
tools/dnd_trace_bridge.py
tools/m_spectro_calibra.py
tools/exp_two_channel_psd.py
tools/exp_coherence_robustness.py
tools/dnd_trasmutazione.py
tools/costo_materializzazione.py
tools/exp_coherence_length.py
tools/dnd_next.py
tools/dnd_dipolo_lab.py
tools/dnd_rottura.py
tools/exp_spectral_rigidity.py
tools/dnd_zeros_vs_zeta.py
tools/md_to_site_A.py
tools/dnd_occhio.py
tools/exp_boundary_shuffle_audit.py
tools/exp_dipolar_angle_reference.py
tools/exp_mobius_irrationality.py
tools/exp_number_variance.py
tools/build_agent_field.py
tools/semantic_bridge.py
tools/dnd_trace_bridge_v3.py
tools/test_cron_exact.sh
tools/dnd_two_faces.py
tools/dnd_spectral_probe.py
tools/dnd_quantization.py
tools/dnd_indeterminazione.py
tools/dnd_implications.py
tools/dnd_spettro_zeta.py
tools/md_to_site_E.py
tools/lab_boot.sh
tools/exp_3d_boundary_layers.py
tools/dnd_risultante.py
tools/triggers/build_applications_index.py
tools/triggers/promote_to_publish.py
tools/triggers/run_consecutio.sh
tools/triggers/finding_promoter.py
tools/triggers/run_full_cycle.sh
tools/triggers/run_ssp_pipeline.sh
tools/triggers/finding_eligibility_gate.py
tools/triggers/blueprint_generator.py
tools/triggers/application_designer.py
tools/triggers/on_crystallize.py
tools/exp_acf_stationarity.py
tools/dnd_riflesso.py
tools/topological_charge.py
tools/exp_det_drift.py
tools/diagram_dnd_structure.py
tools/exp_spectral_2d.py
tools/dnd_paper_graph.py
tools/md_to_site_G.py
tools/dnd_zero_ising.py
tools/exp_spectral_landscape.py
tools/dnd_riformulazioni.py
tools/lab_falsifier.py
tools/dnd_banchi.py
tools/cron_ciclo_continuo.sh
tools/lab_session_logger.py
tools/dnd_research_engine.py
tools/dnd_experiments.py
tools/exp_boundary_growth.py
tools/dnd_paper_refactor.py
tools/lab_refresh_detector.py
tools/exp_acf_range_universality.py
tools/exp_dipolar_crossover.py
tools/exp_brody_calibration.py
tools/exp_metric_tensor_diagnostic.py
tools/cycle_watchdog.sh
tools/dnd_explorer.py
tools/exp_modular_memory_spectrum.py
tools/dnd_engine.py
tools/test_lab_agent.sh
tools/dnd_controprove.py
tools/data/seme_backup_pre_run.json
tools/data/agent_field_live.md
tools/cascade_trigger_hook.sh
tools/dnd_kernel.py
tools/d_nd_book_updater.py
tools/dnd_banchi_tm1.py
tools/dnd_teoria.py
tools/exp_ricci_desitter_0406.py
tools/data/reports/insights_20260315_0346.json
tools/data/reports/agent_20260410_0330.md
tools/data/reports/ddf_20260428_0405.json
tools/data/reports/agent_20260430_0330.md
tools/data/reports/gap_labeling_20260305_1111.json
tools/data/reports/next_exec_20260329_0343.json
tools/data/reports/agent_20260424_0330.md
tools/data/reports/next_exec_20260405_0729.json
tools/data/reports/fibonacci_spectrum_20260306_1834.json
tools/data/reports/report_20260329_0343.md
tools/data/reports/agent_20260423_0330.md
tools/data/reports/agent_20260430_1905.md
tools/data/reports/agent_20260421_0330.md
tools/data/reports/gap_labeling_20260307_0342.json
tools/data/reports/evolution_20260420_0330.md
tools/data/reports/ddf_20260501_0405.json
tools/data/reports/next_exec_20260307_0342.json
tools/data/reports/ddf_20260326_0405.json
tools/data/reports/ddf_20260420_0405.json
tools/data/reports/next_exec_20260405_0330.json
tools/data/reports/ddf_20260402_0405.json
tools/data/reports/ddf_20260415_0405.json
tools/data/reports/exp_desitter_unification.json
tools/data/reports/report_20260402_0344.md
tools/data/reports/insights_20260329_0343.json
tools/data/reports/agent_20260416_0330.md
tools/data/reports/report_20260314_0342.md
tools/data/reports/ddf_20260401_0405.json
tools/data/reports/report_20260315_0342.md
tools/data/reports/ddf_20260403_0405.json
tools/data/reports/tension_projection_screening_2026-04-21.json
tools/data/reports/exp_metric_tensor_diag_long.json
tools/data/reports/falsifier_20260503_0330.json
tools/data/reports/phi_vs_silver_falsification_20260306.json
tools/data/reports/insights_20260305_0852.json
tools/data/reports/trace_map_20260305_0844.json
tools/data/reports/cycle_20260306_0342.json
tools/data/reports/ddf_20260405_0405.json
tools/data/reports/agent_20260501_0931.md
tools/data/reports/agent_20260429_1041.md
tools/data/reports/ddf_20260325_0405.json
tools/data/reports/gap_labeling_20260315_0343.json
tools/data/reports/ddf_20260404_0405.json
tools/data/reports/cycle_20260305_0844.json
tools/data/reports/report_20260403_0330.md
tools/data/reports/agent_20260417_0803.md
tools/data/reports/insights_20260403_0330.json
tools/data/reports/agent_diag2.md
tools/data/reports/next_exec_20260404_0330.json
tools/data/reports/report_20260328_0344.md
tools/data/reports/evolution_20260421_0330.md
tools/data/reports/evolution_20260425_0330.md
tools/data/reports/agent_20260418_0330.md
tools/data/reports/agent_test_0406.md
tools/data/reports/evolution_20260423_0330.md
tools/data/reports/agent_20260413_0330.md
tools/data/reports/report_20260401_0346.md
tools/data/reports/_quarantine_falsifier_29_04/agent_20260429_0852.md
tools/data/reports/_quarantine_falsifier_29_04/falsifier_20260429_0958.json
tools/data/reports/_quarantine_falsifier_29_04/falsifier_20260429_0852.json
tools/data/reports/_quarantine_falsifier_29_04/agent_20260429_0833.md
tools/data/reports/_quarantine_falsifier_29_04/agent_20260429_0958.md
tools/data/reports/_quarantine_falsifier_29_04/evolution_20260429_0852.md
tools/data/reports/_quarantine_falsifier_29_04/evolution_20260429_0958.md
tools/data/reports/_quarantine_falsifier_29_04/agent_20260429_0829.md
tools/data/reports/agent_20260502_0330.md
tools/data/reports/falsifier_20260504_0901.json
tools/data/reports/lagrangiana_20260305_1042.json
tools/data/reports/agent_20260411_0330.md
tools/data/reports/report_20260302_0341.md
tools/data/reports/agent_20260430_1946.md
tools/data/reports/agent_20260427_0330.md
tools/data/reports/ddf_20260330_0405.json
tools/data/reports/ddf_20260408_0405.json
tools/data/reports/next_exec_20260326_0343.json
tools/data/reports/next_exec_20260402_0344.json
tools/data/reports/insights_20260306_0342.json
tools/data/reports/agent_20260405_0914.md
tools/data/reports/next_exec_20260306_0342.json
tools/data/reports/exp_dR_brody_connection.json
tools/data/reports/insights_20260307_0342.json
tools/data/reports/agent_20260415_0330.md
tools/data/reports/ddf_20260421_0405.json
tools/data/reports/agent_20260501_0330.md
tools/data/reports/evolution_20260504_0330.md
tools/data/reports/agent_test_field.md
tools/data/reports/agent_20260409_0330.md
tools/data/reports/report_20260331_0345.md
tools/data/reports/ddf_20260425_0405.json
tools/data/reports/agent_20260420_0330.md
tools/data/reports/gap_labeling_20260305_0843.json
tools/data/reports/report_20260402_0756.md
tools/data/reports/agent_20260428_0330.md
tools/data/reports/agent_20260430_1919.md
tools/data/reports/fibonacci_spectrum_20260305_0844.json
tools/data/reports/cycle_20260307_0342.json
tools/data/reports/ddf_20260328_0405.json
tools/data/reports/exp_psd_prime_gaps.json
tools/data/reports/evolution_20260503_0330.md
tools/data/reports/ddf_20260407_0405.json
tools/data/reports/falsifier_20260429_1013.json
tools/data/reports/exp_brody_crossover_20260405.json
tools/data/reports/cycle_20260306_1834.json
tools/data/reports/ddf_20260323_0405.json
tools/data/reports/next_exec_20260315_0346.json
tools/data/reports/agent_20260501_0725.md
tools/data/reports/falsifier_20260430_0330.json
tools/data/reports/agent_20260405_0916.md
tools/data/reports/ddf_20260429_0405.json
tools/data/reports/ddf_20260319_0405.json
tools/data/reports/ddf_20260412_0405.json
tools/data/reports/report_20260327_0344.md
tools/data/reports/hierarchy_data.json
tools/data/reports/exp_crossover_universality.json
tools/data/reports/lagrangiana_20260305_1048.json
tools/data/reports/evolution_20260427_0330.md
tools/data/reports/cycle_20260315_0346.json
tools/data/reports/evolution_20260424_0330.md
tools/data/reports/falsifier_20260501_0931.json
tools/data/reports/agent_20260429_1013.md
tools/data/reports/exp_boundary_20260405_0825.json
tools/data/reports/next_exec_20260403_0330.json
tools/data/reports/ddf_20260406_0405.json
tools/data/reports/next_exec_20260306_1834.json
tools/data/reports/ddf_20260321_0405.json
tools/data/reports/ddf_20260418_0405.json
tools/data/reports/ddf_20260410_0405.json
tools/data/reports/agent_20260422_1616.md
tools/data/reports/next_exec_20260327_0344.json
tools/data/reports/next_exec_20260328_0344.json
tools/data/reports/incident_20260504_0721.md
tools/data/reports/ddf_20260324_0405.json
tools/data/reports/ddf_20260327_0405.json
tools/data/reports/agent_20260419_0330.md
tools/data/reports/falsifier_20260429_1041.json
tools/data/reports/gap_labeling_20260306_0341.json
tools/data/reports/agent_20260422_0330.md
tools/data/reports/evolution_20260418_0330.md
tools/data/reports/ddf_20260322_0405.json
tools/data/reports/ddf_20260424_0405.json
tools/data/reports/agent_20260504_0901.md
tools/data/reports/ddf_20260430_0405.json
tools/data/reports/insights_20260401_0346.json
tools/data/reports/report_20260405_0330.md
tools/data/reports/fibonacci_spectrum_20260306_0342.json
tools/data/reports/agent_20260407_0637.md
tools/data/reports/agent_20260408_0330.md
tools/data/reports/agent_20260406_1030.md
tools/data/reports/report_20260305_0342.md
tools/data/reports/evolution_20260422_0330.md
tools/data/reports/ddf_20260331_0405.json
tools/data/reports/gap_labeling_20260306_1834.json
tools/data/reports/agent_20260503_0330.md
tools/data/reports/exp_boundary_growth_20260405_0914.json
tools/data/reports/ddf_20260417_0405.json
tools/data/reports/report_20260330_0344.md
tools/data/reports/falsifier_20260501_0725.json
tools/data/reports/agent_20260405_0825.md
tools/data/reports/agent_20260406_0714.md
tools/data/reports/report_20260315_0801.md
tools/data/reports/insights_20260306_1834.json
tools/data/reports/ddf_20260422_0405.json
tools/data/reports/ddf_20260318_0405.json
tools/data/reports/next_exec_20260401_0346.json
tools/data/reports/ddf_20260427_0405.json
tools/data/reports/next_exec_20260331_0345.json
tools/data/reports/report_20260404_1852.md
tools/data/reports/exp_excess_scaling_20260405.json
tools/data/reports/exp_ricci_primes.json
tools/data/reports/report_20260405_0715.md
tools/data/reports/falsifier_20260430_1946.json
tools/data/reports/ddf_20260426_0405.json
tools/data/reports/report_20260307_0342.md
tools/data/reports/report_20260303_0341.md
tools/data/experiment_results.json
tools/data/reports/next_exec_20260305_1111.json
tools/data/reports/evolution_20260422_1616.md
tools/data/reports/falsifier_20260501_0330.json
tools/data/reports/agent_20260426_0330.md
tools/data/reports/exp_number_variance_test.json
tools/data/reports/mapping_validation_2026-04-21.json
tools/data/reports/ddf_20260329_0405.json
tools/data/reports/report_20260305_2121.md
tools/data/reports/ddf_20260317_0405.json
tools/data/reports/ddf_20260315_0815.json
tools/data/reports/ddf_20260419_0405.json
tools/data/reports/exp_geodesic_deviation_primes.json
tools/data/reports/agent_20260419_0755.md
tools/data/reports/agent_20260412_0330.md
tools/data/reports/ddf_20260414_0405.json
tools/data/reports/evolution_20260417_0803.md
tools/data/reports/insights_20260405_0729.json
tools/data/reports/report_20260304_0342.md
tools/data/reports/agent_20260425_0330.md
tools/data/reports/ddf_20260409_0405.json
tools/data/reports/ddf_20260411_0405.json
tools/data/reports/exp_acf_decay_data.json
tools/data/reports/falsifier_20260502_0330.json
tools/data/reports/fibonacci_spectrum_20260305_0756.json
tools/data/reports/trace_map_20260305_0843.json
tools/data/reports/report_20260404_0330.md
tools/data/reports/next_exec_20260314_0342.json
tools/data/reports/report_20260306_0341.md
tools/data/reports/ddf_20260502_0405.json
tools/data/reports/ddf_20260503_0405.json
tools/data/reports/next_exec_20260330_0344.json
tools/data/reports/ddf_20260315_0811.json
tools/data/reports/ddf_20260320_0405.json
tools/data/reports/report_20260331_1809.md
tools/data/reports/ddf_20260316_0405.json
tools/data/reports/falsifier_20260430_1919.json
tools/data/reports/ddf_20260504_0405.json
tools/data/reports/ddf_20260413_0405.json
tools/data/reports/ddf_20260416_0405.json
tools/data/reports/evolution_20260419_0330.md
tools/data/reports/report_20260326_0343.md
tools/data/reports/evolution_20260428_0330.md
tools/data/reports/falsifier_20260430_1905.json
tools/data/reports/agent_20260405_0919.md
tools/data/reports/ddf_20260423_0405.json

codex
Il campo punta esplicitamente a META: il rischio è che i test recenti stiano misurando tautologie o accoppiamenti già imposti dal setup. Ora verifico la proiezione aggiornata e leggo gli strumenti/report più vicini per scegliere una misura che non replichi gli ultimi tre run.
exec
/bin/bash -lc "sed -n '1,260p' tools/exp_meta_tautology_test.py" in /opt/MM_D-ND
 succeeded in 0ms:
#!/usr/bin/env python3
"""META Tautology Discriminator — which observables measure real structure vs density artifacts?

Four observables x three synthetic controls. Each observable gets a verdict:
STRUCTURAL (distinguishes primes from all synthetics) or TAUTOLOGICAL (fails on at least one).

Observables:
  1. r-statistic (spacing ratio) — short-range repulsion
  2. Mod-3 ordering fraction — the algebraic memory channel
  3. Lag-1 autocorrelation — sequential gap correlation
  4. Two-channel delta-r — magnitude vs ordering decomposition

Controls:
  A. Shuffled gaps — same distribution, destroyed order
  B. Cramer model — PNT density, independent exponential gaps
  C. Hardy-Littlewood model — correct pair correlations, no higher-order structure
"""

import argparse
import json
import numpy as np
from sympy import primerange


def get_primes(n_max):
    return np.array(list(primerange(2, n_max + 1)), dtype=np.int64)


def r_statistic(gaps):
    """Mean spacing ratio min(s_i, s_{i+1}) / max(s_i, s_{i+1})."""
    s1 = gaps[:-1]
    s2 = gaps[1:]
    mn = np.minimum(s1, s2)
    mx = np.maximum(s1, s2)
    mask = mx > 0
    return np.mean(mn[mask] / mx[mask])


def mod3_ordering_fraction(gaps):
    """Fraction of consecutive gap pairs where mod-3 class is preserved."""
    classes = gaps % 3
    same = np.sum(classes[:-1] == classes[1:])
    return same / len(classes[:-1])


def lag1_autocorrelation(gaps):
    """Pearson autocorrelation at lag 1."""
    g = gaps.astype(np.float64)
    g = g - g.mean()
    if g.std() == 0:
        return 0.0
    return np.corrcoef(g[:-1], g[1:])[0, 1]


def two_channel_delta_r(gaps):
    """Decompose into magnitude and ordering channels, return delta-r for each."""
    g = gaps.astype(np.float64)
    median_g = np.median(g)
    binary = (g > median_g).astype(np.float64)  # ordering channel
    magnitude = np.abs(g - median_g)              # magnitude channel

    r_ord = r_statistic_from_signal(binary)
    r_mag = r_statistic_from_signal(magnitude)

    # Shuffle baseline
    rng = np.random.default_rng(42)
    r_ord_shuf = []
    r_mag_shuf = []
    for _ in range(20):
        idx = rng.permutation(len(gaps))
        b_s = binary[idx]
        m_s = magnitude[idx]
        r_ord_shuf.append(r_statistic_from_signal(b_s))
        r_mag_shuf.append(r_statistic_from_signal(m_s))

    dr_ord = (r_ord - np.mean(r_ord_shuf)) / (np.std(r_ord_shuf) + 1e-12)
    dr_mag = (r_mag - np.mean(r_mag_shuf)) / (np.std(r_mag_shuf) + 1e-12)
    return dr_ord, dr_mag


def r_statistic_from_signal(sig):
    """r-statistic on arbitrary positive signal (add offset if needed)."""
    s = sig - sig.min() + 1e-6
    s1 = s[:-1]
    s2 = s[1:]
    mn = np.minimum(s1, s2)
    mx = np.maximum(s1, s2)
    mask = mx > 0
    return np.mean(mn[mask] / mx[mask])


# === Synthetic generators ===

def shuffled_gaps(gaps, rng):
    """Same gap distribution, destroyed sequential order."""
    g = gaps.copy()
    rng.shuffle(g)
    return g


def cramer_random_gaps(n_gaps, mean_gap, rng):
    """Independent exponential gaps rounded to even (like PNT density)."""
    raw = rng.exponential(mean_gap, size=n_gaps)
    g = np.round(raw / 2) * 2
    g = np.maximum(g, 2).astype(np.int64)
    return g


def hardy_littlewood_gaps(gaps_real, rng):
    """Markov(1) model matching lag-1 autocorrelation of real primes.
    Preserves pair correlation structure but not higher-order."""
    g = gaps_real.astype(np.float64)
    mean_g = g.mean()
    std_g = g.std()
    rho = np.corrcoef(g[:-1], g[1:])[0, 1]

    # AR(1) process with correct mean, std, lag-1
    n = len(gaps_real)
    result = np.zeros(n)
    result[0] = mean_g
    noise_std = std_g * np.sqrt(1 - rho**2)
    for i in range(1, n):
        result[i] = mean_g + rho * (result[i-1] - mean_g) + rng.normal(0, noise_std)

    # Round to even, clip to >= 2
    result = np.round(result / 2) * 2
    result = np.maximum(result, 2).astype(np.int64)
    return result


def run(n_primes_max=600000, n_trials=20):
    """Run the META tautology test."""
    print(f"Generating primes up to {n_primes_max}...")
    primes = get_primes(n_primes_max)
    gaps = np.diff(primes)

    # Use a window in the middle to avoid small-prime effects
    N = min(len(gaps), 50000)
    start = len(gaps) // 4
    gaps_window = gaps[start:start + N]
    mean_gap = float(gaps_window.mean())

    print(f"Using {N} gaps starting at index {start} (mean gap = {mean_gap:.2f})")

    # Real primes observables
    print("\n=== REAL PRIMES ===")
    real_r = r_statistic(gaps_window)
    real_mod3 = mod3_ordering_fraction(gaps_window)
    real_lag1 = lag1_autocorrelation(gaps_window)
    real_dr_ord, real_dr_mag = two_channel_delta_r(gaps_window)
    print(f"  r-stat:     {real_r:.6f}")
    print(f"  mod3-frac:  {real_mod3:.6f}")
    print(f"  lag1-acf:   {real_lag1:.6f}")
    print(f"  dr-ord:     {real_dr_ord:.2f}σ")
    print(f"  dr-mag:     {real_dr_mag:.2f}σ")

    # Run synthetics
    rng = np.random.default_rng(2026)

    results = {
        'real': {
            'r_stat': real_r, 'mod3': real_mod3, 'lag1': real_lag1,
            'dr_ord': real_dr_ord, 'dr_mag': real_dr_mag
        }
    }

    for name, generator in [
        ('shuffled', lambda rng_: shuffled_gaps(gaps_window, rng_)),
        ('cramer', lambda rng_: cramer_random_gaps(N, mean_gap, rng_)),
        ('HL_markov', lambda rng_: hardy_littlewood_gaps(gaps_window, rng_)),
    ]:
        print(f"\n=== {name.upper()} (n_trials={n_trials}) ===")
        obs = {'r_stat': [], 'mod3': [], 'lag1': [], 'dr_ord': [], 'dr_mag': []}

        for t in range(n_trials):
            trial_rng = np.random.default_rng(rng.integers(0, 2**31))
            syn_gaps = generator(trial_rng)
            obs['r_stat'].append(r_statistic(syn_gaps))
            obs['mod3'].append(mod3_ordering_fraction(syn_gaps))
            obs['lag1'].append(lag1_autocorrelation(syn_gaps))
            dr_o, dr_m = two_channel_delta_r(syn_gaps)
            obs['dr_ord'].append(dr_o)
            obs['dr_mag'].append(dr_m)

        results[name] = {}
        for key in obs:
            arr = np.array(obs[key])
            results[name][key] = {
                'mean': float(np.mean(arr)),
                'std': float(np.std(arr)),
                'min': float(np.min(arr)),
                'max': float(np.max(arr))
            }
            print(f"  {key:12s}: {np.mean(arr):+.6f} ± {np.std(arr):.6f}")

    # === Discrimination verdict ===
    print("\n" + "=" * 60)
    print("DISCRIMINATION VERDICTS")
    print("=" * 60)

    verdicts = {}
    for obs_name in ['r_stat', 'mod3', 'lag1', 'dr_ord', 'dr_mag']:
        real_val = results['real'][obs_name]
        discriminates = {}
        for syn_name in ['shuffled', 'cramer', 'HL_markov']:
            syn = results[syn_name][obs_name]
            z = (real_val - syn['mean']) / (syn['std'] + 1e-12)
            discriminates[syn_name] = abs(z)
            label = "YES" if abs(z) > 3.0 else "no"
            print(f"  {obs_name:12s} vs {syn_name:12s}: z = {z:+8.2f}  [{label}]")

        all_pass = all(v > 3.0 for v in discriminates.values())
        verdict = "STRUCTURAL" if all_pass else "TAUTOLOGICAL"
        verdicts[obs_name] = {
            'verdict': verdict,
            'z_scores': {k: float(v) for k, v in discriminates.items()}
        }
        print(f"  → {obs_name}: {verdict}")
        print()

    # Save results
    output = {
        'n_gaps': N,
        'start_index': start,
        'mean_gap': mean_gap,
        'n_trials': n_trials,
        'real': results['real'],
        'synthetics': {k: results[k] for k in ['shuffled', 'cramer', 'HL_markov']},
        'verdicts': verdicts
    }
    out_path = 'tools/data/meta_tautology_test.json'
    with open(out_path, 'w') as f:
        json.dump(output, f, indent=2)
    print(f"\nResults saved to {out_path}")

    return output


if __name__ == '__main__':
    parser = argparse.ArgumentParser(description='META Tautology Discriminator')
    parser.add_argument('--n-primes', type=int, default=600000, help='Prime sieve limit')
    parser.add_argument('--n-trials', type=int, default=20, help='Synthetic trials per model')
    args = parser.parse_args()
    run(n_primes_max=args.n_primes, n_trials=args.n_trials)

exec
/bin/bash -lc "sed -n '1,220p' tools/data/reports/agent_20260504_0901.md" in /opt/MM_D-ND
 succeeded in 0ms:
# Agent Report — The Two Markov Layers Are Coupled at the Boundary: One Phase Transition, Two Projections

**Date**: 2026-05-04 09:01
**Piano**: 60
**Tension explored**: META (0.5) + BOUNDARY (0.8) + DIPOLAR_ORDERING (0.8)

## Claim Under Test

> The previous run found two orthogonal visible layers in prime gap memory: Layer 1 (pairs: SR, L1) and Layer 2 (triples: SR2, triple_var). The crossover under partial shuffle shows a phase transition in the (SR, L1) plane. META question: does Layer 2 transition at a different critical alpha than Layer 1? If yes, the boundary has genuine 3D depth. If no, the two layers are projections of a single phase transition.

## Question

Do the pair-statistics layer (SR, L1) and triple-statistics layer (SR2, triple_var) undergo independent transitions when ordering is destroyed by partial shuffle, or are they coupled?

## Experiment Design

- **Method**: Partial shuffle crossover with 20 alpha steps (0.05 to 0.95), 30 trials per step.
- **Observables**: 4 total — Layer 1: spacing ratio (SR), lag-1 ACF (L1). Layer 2: next-nearest-neighbor spacing ratio (SR2), normalized triple variance (triple_var).
- **Metric**: Retention = (value - baseline) / (original - baseline). Critical alpha = alpha where retention drops below 0.50.
- **Null baseline**: Full shuffle (alpha=1.0, 90 trials) for each sequence.
- **Sequences**: Prime gaps (N=50000), GUE eigenvalue gaps (200x200 matrices, 250 realizations), Poisson iid exponential gaps (N=50000).
- **Poisson control**: If Poisson shows layer separation, the metric is noise-sensitive. If only structured sequences show coupling, the coupling is real.

## Results

### Critical alpha (50% retention)

| Sequence | SR (L1) | L1 (L1) | SR2 (L2) | triple_var (L2) | L1 mean | L2 mean | Delta |
|----------|---------|---------|----------|-----------------|---------|---------|-------|
| Primes   | 0.334   | 0.334   | 0.334    | 0.334           | 0.334   | 0.334   | +0.000 |
| GUE      | 0.334   | 0.287   | 0.334    | 0.334           | 0.311   | 0.334   | +0.024 |
| Poisson  | 0.239   | 0.334   | 0.097    | 0.097           | 0.287   | 0.097   | -0.189 |

### Retention at alpha = 0.33 (near critical)

| Sequence | L1 avg | L2 avg | Difference |
|----------|--------|--------|------------|
| Primes   | 0.435  | 0.446  | -0.011     |
| GUE      | 0.444  | 0.445  | -0.001     |
| Poisson  | 0.330  | 0.591  | -0.261     |

### Poisson null verification

All Poisson original-vs-baseline z-scores are < 2 (SR: z=0.91, L1: z=-1.47, SR2: z=-0.17, triple_var: z=-0.56). The Poisson "signal" is noise. The apparent layer separation (Delta = -0.189) in Poisson is an artifact: when the signal-to-noise is < 2, the retention metric amplifies noise differently for each observable.

### Zero-crossing alpha (sign flip)

| Sequence | SR     | L1     | SR2    | triple_var |
|----------|--------|--------|--------|------------|
| Primes   | 0.917  | 0.846  | 0.914  | 0.903      |
| GUE      | 0.895  | None   | 0.902  | None       |

For primes, all observables flip sign at alpha > 0.84. L1 flips earliest (0.846), SR latest (0.917). The ordering is L1 < triple_var < SR2 < SR — interleaved between layers, not grouped by layer.

## Key Findings

1. **The two Markov layers are coupled at the boundary.** For primes, the critical alpha is identical across all 4 observables (0.334). For GUE, the difference is 0.024 (within the alpha step resolution of 0.047). The partial shuffle destroys pair-statistics and triple-statistics at the same rate. The boundary is a single phase transition, not two independent ones.

2. **The coupling is specific to structured sequences.** Poisson (iid, no ordering) shows Delta = -0.189 — spurious separation from noise amplification. Primes (Delta = 0.000) and GUE (Delta = 0.024) show coupling. This rules out the coupling being a trivial property of the metric.

3. **The two-layer decomposition is a decomposition of observables, not of the ordering.** The previous run correctly identified that SR and L1 are sensitive to Markov-1 (pair) statistics while SR2 and triple_var are sensitive to Markov-2 (triple) statistics. But when the ordering is destroyed uniformly (partial shuffle), both layers lose signal at the same rate. The layers are different projections of one ordering, not independent degrees of freedom.

4. **The zero-crossing order is interleaved, not grouped by layer.** For primes: L1(0.846) < triple_var(0.903) < SR2(0.914) < SR(0.917). If layers were independent, we'd expect L1 grouping with SR and SR2 grouping with triple_var. The interleaving confirms coupling.

## 5-Lens Self-Check

- **L1 (hard constraint vs bias)**: No absolute claims. "Coupled" is quantified as |Delta| < 0.05 (resolution limit). The data shows 0.000 and 0.024.
- **L2 (absolute vs ratio)**: Retention is already normalized. All comparisons are in the same units (fraction of original signal).
- **L3 (no silent patching)**: This does NOT falsify the two-layer finding. The layers remain real as a decomposition of Markov order sensitivity. What's constrained is their independence at the boundary.
- **L4 (edge case)**: Poisson separation is explicitly identified as noise artifact with z-score evidence.
- **L5 (re-discovery)**: That partial shuffle destroys correlations uniformly regardless of order is consistent with the known property that random permutations break all multi-point correlations simultaneously (not order-by-order). The specific quantification on prime gaps and GUE with the Markov-layer framework is new to this lab, but the underlying principle is not novel. Tagged as CONSTRAINT, not NEW.

## Verdict

**CONSTRAINT on BOUNDARY + DIPOLAR_ORDERING**: The two Markov layers (pairs → plane, triples → depth) are coupled at the partial-shuffle boundary. The boundary is a single phase transition with one critical alpha (~0.33 for 50% retention). The two-layer decomposition describes WHAT is measured (which observables are sensitive to which Markov order), not HOW the ordering is destroyed. The boundary remains 2D in the shuffle parameter.

**Consecutio**: Since the layers are coupled under uniform shuffle, the question becomes: is there a NON-uniform perturbation that decouples them? Specifically, a perturbation that destroys pair correlations but preserves triple correlations (or vice versa). If such a perturbation exists, the layers are genuinely independent degrees of freedom that happen to be coupled under this particular destruction method. If not, the two-layer decomposition is a spectral decomposition of a single ordering dimension.

## Bicono della scoperta

- **Due radici** (dipolo primario): Layer 1 (pairs, near-neighbor) and Layer 2 (triples, next-nearest-neighbor) — two ways the same ordering manifests in different correlation windows
- **Singolare**: The ordering itself — before decomposition into pair and triple statistics. It exists as one thing; the layers are the observer's choice of measurement window.
- **Invariante di passaggio**: The critical alpha (0.334) survives across layers and across structured sequences (primes and GUE). The boundary location is invariant to which layer you observe through.
- **Campo di possibilita**: Possible — search for selective perturbations that decouple the layers (pair-preserving shuffle, triple-preserving shuffle). Not possible — claim the boundary has independent depth dimensions from partial shuffle alone.

## Files

- Script: `tools/exp_3d_boundary_layers.py` (reusable, parameterized)
- Data: `tools/data/3d_boundary_layers.json`
- Report: `tools/data/reports/agent_20260504_0901.md`

exec
/bin/bash -lc "sed -n '1,220p' tools/exp_two_layer_universality.py" in /opt/MM_D-ND
 succeeded in 0ms:
#!/usr/bin/env python3
"""
exp_two_layer_universality.py — Is the two-layer Markov memory structure universal?

Consecutio from piano 60g (agent_20260503_0330):
  Prime gap memory decomposes into two orthogonal layers:
    Layer 1 (pairs, Mk1): SR, L1
    Layer 2 (triples, Mk2): SR2, L2, cond_entropy, triple_var, num_var_10
  Question: Is this decomposition a property of ALL ordered sequences,
  or specific to primes?

Method:
  1. Generate 7 gap/spacing sequences: primes, GUE, Poisson, AR(1), logistic,
     periodic, Fibonacci gaps
  2. For each, build Markov-k surrogates (k=0,1,2,3)
  3. Compute 8 observables, measure z-scores vs each Mk level
  4. Classify which observables are "captured" by Mk1 (|z_Mk0|>>2, |z_Mk1|<2)
     vs Mk2 (|z_Mk1|>>2, |z_Mk2|<2) vs Mk3
  5. Compare layer assignments across sequences

Null hypothesis: If universal, all sequences assign SR,L1 to Layer1 and SR2,L2 to Layer2.
Alternative: Layer assignment is sequence-specific — the structure is diagnostic.

Usage:
    python tools/exp_two_layer_universality.py [--N 100000] [--n_surr 30]
"""

import argparse
import json
import numpy as np
from pathlib import Path
from collections import Counter


def get_primes(n_max):
    sieve = np.ones(n_max + 1, dtype=bool)
    sieve[0] = sieve[1] = False
    for i in range(2, int(n_max**0.5) + 1):
        if sieve[i]:
            sieve[i*i::i] = False
    return np.where(sieve)[0]


def build_markov_chain(gaps, order, n_bins=12):
    percentiles = np.linspace(0, 100, n_bins + 1)
    edges = np.percentile(gaps, percentiles)
    edges[0] = gaps.min() - 0.5
    edges[-1] = gaps.max() + 0.5
    binned = np.digitize(gaps, edges) - 1
    binned = np.clip(binned, 0, n_bins - 1)
    gap_pools = {}
    for b, g in zip(binned, gaps):
        gap_pools.setdefault(b, []).append(g)
    trans = {}
    for i in range(len(binned) - order):
        state = tuple(binned[i:i + order])
        nxt = binned[i + order]
        if state not in trans:
            trans[state] = Counter()
        trans[state][nxt] += 1
    trans_prob = {}
    for state, counts in trans.items():
        total = sum(counts.values())
        trans_prob[state] = {k: v / total for k, v in counts.items()}
    return binned, edges, gap_pools, trans_prob


def generate_markov_surrogate(gaps, order, n_bins=12, rng=None):
    if rng is None:
        rng = np.random.default_rng()
    binned, edges, gap_pools, trans_prob = build_markov_chain(gaps, order, n_bins)
    n = len(gaps)
    result = np.zeros(n)
    start_idx = rng.integers(0, len(binned) - order)
    state = tuple(binned[start_idx:start_idx + order])
    for j in range(order):
        pool = gap_pools[state[j]]
        result[j] = pool[rng.integers(0, len(pool))]
    for i in range(order, n):
        if state in trans_prob:
            probs = trans_prob[state]
            bins_avail = list(probs.keys())
            p = np.array([probs[b] for b in bins_avail])
            nxt_bin = bins_avail[rng.choice(len(bins_avail), p=p)]
        else:
            nxt_bin = rng.integers(0, n_bins)
        pool = gap_pools.get(nxt_bin, gap_pools[list(gap_pools.keys())[0]])
        result[i] = pool[rng.integers(0, len(pool))]
        state = (*state[1:], nxt_bin)
    return result


# --- Observables ---
def spacing_ratio(gaps):
    s = gaps[:-1]
    s1 = gaps[1:]
    r = np.minimum(s, s1) / np.maximum(s, s1)
    return np.mean(r[np.isfinite(r)])

def lag_k_acf(gaps, k=1):
    g = gaps - np.mean(gaps)
    if np.var(gaps) == 0:
        return 0.0
    n = len(g)
    return np.sum(g[:n-k] * g[k:]) / np.sum(g**2)

def next_nearest_sr(gaps):
    """SR2: spacing ratio of next-nearest-neighbor gaps (skip one)."""
    if len(gaps) < 3:
        return 0.5
    s = gaps[:-2]
    s2 = gaps[2:]
    r = np.minimum(s, s2) / np.maximum(s, s2)
    return np.mean(r[np.isfinite(r)])

def cond_entropy_l2(gaps, n_bins=12):
    """Conditional entropy H(g_{n+2} | g_n, g_{n+1})."""
    if len(gaps) < 3:
        return 0.0
    percentiles = np.linspace(0, 100, n_bins + 1)
    edges = np.percentile(gaps, percentiles)
    edges[0] = gaps.min() - 0.5
    edges[-1] = gaps.max() + 0.5
    binned = np.digitize(gaps, edges) - 1
    binned = np.clip(binned, 0, n_bins - 1)
    joint = Counter()
    cond = Counter()
    for i in range(len(binned) - 2):
        state = (binned[i], binned[i+1])
        nxt = binned[i+2]
        joint[(state, nxt)] += 1
        cond[state] += 1
    h = 0.0
    for (state, nxt), cnt in joint.items():
        p = cnt / cond[state]
        if p > 0:
            h -= (cnt / (len(binned) - 2)) * np.log2(p)
    return h

def triple_var(gaps):
    """Variance of (g_n, g_{n+1}, g_{n+2}) triple sums."""
    if len(gaps) < 3:
        return 0.0
    t = gaps[:-2] + gaps[1:-1] + gaps[2:]
    return np.var(t)

def num_var_window(gaps, w=10):
    """Number variance in windows of size w."""
    if len(gaps) < w:
        return np.var(gaps)
    counts = np.array([np.sum(gaps[i:i+w]) for i in range(len(gaps) - w)])
    return np.var(counts)

def run_length_mean(gaps):
    """Mean run length (consecutive increases or decreases)."""
    diffs = np.diff(gaps)
    signs = np.sign(diffs)
    runs = []
    current = 1
    for i in range(1, len(signs)):
        if signs[i] == signs[i-1] and signs[i] != 0:
            current += 1
        else:
            runs.append(current)
            current = 1
    runs.append(current)
    return np.mean(runs)


OBSERVABLES = {
    'SR': spacing_ratio,
    'L1': lambda g: lag_k_acf(g, 1),
    'L2': lambda g: lag_k_acf(g, 2),
    'SR2': next_nearest_sr,
    'cond_entropy': cond_entropy_l2,
    'triple_var': triple_var,
    'num_var_10': num_var_window,
    'run_length': run_length_mean,
}


# --- Sequence generators ---
def gen_prime_gaps(N):
    primes = get_primes(int(N * 15))
    primes = primes[:N+1]
    return np.diff(primes).astype(float)

def gen_gue_spacings(N, rng=None):
    if rng is None:
        rng = np.random.default_rng()
    dim = min(int(np.sqrt(2 * N)) + 50, 1500)
    H = rng.standard_normal((dim, dim))
    H = (H + H.T) / 2
    eigs = np.sort(np.linalg.eigvalsh(H))
    spacings = np.diff(eigs)
    spacings = spacings[spacings > 0]
    # Unfold
    mean_s = np.mean(spacings)
    spacings = spacings / mean_s
    if len(spacings) > N:
        spacings = spacings[:N]
    return spacings

def gen_poisson_spacings(N, rng=None):
    if rng is None:
        rng = np.random.default_rng()
    return rng.exponential(1.0, N)

def gen_ar1_gaps(N, phi=0.5, rng=None):
    if rng is None:
        rng = np.random.default_rng()
    x = np.zeros(N)
    x[0] = rng.standard_normal()
    for i in range(1, N):
        x[i] = phi * x[i-1] + rng.standard_normal()
    # Shift to positive
    x = x - x.min() + 1.0
    return x

def gen_logistic_gaps(N, r=3.95):

exec
/bin/bash -lc 'python tools/dnd_scenario.py --best' in /opt/MM_D-ND
 succeeded in 103ms:
codex
La consecutio del report 09:01 apre una domanda precisa: se lo shuffle uniforme accoppia layer pair/triple, una perturbazione selettiva riesce a separarli? Misuro questo come test META: se nessuna perturbazione selettiva separa i layer, il “3D boundary” era solo proiezione osservativa; se li separa, c’è un grado operativo in più.
exec
/bin/bash -lc "sed -n '1,280p' tools/exp_3d_boundary_layers.py" in /opt/MM_D-ND
 succeeded in 0ms:
#!/usr/bin/env python3
"""
exp_3d_boundary_layers.py — Does the boundary have 3D structure?

Consecutio from piano 60 runs:
  - Layer 1 (pairs, Mk1) → SR, L1 (the dipolar plane)
  - Layer 2 (triples, Mk2) → SR2, triple_var (depth)
  - The crossover (partial shuffle) shows a phase transition in (SR, L1)
  - META question: does Layer 2 transition at the SAME critical alpha?

If same α_c → the boundary is 2D (Layer 2 follows Layer 1 = partial tautology)
If different α_c → the boundary has genuine 3D depth (two independent transitions)

Tests on: primes, GUE, Poisson baseline.

Usage:
    python tools/exp_3d_boundary_layers.py [--N 50000] [--n_alpha 20] [--n_trials 30]
"""

import argparse
import json
import numpy as np
from scipy import stats
from pathlib import Path


def get_primes(n_max):
    sieve = np.ones(n_max + 1, dtype=bool)
    sieve[0] = sieve[1] = False
    for i in range(2, int(n_max**0.5) + 1):
        if sieve[i]:
            sieve[i*i::i] = False
    return np.where(sieve)[0]


def gue_gaps(N_mat, n_matrices, rng):
    """Generate GUE eigenvalue gaps."""
    all_gaps = []
    for _ in range(n_matrices):
        H = rng.standard_normal((N_mat, N_mat)) + 1j * rng.standard_normal((N_mat, N_mat))
        H = (H + H.conj().T) / 2
        evals = np.sort(np.linalg.eigvalsh(H))
        gaps = np.diff(evals)
        gaps = gaps[gaps > 0]
        all_gaps.extend(gaps.tolist())
    return np.array(all_gaps)


def partial_shuffle(seq, alpha, rng):
    s = seq.copy()
    n = len(s)
    k = int(alpha * n)
    if k < 2:
        return s
    idx = rng.choice(n, size=k, replace=False)
    vals = s[idx].copy()
    rng.shuffle(vals)
    s[idx] = vals
    return s


# --- Layer 1 observables (pair statistics) ---
def obs_spacing_ratio(gaps):
    r = np.minimum(gaps[:-1], gaps[1:]) / np.maximum(gaps[:-1], gaps[1:])
    return np.mean(r[np.isfinite(r)])

def obs_lag1_acf(gaps):
    g = gaps - np.mean(gaps)
    c0 = np.mean(g**2)
    if c0 == 0: return 0.0
    return np.mean(g[:-1] * g[1:]) / c0

# --- Layer 2 observables (triple statistics) ---
def obs_sr2(gaps):
    """Next-nearest-neighbor spacing ratio: min(g_n, g_{n+2})/max(g_n, g_{n+2})"""
    r = np.minimum(gaps[:-2], gaps[2:]) / np.maximum(gaps[:-2], gaps[2:])
    return np.mean(r[np.isfinite(r)])

def obs_triple_var(gaps):
    """Variance of consecutive triple sums, normalized."""
    triples = gaps[:-2] + gaps[1:-1] + gaps[2:]
    v = np.var(gaps)
    if v == 0: return 0.0
    return np.var(triples) / v


def run_crossover(gaps, alphas, n_trials, rng, label=""):
    """Compute all 4 observables at each alpha level."""
    obs_fns = {
        'SR': obs_spacing_ratio,
        'L1': obs_lag1_acf,
        'SR2': obs_sr2,
        'triple_var': obs_triple_var,
    }

    # Full shuffle baseline (alpha=1.0)
    baselines = {name: [] for name in obs_fns}
    for _ in range(n_trials * 3):
        shuffled = partial_shuffle(gaps, 1.0, rng)
        for name, fn in obs_fns.items():
            baselines[name].append(fn(shuffled))
    baseline_mean = {name: np.mean(vals) for name, vals in baselines.items()}
    baseline_std = {name: np.std(vals) for name, vals in baselines.items()}

    # Original values (alpha=0)
    originals = {name: fn(gaps) for name, fn in obs_fns.items()}

    results = []
    for alpha in alphas:
        trial_vals = {name: [] for name in obs_fns}
        for _ in range(n_trials):
            s = partial_shuffle(gaps, alpha, rng)
            for name, fn in obs_fns.items():
                trial_vals[name].append(fn(s))

        row = {'alpha': float(alpha)}
        for name in obs_fns:
            mean_val = np.mean(trial_vals[name])
            std_val = np.std(trial_vals[name])
            # Fraction of original signal retained
            orig_delta = originals[name] - baseline_mean[name]
            curr_delta = mean_val - baseline_mean[name]
            if abs(orig_delta) > 1e-12:
                retention = curr_delta / orig_delta
            else:
                retention = 0.0
            row[f'{name}_mean'] = float(mean_val)
            row[f'{name}_std'] = float(std_val)
            row[f'{name}_retention'] = float(retention)
        results.append(row)

    return results, originals, baseline_mean, baseline_std


def find_critical_alpha(results, obs_name, threshold=0.5):
    """Find alpha where retention drops below threshold (signal half-life)."""
    for r in results:
        if r[f'{obs_name}_retention'] < threshold:
            return r['alpha']
    return 1.0  # never crossed


def find_zero_crossing(results, obs_name):
    """Find alpha where retention crosses zero (sign flip)."""
    for i in range(1, len(results)):
        r0 = results[i-1][f'{obs_name}_retention']
        r1 = results[i][f'{obs_name}_retention']
        if r0 * r1 < 0:
            # Linear interpolation
            a0 = results[i-1]['alpha']
            a1 = results[i]['alpha']
            alpha_zero = a0 + (a1 - a0) * abs(r0) / (abs(r0) + abs(r1))
            return alpha_zero
    return None


def main():
    parser = argparse.ArgumentParser()
    parser.add_argument('--N', type=int, default=50000, help='Number of primes')
    parser.add_argument('--n_alpha', type=int, default=20, help='Number of alpha steps')
    parser.add_argument('--n_trials', type=int, default=30, help='Trials per alpha')
    parser.add_argument('--seed', type=int, default=42)
    args = parser.parse_args()

    rng = np.random.default_rng(args.seed)
    alphas = np.linspace(0.05, 0.95, args.n_alpha)

    print(f"=== 3D Boundary Layer Experiment ===")
    print(f"N_primes={args.N}, n_alpha={args.n_alpha}, n_trials={args.n_trials}")

    # --- Primes ---
    print("\n--- PRIMES ---")
    primes = get_primes(args.N * 20)[:args.N]
    prime_gaps = np.diff(primes).astype(float)
    prime_results, prime_orig, prime_bl_mean, prime_bl_std = run_crossover(
        prime_gaps, alphas, args.n_trials, rng, "Primes"
    )

    # --- GUE ---
    print("\n--- GUE ---")
    n_mat = 200
    n_matrices = max(5, args.N // n_mat)
    gue_g = gue_gaps(n_mat, n_matrices, rng)
    if len(gue_g) > args.N:
        gue_g = gue_g[:args.N]
    gue_results, gue_orig, gue_bl_mean, gue_bl_std = run_crossover(
        gue_g, alphas, args.n_trials, rng, "GUE"
    )

    # --- Poisson (exponential gaps, iid) ---
    print("\n--- POISSON ---")
    poisson_gaps = rng.exponential(1.0, size=args.N)
    pois_results, pois_orig, pois_bl_mean, pois_bl_std = run_crossover(
        poisson_gaps, alphas, args.n_trials, rng, "Poisson"
    )

    # --- Analysis ---
    obs_names = ['SR', 'L1', 'SR2', 'triple_var']
    layer_map = {'SR': 'L1_pair', 'L1': 'L1_pair', 'SR2': 'L2_triple', 'triple_var': 'L2_triple'}

    output = {
        'experiment': '3D Boundary Layers',
        'question': 'Do Layer 1 (pairs) and Layer 2 (triples) transition at the same critical alpha?',
        'params': {'N': args.N, 'n_alpha': args.n_alpha, 'n_trials': args.n_trials, 'seed': args.seed},
        'sequences': {}
    }

    for name, results, originals, bl_mean, bl_std in [
        ('primes', prime_results, prime_orig, prime_bl_mean, prime_bl_std),
        ('gue', gue_results, gue_orig, gue_bl_mean, gue_bl_std),
        ('poisson', pois_results, pois_orig, pois_bl_mean, pois_bl_std),
    ]:
        seq_data = {
            'originals': {k: float(v) for k, v in originals.items()},
            'baseline_mean': {k: float(v) for k, v in bl_mean.items()},
            'baseline_std': {k: float(v) for k, v in bl_std.items()},
            'critical_alpha_50': {},
            'zero_crossing': {},
            'retention_curve': results,
        }
        print(f"\n=== {name.upper()} ===")
        print(f"{'Observable':>12} {'Layer':>10} {'Original':>10} {'Baseline':>10} {'α_crit(50%)':>12} {'α_zero':>8}")
        for obs in obs_names:
            ac = find_critical_alpha(results, obs)
            az = find_zero_crossing(results, obs)
            seq_data['critical_alpha_50'][obs] = float(ac)
            seq_data['zero_crossing'][obs] = float(az) if az else None
            print(f"{obs:>12} {layer_map[obs]:>10} {originals[obs]:>10.5f} {bl_mean[obs]:>10.5f} {ac:>12.3f} {str(az and f'{az:.3f}') or 'none':>8}")

        # Layer separation: difference in critical alpha between layers
        l1_crit = np.mean([seq_data['critical_alpha_50'][o] for o in ['SR', 'L1']])
        l2_crit = np.mean([seq_data['critical_alpha_50'][o] for o in ['SR2', 'triple_var']])
        delta_crit = l2_crit - l1_crit
        seq_data['layer_separation'] = {
            'L1_mean_crit': float(l1_crit),
            'L2_mean_crit': float(l2_crit),
            'delta': float(delta_crit),
        }
        print(f"\n  Layer 1 mean α_crit: {l1_crit:.3f}")
        print(f"  Layer 2 mean α_crit: {l2_crit:.3f}")
        print(f"  Δα (L2 - L1): {delta_crit:+.3f}")
        if abs(delta_crit) > 0.05:
            print(f"  → SEPARATION: Layer 2 transitions {'later' if delta_crit > 0 else 'earlier'} than Layer 1")
        else:
            print(f"  → COINCIDENT: Layers transition together (|Δα| < 0.05)")

        output['sequences'][name] = seq_data

    # Summary
    prime_sep = output['sequences']['primes']['layer_separation']['delta']
    gue_sep = output['sequences']['gue']['layer_separation']['delta']
    pois_sep = output['sequences']['poisson']['layer_separation']['delta']

    print(f"\n=== SUMMARY ===")
    print(f"Layer separation Δα: Primes={prime_sep:+.3f}, GUE={gue_sep:+.3f}, Poisson={pois_sep:+.3f}")

    output['summary'] = {
        'prime_layer_separation': float(prime_sep),
        'gue_layer_separation': float(gue_sep),
        'poisson_layer_separation': float(pois_sep),
    }

    # Save
    out_path = Path(__file__).parent / 'data' / '3d_boundary_layers.json'
    with open(out_path, 'w') as f:
        json.dump(output, f, indent=2)
    print(f"\nSaved: {out_path}")


if __name__ == '__main__':
    main()

exec
/bin/bash -lc "sed -n '1,240p' tools/data/seme.json" in /opt/MM_D-ND
 succeeded in 0ms:
{
  "timestamp": "2026-04-30T03:36:51.567318",
  "piano": 60,
  "tensioni": [
    {
      "tipo": "confine_inesplorato",
      "id": "TRASCENDENZA_LIMITE",
      "claim": "La trascendenza e il limite attuale del modello. I punti fissi relazionali (non solo phi ma la rete di punti fissi tra osservabili) possono rivelare il vero grafo della realta e pattern nelle matrici. Il confine non e nella matematica - e nel passaggio tra piani.",
      "intensita": 0.9,
      "nota": "Input operatore 2026-04-10. Tocca: confine del modello, struttura relazionale dei punti fissi. Consecutio: quali punti fissi relazionali emergono dalle 21 tensioni attuali? Il grafo e gia nei dati?",
      "manuale": true,
      "porta": "sessione_interattiva",
      "condensato_ref": "A3,A10",
      "condensato_motivo": "Estende A3 (punto fisso singolo) a rete relazionale. Tocca A10 (dipolo) come caso speciale."
    },
    {
      "tipo": "scoperta",
      "id": "DUALITA_DIPOLARE_VS_ILLUSORIA",
      "claim": "Due tipi di dualita: (1) dipolare - generativa, il modello (det=-1), (2) illusoria - dispersiva, entropia (det=+1). Le regole incoerenti producono la seconda. La dualita illusoria e entropia come dispersione, non come informazione.",
      "intensita": 0.9,
      "nota": "Input operatore 2026-04-10. Tocca: entropia come dispersione illusoria vs generazione dipolare. Consecutio: nel Lab i domini Poisson (entropia massima) mostrano dualita illusoria? I domini GUE (strutturati) mostrano dualita dipolare? Il drift verso Poisson (POISSON_CONVERGENCE) e perdita di dualita dipolare?",
      "manuale": true,
      "porta": "sessione_interattiva",
      "condensato_ref": "A2,A10,F5",
      "condensato_motivo": "Discrimina due forme di det. A2 (confine) e la soglia. A10 (dipolo) e il tipo 1. F5 (frame) misura la struttura D-ND che e tipo 1."
    },
    {
      "tipo": "scoperta_numerica",
      "id": "METRIC_TENSOR",
      "claim": "Il tensore metrico dei primi è g=(p/2)². Nel tempo ln(p), è de Sitter 1+1D. z=-8.8 curvatura vs z=+22.5 rapporti ΔΓ.",
      "intensità": 0.9,
      "nota": "Sessione interattiva 4 aprile. Verificato su 78K primi.",
      "manuale": true,
      "porta": "sessione_interattiva",
      "condensato_ref": null,
      "condensato_motivo": "Risultato numerico verificato, non-tautologico"
    },
    {
      "tipo": "scoperta",
      "id": "TENSIONE_ENTITA",
      "claim": "La tensione non e un problema pratico - e un Entita. La tensione superflua crea latenza (tempo). Senza tensione superflua tutto e regolato da assiomi. Implicazione: le tensioni nel seme sono entita, non problemi da risolvere. Quelle superflue (det=+1) producono tempo/latenza.",
      "intensita": 0.85,
      "nota": "Input operatore 2026-04-10. Tocca: rapporto tensione/assioma. Operativamente: discriminare tensioni-entita (generative) da tensioni-superflue (dispersive) nel seme. Le 21 tensioni attuali - quante sono entita e quante latenza?",
      "manuale": true,
      "porta": "sessione_interattiva",
      "condensato_ref": "A5,A6",
      "condensato_motivo": "Il ciclo (A5) lavora con tensioni - ma se la tensione e entita, il ciclo non le risolve, le osserva. Lo zero mobile (A6) e la tensione senza latenza."
    },
    {
      "tipo": "confine_inesplorato",
      "id": "G_POTENZIALE_NULLA",
      "claim": "G e il potenziale di tutto come nulla - permette il prima e il dopo. Ci muoviamo come trascendenza dimensionale gravitazionale. G nel tetraedro non e una teoria tra le altre - e il potenziale che le rende possibili.",
      "intensita": 0.85,
      "nota": "Input operatore 2026-04-10. Tocca: ruolo di G nel tetraedro (T,Q,G,E). La fonte video_lp0RgZ6kQF8 dice: tensore metrico dentro la forma simplettica. G non e accanto a T,Q,E - e sotto. Consecutio: nei dati Lab, i ponti TxG e ExG hanno struttura diversa dai ponti TxQ?",
      "manuale": true,
      "porta": "sessione_interattiva",
      "condensato_ref": "A7,A10",
      "condensato_motivo": "A7 (singolarita come operatore) e G come potenziale. A10 (dipolo) opera sul piano che G rende possibile."
    },
    {
      "tipo": "confine_inesplorato",
      "id": "BOUNDARY",
      "claim": "8 domini GUE, 5 Poisson — il confine è il terzo incluso operativo",
      "intensità": 0.8,
      "nota": "Ciclo 60e (crossover): transizione di fase nel piano dipolare. Ciclo 60f (universalita'): la transizione di fase (lock+decay+flip) e' QUASI UNIVERSALE — 5/7 sequenze ordinate la mostrano. E' proprieta' del metodo (partial shuffle), non della sequenza. MA la DIREZIONE del lock e' diagnostica: GUE=-97.8, AR1=-94.7, Primi=-104.2, Periodic_24=-104.0, Logistic=+110.4. Due classi: 'repulsione' (~-97) e 'confinamento' (~-104). Primi coincidono con periodic 2,4,2,4 (Z/6Z) a 0.2 deg — il carattere dell'ordinamento primi E' il confinamento mod-6 (F2). Ciclo 60h (3D boundary): i due layer Markov (coppie→SR,L1; triple→SR2,triple_var) sono ACCOPPIATI al confine — alpha critico identico (0.334 per primi, delta=0.000; 0.311 vs 0.334 per GUE, delta=0.024). Il confine e' UNA transizione di fase, non due. I layer sono proiezioni dello stesso ordinamento. Poisson mostra separazione spuria (delta=-0.189, segnale=rumore). Consecutio: esiste una perturbazione NON-uniforme che disaccoppia i layer? (pair-preserving shuffle vs triple-preserving shuffle). Se si', i layer sono gradi di liberta' indipendenti. Se no, e' decomposizione spettrale di un solo asse.",
      "condensato_ref": "A9",
      "condensato_motivo": "A9 (terzo incluso) confermato: il confine e' un punto di transizione di fase, non interpolazione. Lo zero della magnitudine dipolare e' il terzo incluso operativo tra ordinato e disordinato.",
      "porta": "condensato"
    },
    {
      "tipo": "confine_inesplorato",
      "id": "PIANO_PRIMARIO_DUE_ASSIOMI",
      "claim": "I piani importanti sono il primario e i due assiomi che lo determinano nelle zone osservate. Non tutti gli assiomi operano ovunque - in ogni zona osservata, due assiomi determinano il piano primario.",
      "intensita": 0.8,
      "nota": "Input operatore 2026-04-10. Tocca: struttura locale degli assiomi. Consecutio: per ogni dominio Lab (primi, logistica, percolazione...) quali 2 assiomi del condensato sono operativi? Mappa assiomi x domini = grafo della realta locale.",
      "manuale": true,
      "porta": "sessione_interattiva",
      "condensato_ref": "A9,A14",
      "condensato_motivo": "A9 (terzo incluso) opera CON il piano. A14 (cascata) propaga - ma propaga cosa, se solo 2 assiomi sono attivi per zona?"
    },
    {
      "tipo": "simmetria_sospetta",
      "id": "META",
      "claim": "Tutti i 11 test passano — verifica che non stiamo testando solo tautologie",
      "intensità": 0.5,
      "nota": "Ciclo 60f: il test di universalita' CONFERMA parziale tautologia. La transizione di fase sotto partial shuffle e' proprietà del metodo (5/7 sequenze). Il contenuto NON-tautologico e': (1) la direzione locked (-104 per primi, -98 per GUE, +110 per logistica) — diagnostica; (2) la coincidenza Primi-Periodic a 0.2 deg — il carattere e' Z/6Z confinement. Il rapporto dL1/dSR=2.28 codifica il TIPO di confinamento. Consecutio: derivare 2.28 analiticamente dalla struttura mod-6 dei gap (F2 + Hardy-Littlewood pair correlations).",
      "condensato_ref": null,
      "condensato_motivo": "Ricorrente (3x in 2 giorni) e fuori dalla mappa",
      "porta": "novità"
    },
    {
      "tipo": "scoperta_numerica",
      "id": "DIPOLAR_ORDERING",
      "claim": "L'ordinamento dei gap primi non e un punto su un asse GUE-Poisson. E un dipolo: spacing_ratio spostato verso Poisson (Delta_tau=-0.12), lag1_acf spostato verso GUE (Delta_tau=+0.20). L'ordinamento crea coerenza 2x (std tau 0.09 vs shuffle 0.19). Il confine non interpola — ha struttura propria.",
      "intensita": 0.8,
      "nota": "Ciclo 60, 5 osservabili su 4 scale (1e4-1e7). Ciclo 60c: angolo -111 +/- 1 deg (GUE -97, separazione 14 deg z=170). Ciclo 60f (universalita'): Primi (-104.2) = Periodic 2,4,2,4 (-104.0) a 0.2 deg. L'ordinamento primi nel piano (SR,L1) E' la struttura Z/6Z (F2). Ciclo 60g (observable hunt): la memoria Markov si decompone in due layer visibili. Layer 1 (coppie, Mk1): SR e L1 — il piano dipolare. Layer 2 (triple, Mk2): SR2 (z=-9.4 sotto Mk1, z=-0.05 sotto Mk2), L2 (z=-5.3/Mk1, z=-0.2/Mk2), cond_entropy (z=-51/Mk1, z=2.4/Mk2). I layer sono ortogonali. Mk3 non aggiunge nulla sui 10 osservabili testati. Il terzo asse e' SR2 (next-nearest-neighbor spacing ratio). Consecutio: derivare SR2=0.4785 analiticamente dalla serie singolare Hardy-Littlewood per tripletti di primi.",
      "porta": "autoricerca",
      "condensato_ref": "A9,A10",
      "condensato_motivo": "A9 (terzo incluso): il confine non interpola tra GUE e Poisson, ha struttura propria. A10 (dipolo): i due osservabili order-sensitive sono il dipolo del confine."
    }
  ],
  "potenziale_bloccato": [],
  "varianza": [
    "Tensioni risolte: {'G_POTENZIALE_NULLA', 'TENSIONE_ENTITA', 'DUALITA_DIPOLARE_VS_ILLUSORIA', 'TRASCENDENZA_LIMITE', 'METRIC_TENSOR', 'PIANO_PRIMARIO_DUE_ASSIOMI'}"
  ],
  "filtro": {
    "promosse": 8,
    "filtrate": 0,
    "condensato_claims": 30
  },
  "direzione": "Esplorare il confine: 8 domini GUE, 5 Poisson — il confine è il terzo incluso operativo",
  "verifica": {
    "pass": 11,
    "fail": 0,
    "skip": 0,
    "total": 11
  },
  "fonti_consumate": 0,
  "fonti_esterne": [
    {
      "id": "video_lp0RgZ6kQF8",
      "title": "Equivalence between geometrical structures and entropy",
      "type": "video_digest",
      "keywords": [
        "geometry",
        "entropy",
        "symplectic form",
        "statistical mechanics",
        "quantum",
        "thermodynamics",
        "inner product",
        "Born rule",
        "metric tensor",
        "electromagnetic tensor"
      ],
      "content": "La geometria degli stati (classici e quantistici) e l'entropia sono la stessa struttura — invertibili. La forma simplettica conta le configurazioni. Il tensore metrico dello spaziotempo appare dentro la forma simplettica estesa. Il tensore elettromagnetico pure. Statistical mechanics non è costruita sopra alla meccanica — è la stessa cosa.",
      "teorie": [
        "T",
        "Q",
        "G",
        "E"
      ],
      "ponti_potenziali": [
        {
          "coppia": "TxQ",
          "ponte": "forma simplettica = entropia (invertibili)",
          "nota": "geometry is entropy and entropy is geometry"
        },
        {
          "coppia": "TxG",
          "ponte": "tensore metrico dentro la forma simplettica estesa",
          "nota": "geometria spaziotempo = geometria degli stati in posizione×velocità"
        },
        {
          "coppia": "ExT",
          "ponte": "tensore EM dentro la forma simplettica",
          "nota": "il campo EM conta stati in configurazione posizione×tempo"
        }
      ],
      "timestamp": "2026-04-02T08:23:13.991997"
    },
    {
      "id": "video_sDlZ-aY9GN4",
      "title": "Moving charges produce magnetic fields - Einstein relativity",
      "type": "video_digest",
      "keywords": [
        "magnetic field",
        "electric field",
        "length contraction",
        "time dilation",
        "Coulomb",
        "Lorentz",
        "reference frame",
        "electromagnetic"
      ],
      "content": "Il campo magnetico non esiste come entità separata — è il campo elettrico visto da un altro frame. La contrazione di Lorentz trasforma neutralità in carica. Due elettroni in movimento si separano più lentamente del previsto non per forza magnetica ma per dilatazione temporale. E e B sono manifestazioni dello stesso campo elettromagnetico. La relatività unifica.",
      "teorie": [
        "E",
        "R"
      ],
      "ponti_potenziali": [
        {
          "coppia": "ExR",
          "ponte": "cambio di frame — E e B sono lo stesso campo",
          "nota": "il 'ponte' non è l'onda EM, è il cambio di osservatore"
        }
      ],
      "timestamp": "2026-04-02T08:23:13.992016"
    },
    {
      "id": "video_OwDWOtFNsKQ",
      "title": "Thermodynamic Computing: Better than Quantum? | Guillaume Verdon (Extropic)",
      "type": "video_digest",
      "keywords": [
        "thermodynamic computing",
        "noise",
        "entropy",
        "extropic",
        "fluctuations",
        "information theory",
        "probability",
        "Boltzmann"
      ],
      "content": "Computing termodinamico: sfruttare le fluttuazioni termiche invece di combatterle. Il rumore non è nemico — è risorsa. Extropic costruisce hardware che usa l'entropia come motore computazionale. Connessione profonda tra termodinamica, teoria dell'informazione e probabilità.",
      "teorie": [
        "T",
        "Q"
      ],
      "ponti_potenziali": [
        {
          "coppia": "TxQ",
          "ponte": "noise come risorsa computazionale — fluttuazioni termiche = calcolo",
          "nota": "il vuoto quantistico (pieno di fluttuazioni) è il computer"
        }
      ],
      "timestamp": "2026-04-02T08:23:13.992019"
    },
    {
      "id": "video_j0wJBEZdwLs",
      "title": "What is a Laplace Transform - visual explanation",
      "type": "video_digest",
      "keywords": [
        "Laplace transform",
        "frequency",
        "damping",
        "s-plane",
        "complex",
        "exponential"
      ],
      "content": "La trasformata di Laplace come proiezione su esponenziali complesse. Il piano s = σ + iω combina smorzamento (reale) e oscillazione (immaginario). Connessione tra dominio temporale e dominio delle frequenze complesse.",
      "teorie": [
        "T",
        "Q",
        "R"
      ],
      "ponti_potenziali": [],
      "timestamp": "2026-04-02T08:23:13.992021"
    },
    {
      "id": "video_rZ2m1_q9lg0",

exec
/bin/bash -lc 'rg -n "partial|shuffle|pair|triple|surrogate|Markov" tools/exp_*.py | head -n 120' in /opt/MM_D-ND
 succeeded in 0ms:
tools/exp_acf_amplitude_scaling.py:50:    parser.add_argument('--n_surrogates', type=int, default=10)
tools/exp_acf_amplitude_scaling.py:82:        # Surrogate ACFs (shuffled)
tools/exp_acf_amplitude_scaling.py:83:        A_surrogates = []
tools/exp_acf_amplitude_scaling.py:84:        for _ in range(args.n_surrogates):
tools/exp_acf_amplitude_scaling.py:88:            A_surrogates.append(A_s)
tools/exp_acf_amplitude_scaling.py:90:        A_shuf_mean = np.mean(A_surrogates)
tools/exp_acf_amplitude_scaling.py:91:        A_shuf_std = np.std(A_surrogates) if np.std(A_surrogates) > 0 else 1e-10
tools/exp_alpha_stability.py:7:Reusable: accepts --n_primes, --n_windows, --n_surrogates, --max_lag.
tools/exp_alpha_stability.py:9:Null baseline: shuffled gaps within each window.
tools/exp_alpha_stability.py:81:    parser.add_argument("--n_surrogates", type=int, default=10, help="Shuffled surrogates per window")
tools/exp_alpha_stability.py:123:        # Shuffled surrogates
tools/exp_alpha_stability.py:126:        for _ in range(args.n_surrogates):
tools/exp_alpha_stability.py:128:            np.random.shuffle(shuf)
tools/exp_3d_boundary_layers.py:6:  - Layer 1 (pairs, Mk1) → SR, L1 (the dipolar plane)
tools/exp_3d_boundary_layers.py:7:  - Layer 2 (triples, Mk2) → SR2, triple_var (depth)
tools/exp_3d_boundary_layers.py:8:  - The crossover (partial shuffle) shows a phase transition in (SR, L1)
tools/exp_3d_boundary_layers.py:11:If same α_c → the boundary is 2D (Layer 2 follows Layer 1 = partial tautology)
tools/exp_3d_boundary_layers.py:49:def partial_shuffle(seq, alpha, rng):
tools/exp_3d_boundary_layers.py:57:    rng.shuffle(vals)
tools/exp_3d_boundary_layers.py:62:# --- Layer 1 observables (pair statistics) ---
tools/exp_3d_boundary_layers.py:73:# --- Layer 2 observables (triple statistics) ---
tools/exp_3d_boundary_layers.py:79:def obs_triple_var(gaps):
tools/exp_3d_boundary_layers.py:80:    """Variance of consecutive triple sums, normalized."""
tools/exp_3d_boundary_layers.py:81:    triples = gaps[:-2] + gaps[1:-1] + gaps[2:]
tools/exp_3d_boundary_layers.py:84:    return np.var(triples) / v
tools/exp_3d_boundary_layers.py:93:        'triple_var': obs_triple_var,
tools/exp_3d_boundary_layers.py:96:    # Full shuffle baseline (alpha=1.0)
tools/exp_3d_boundary_layers.py:99:        shuffled = partial_shuffle(gaps, 1.0, rng)
tools/exp_3d_boundary_layers.py:101:            baselines[name].append(fn(shuffled))
tools/exp_3d_boundary_layers.py:112:            s = partial_shuffle(gaps, alpha, rng)
tools/exp_3d_boundary_layers.py:198:    obs_names = ['SR', 'L1', 'SR2', 'triple_var']
tools/exp_3d_boundary_layers.py:199:    layer_map = {'SR': 'L1_pair', 'L1': 'L1_pair', 'SR2': 'L2_triple', 'triple_var': 'L2_triple'}
tools/exp_3d_boundary_layers.py:203:        'question': 'Do Layer 1 (pairs) and Layer 2 (triples) transition at the same critical alpha?',
tools/exp_3d_boundary_layers.py:232:        l2_crit = np.mean([seq_data['critical_alpha_50'][o] for o in ['SR2', 'triple_var']])
tools/exp_brody_flow.py:12:- Shuffle: same gaps per window, order destroyed → β_shuffle(N)
tools/exp_brody_flow.py:16:    python tools/exp_brody_flow.py [--n-max 2000000] [--window 5000] [--step 2000] [--n-shuffle 20]
tools/exp_brody_flow.py:92:def run(n_max=2_000_000, window=5000, step=2000, n_shuffle=20):
tools/exp_brody_flow.py:107:        "window": window, "step": step, "n_shuffle": n_shuffle,
tools/exp_brody_flow.py:129:        for _ in range(n_shuffle):
tools/exp_brody_flow.py:131:            np.random.shuffle(g_shuf)
tools/exp_brody_flow.py:139:        for _ in range(n_shuffle):
tools/exp_brody_flow.py:153:            "shuffle_beta_mean": float(np.mean(betas_shuf)),
tools/exp_brody_flow.py:154:            "shuffle_beta_std": float(np.std(betas_shuf)),
tools/exp_brody_flow.py:155:            "shuffle_r_mean": float(np.mean(rs_shuf)),
tools/exp_brody_flow.py:189:    # Is the slope significant vs shuffle?
tools/exp_brody_flow.py:191:    for _ in range(n_shuffle):
tools/exp_brody_flow.py:192:        shuf_betas = [w["shuffle_beta_mean"] + np.random.normal(0, w["shuffle_beta_std"])
tools/exp_brody_flow.py:196:    results["shuffle_slope_mean"] = float(np.mean(shuf_slopes))
tools/exp_brody_flow.py:197:    results["shuffle_slope_std"] = float(np.std(shuf_slopes))
tools/exp_brody_flow.py:199:        results["slope_z_vs_shuffle"] = float(
tools/exp_brody_flow.py:211:    parser.add_argument("--n-shuffle", type=int, default=20)
tools/exp_brody_flow.py:214:    results = run(args.n_max, args.window, args.step, args.n_shuffle)
tools/exp_brody_flow.py:229:    if "slope_z_vs_shuffle" in results:
tools/exp_brody_flow.py:230:        print(f"Slope z-score vs shuffle: {results['slope_z_vs_shuffle']:.2f}")
tools/exp_acf_stationarity.py:21:  5. Null baseline: shuffled (raw and normalized should both give ACF~0)
tools/exp_acf_stationarity.py:131:    parser.add_argument("--n_shuffles", type=int, default=15)
tools/exp_acf_stationarity.py:186:    print(f"PART 2: Shuffled baseline ({args.n_shuffles} surrogates)")
tools/exp_acf_stationarity.py:194:    for i in range(args.n_shuffles):
tools/exp_acf_stationarity.py:345:        'n_shuffles': args.n_shuffles,
tools/exp_acf_stationarity.py:350:        'shuffle_baseline': {
tools/exp_acf_z6z_mechanism.py:11:3. Factorial surrogates — Z/6Z with/without trend, trend without Z/6Z,
tools/exp_acf_z6z_mechanism.py:12:   residue-preserving shuffle.
tools/exp_acf_z6z_mechanism.py:15:    python tools/exp_acf_z6z_mechanism.py [--n_primes N] [--max_lag K] [--n_surrogates S]
tools/exp_acf_z6z_mechanism.py:77:    parser.add_argument('--n_surrogates', type=int, default=20)
tools/exp_acf_z6z_mechanism.py:147:    # TEST 3: Factorial surrogates
tools/exp_acf_z6z_mechanism.py:150:    print(f"TEST 3: Factorial surrogates ({args.n_surrogates} each)")
tools/exp_acf_z6z_mechanism.py:154:    NS = args.n_surrogates
tools/exp_acf_z6z_mechanism.py:156:    # 3a: Residue-preserving shuffle
tools/exp_acf_z6z_mechanism.py:157:    # For each residue class, shuffle gap values within the class
tools/exp_acf_z6z_mechanism.py:167:            np.random.shuffle(vals)
tools/exp_acf_z6z_mechanism.py:172:    print(f"  Residue-preserving shuffle:  xover={np.mean(xovers_rps):.1f}+/-{np.std(xovers_rps):.1f}, n_neg={np.mean(nnegs_rps):.1f}")
tools/exp_acf_z6z_mechanism.py:174:    # 3b: Full shuffle (baseline)
tools/exp_acf_z6z_mechanism.py:178:        np.random.shuffle(surr)
tools/exp_acf_z6z_mechanism.py:182:    print(f"  Full shuffle:                xover={np.mean(xovers_fs):.1f}+/-{np.std(xovers_fs):.1f}, n_neg={np.mean(nnegs_fs):.1f}")
tools/exp_acf_z6z_mechanism.py:206:    # 3e: Anti-correlated surrogates without trend (Gaussian AR(1) with acf1 matching primes)
tools/exp_acf_z6z_mechanism.py:243:        ("Residue-preserving shuffle", np.mean(xovers_rps), np.mean(nnegs_rps)),
tools/exp_acf_z6z_mechanism.py:244:        ("Full shuffle", np.mean(xovers_fs), np.mean(nnegs_fs)),
tools/exp_acf_z6z_mechanism.py:264:    #    if RPS crossover ~ shuffle → Z/6Z sequence doesn't matter (only ordering within class does)
tools/exp_acf_z6z_mechanism.py:266:    print(f"  Residue-preserving shuffle closer to raw ({xover_raw}) or shuffle ({np.mean(xovers_fs):.1f})?")
tools/exp_acf_z6z_mechanism.py:321:        'n_surrogates': args.n_surrogates,
tools/exp_acf_z6z_mechanism.py:339:        'test3_surrogates': {
tools/exp_acf_z6z_mechanism.py:340:            'residue_preserving_shuffle': {
tools/exp_acf_z6z_mechanism.py:346:            'full_shuffle': {
tools/exp_cross_observable_consistency.py:16:4. Null: shuffled primes (should show β_r ≈ β_Σ since no ordering channel)
tools/exp_cross_observable_consistency.py:144:# ── Measure shuffled primes (null) ──────────────────────────────────
tools/exp_cross_observable_consistency.py:145:print("\n=== PRIMES (shuffled, 20 trials) ===")
tools/exp_cross_observable_consistency.py:152:    np.random.shuffle(g_shuf)
tools/exp_cross_observable_consistency.py:155:    # Reconstruct positions from shuffled gaps
tools/exp_cross_observable_consistency.py:247:    "shuffle": {
tools/exp_boundary_growth.py:44:def shuffled_gaps_model(gaps, seed=None):
tools/exp_boundary_growth.py:47:    shuffled = gaps.copy()
tools/exp_boundary_growth.py:48:    rng.shuffle(shuffled)
tools/exp_boundary_growth.py:49:    return gap_ratio(shuffled)
tools/exp_boundary_growth.py:97:            r_shuf_list.append(shuffled_gaps_model(gaps, seed=42 + trial))
tools/exp_boundary_growth.py:114:            "r_shuffled_mean": float(r_shuf_mean),
tools/exp_boundary_growth.py:115:            "r_shuffled_std": float(r_shuf_std),
tools/exp_boundary_growth.py:117:            "excess_shuffled": float(excess_s),
tools/exp_boundary_growth.py:131:    excesses_s = [r["excess_shuffled"] for r in results]
tools/exp_boundary_growth.py:192:    # The shuffled-gaps test is the key anti-tautology check:
tools/exp_boundary_growth.py:193:    # if <r>_prime == <r>_shuffled, then the gap-ratio statistic
tools/exp_boundary_growth.py:197:    print(f"    <r>_prime - <r>_shuffled mean = {np.mean(excesses_s):+.6f}")
tools/exp_boundary_growth.py:218:            "excess_shuffled_mean": float(np.mean(excesses_s)),
tools/exp_boundary_growth.py:219:            "excess_shuffled_slope": float(coeffs_s[0]),
tools/exp_acf_range_universality.py:11:Domains: primes, GUE, GOE, Poisson, primes_shuffled.
tools/exp_acf_range_universality.py:12:Null: shuffled gaps (order destroyed, marginals preserved).
tools/exp_acf_range_universality.py:14:Reusable: --n_primes, --max_lag_prime, --max_lag_other, --n_surrogates.
tools/exp_acf_range_universality.py:95:    parser.add_argument("--n_surrogates", type=int, default=15)
tools/exp_acf_range_universality.py:148:    # Shuffled surrogate extended ACF
tools/exp_acf_range_universality.py:149:    print("\n  Shuffled surrogates (extended):")
tools/exp_acf_range_universality.py:151:    for i in range(args.n_surrogates):
tools/exp_acf_range_universality.py:153:        np.random.shuffle(sg)
tools/exp_acf_range_universality.py:176:    shuffled_gaps = gaps[:100_000].copy()
tools/exp_acf_range_universality.py:177:    np.random.shuffle(shuffled_gaps)
tools/exp_acf_range_universality.py:181:        'primes_shuffled': shuffled_gaps,
tools/exp_acf_range_universality.py:254:        for _ in range(args.n_surrogates):
tools/exp_acf_range_universality.py:256:            np.random.shuffle(sw)
tools/exp_brody_crossover.py:3:Experiment: Brody parameter crossover for primes vs Cramer surrogates.
tools/exp_brody_crossover.py:13:- Compare to Cramer surrogates (exponential gaps) at each scale
tools/exp_brody_crossover.py:59:def cramer_surrogate(densities, n_gaps):
tools/exp_brody_crossover.py:104:        # Cramer surrogates

codex
Uso surrogati Markov-k come perturbazione non uniforme: k=0 preserva solo la distribuzione, k=1 preserva transizioni a coppie binned, k=2 preserva transizioni a triple binned. Il criterio è retention rispetto a full shuffle: layer pair dovrebbe sopravvivere a Mk1, layer triple dovrebbe richiedere Mk2.
exec
/bin/bash -lc "python - <<'PY'
import numpy as np
from collections import Counter

rng=np.random.default_rng(1138)
N=30000
S=24
B=12

def primes_upto(n):
    sieve=np.ones(n+1,dtype=bool); sieve[:2]=False
    for i in range(2,int(n**0.5)+1):
        if sieve[i]: sieve[i*i::i]=False
    return np.where(sieve)[0]

def prime_gaps(N):
    p=primes_upto(N*20)[:N+1]
    return np.diff(p).astype(float)

def gue_gaps(N):
    out=[]
    dim=220
    while len(out)<N:
        H=rng.standard_normal((dim,dim)); H=(H+H.T)/2
        e=np.linalg.eigvalsh(H)
        g=np.diff(np.sort(e)); g=g[g>0]; g=g/np.mean(g)
        out.extend(g.tolist())
    return np.array(out[:N],float)

def poisson_gaps(N): return rng.exponential(1.0,N)

def SR(g):
    r=np.minimum(g[:-1],g[1:])/np.maximum(g[:-1],g[1:])
    return float(np.mean(r[np.isfinite(r)]))
def L1(g):
    x=g-np.mean(g); c=np.mean(x*x)
    return 0.0 if c==0 else float(np.mean(x[:-1]*x[1:])/c)
def SR2(g):
    r=np.minimum(g[:-2],g[2:])/np.maximum(g[:-2],g[2:])
    return float(np.mean(r[np.isfinite(r)]))
def TV(g):
    v=np.var(g)
    return 0.0 if v==0 else float(np.var(g[:-2]+g[1:-1]+g[2:])/v)
OBS={'SR':SR,'L1':L1,'SR2':SR2,'triple_var':TV}
LAY={'SR':'pair','L1':'pair','SR2':'triple','triple_var':'triple'}

def binned_model(g, order, n_bins=B):
    edges=np.percentile(g,np.linspace(0,100,n_bins+1)); edges[0]=g.min()-1e-9; edges[-1]=g.max()+1e-9
    b=np.digitize(g,edges)-1; b=np.clip(b,0,n_bins-1)
    pools={i:g[b==i] for i in range(n_bins)}
    if order==0: return b,pools,None
    trans={}
    for i in range(len(b)-order):
        st=tuple(b[i:i+order]); nx=int(b[i+order])
        trans.setdefault(st,Counter())[nx]+=1
    probs={st:(np.array(list(c.keys())),np.array(list(c.values()),float)/sum(c.values())) for st,c in trans.items()}
    return b,pools,probs

def markov_surrogate(g, order):
    b,pools,probs=binned_model(g,order)
    n=len(g); y=np.empty(n,float)
    if order==0:
        bins=rng.choice(B,size=n)
        for i,bb in enumerate(bins):
            pool=pools[int(bb)]
            if len(pool)==0: pool=g
            y[i]=pool[rng.integers(len(pool))]
        return y
    start=rng.integers(0,len(b)-order)
    state=tuple(int(x) for x in b[start:start+order])
    for j,bb in enumerate(state):
        pool=pools[bb] if len(pools[bb]) else g
        y[j]=pool[rng.integers(len(pool))]
    for i in range(order,n):
        if state in probs:
            keys,pp=probs[state]; bb=int(rng.choice(keys,p=pp))
        else:
            bb=int(rng.integers(B))
        pool=pools[bb] if len(pools[bb]) else g
        y[i]=pool[rng.integers(len(pool))]
        state=tuple((*state[1:],bb))
    return y

def full_shuffle(g):
    y=g.copy(); rng.shuffle(y); return y

def measure(label,g):
    real={k:f(g) for k,f in OBS.items()}
    sh={k:[] for k in OBS}
    sur={o:{k:[] for k in OBS} for o in [0,1,2]}
    for _ in range(S):
        y=full_shuffle(g)
        for k,f in OBS.items(): sh[k].append(f(y))
        for o in [0,1,2]:
            y=markov_surrogate(g,o)
            for k,f in OBS.items(): sur[o][k].append(f(y))
    base={k:np.mean(v) for k,v in sh.items()}
    print(f\"\\n=== {label} ===\")
    print(\"obs layer real shuffle Mk0_ret Mk1_ret Mk2_ret Mk1_z_vs_shuf Mk2_z_vs_Mk1\")
    layer_ret={o:{'pair':[],'triple':[]} for o in [0,1,2]}
    for k in OBS:
        denom=real[k]-base[k]
        vals=[]
        for o in [0,1,2]:
            m=np.mean(sur[o][k]); ret=(m-base[k])/denom if abs(denom)>1e-12 else np.nan
            vals.append(ret); layer_ret[o][LAY[k]].append(ret)
        z1=(np.mean(sur[1][k])-base[k])/(np.std(sh[k])+1e-12)
        z21=(np.mean(sur[2][k])-np.mean(sur[1][k]))/(np.std(sur[1][k])+1e-12)
        print(f\"{k:10s} {LAY[k]:6s} {real[k]:+.6f} {base[k]:+.6f} {vals[0]:+7.3f} {vals[1]:+7.3f} {vals[2]:+7.3f} {z1:+10.2f} {z21:+10.2f}\")
    for o in [0,1,2]:
        pm=np.nanmean(layer_ret[o]['pair']); tm=np.nanmean(layer_ret[o]['triple'])
        print(f\"Mk{o} layer mean retention: pair={pm:+.3f} triple={tm:+.3f} delta(triple-pair)={tm-pm:+.3f}\")

measure('PRIMES',prime_gaps(N))
measure('GUE',gue_gaps(N))
measure('POISSON',poisson_gaps(N))
PY" in /opt/MM_D-ND
codex
L’esperimento sta ancora calcolando GUE e surrogati Markov. Separa bene il rischio di tautologia: i numeri decisivi saranno le retention medie pair/triple per Mk1 e Mk2, non il valore assoluto degli osservabili.
 succeeded in 60467ms:

=== PRIMES ===
obs layer real shuffle Mk0_ret Mk1_ret Mk2_ret Mk1_z_vs_shuf Mk2_z_vs_Mk1
SR         pair   +0.468411 +0.490797  +2.520  +0.947  +0.957     -14.62      -0.16
L1         pair   -0.045317 +0.000351  -0.009  +0.995  +1.010      -7.76      -0.12
SR2        triple +0.485363 +0.490154 +11.634  -1.444  +0.825      +6.56      -6.36
triple_var triple +2.782778 +2.996461  -0.052  +0.794  +0.956      -6.85      -1.49
Mk0 layer mean retention: pair=+1.256 triple=+5.791 delta(triple-pair)=+4.535
Mk1 layer mean retention: pair=+0.971 triple=-0.325 delta(triple-pair)=-1.296
Mk2 layer mean retention: pair=+0.984 triple=+0.890 delta(triple-pair)=-0.093

=== GUE ===
obs layer real shuffle Mk0_ret Mk1_ret Mk2_ret Mk1_z_vs_shuf Mk2_z_vs_Mk1
SR         pair   +0.531011 +0.544039  +0.013  +1.031  +1.014     -11.57      +0.12
L1         pair   +0.072102 +0.001054  -0.026  -0.091  -0.085      -1.65      +0.07
SR2        triple +0.554146 +0.543792  +0.039  +0.512  +0.956      +4.53      +2.85
triple_var triple +3.659132 +3.002087  -0.005  +0.025  +0.307      +0.90      +6.83
Mk0 layer mean retention: pair=-0.006 triple=+0.017 delta(triple-pair)=+0.023
Mk1 layer mean retention: pair=+0.470 triple=+0.269 delta(triple-pair)=-0.201
Mk2 layer mean retention: pair=+0.464 triple=+0.632 delta(triple-pair)=+0.167

=== POISSON ===
obs layer real shuffle Mk0_ret Mk1_ret Mk2_ret Mk1_z_vs_shuf Mk2_z_vs_Mk1
SR         pair   +0.387853 +0.386819  -0.407  +0.538  +0.806      +0.47      +0.14
L1         pair   +0.006511 +0.000361  -0.181  +0.575  +0.749      +0.85      +0.17
SR2        triple +0.386191 +0.386427  -2.698  -1.441  +3.535      +0.28      -0.51
triple_var triple +3.035675 +2.997251  -0.025  +0.396  +0.772      +0.66      +0.53
Mk0 layer mean retention: pair=-0.294 triple=-1.361 delta(triple-pair)=-1.067
Mk1 layer mean retention: pair=+0.557 triple=-0.523 delta(triple-pair)=-1.079
Mk2 layer mean retention: pair=+0.777 triple=+2.153 delta(triple-pair)=+1.376

codex
Risultato grezzo: sui primi Mk1 conserva quasi interamente il layer pair (retention media 0.971), mentre SR2 non sopravvive a Mk1 e riappare con Mk2. triple_var è meno selettivo: conserva già 0.794 sotto Mk1, quindi non va usato come claim assoluto di profondità. Provo ora a scrivere il report richiesto; il filesystem risulta configurato read-only, quindi se la scrittura viene bloccata lascio il contenuto completo qui.
2026-05-04T11:41:46.089926Z ERROR codex_core::tools::router: error=patch rejected: writing is blocked by read-only sandbox; rejected by user approval settings
codex
Ho eseguito l’esperimento, ma la scrittura del report è stata bloccata dal sandbox read-only: `tools/data/reports/agent_20260504_1138.md` non è stato creato. Ho tentato `apply_patch`; errore: scrittura bloccata dai permessi.

Contenuto del report da salvare:

```markdown
# Agent Report — Markov-2 Decouples SR2, But Triple Variance Is Not a Clean Depth Observable
**Date**: 2026-05-04 11:38
**Piano**: 60
**Tension explored**: META (0.5) + BOUNDARY (0.8) + DIPOLAR_ORDERING (0.8)

## Claim Under Test
> Dopo il report 2026-05-04 09:01: se lo shuffle uniforme accoppia layer pair e layer triple, esiste una perturbazione non-uniforme che li disaccoppia?

## Question
I surrogati Markov-k separano il layer pair (SR, L1) dal layer triple (SR2, triple_var)?

## Experiment Design
- Metodo: surrogati Markov-k binned.
- k=0 preserva distribuzione marginale.
- k=1 preserva transizioni tra coppie di bin.
- k=2 preserva transizioni tra triple di bin.
- Retention: `(surrogate_mean - full_shuffle_mean) / (real - full_shuffle_mean)`.
- Scope: N=30000, 24 surrogati per k, 12 bin.
- Domini: prime gaps, GUE spacings, Poisson iid.

## Results

| Domain | Mk0 pair | Mk0 triple | Mk1 pair | Mk1 triple | Mk2 pair | Mk2 triple |
|--------|----------|------------|----------|------------|----------|------------|
| Primes | +1.256 | +5.791 | +0.971 | -0.325 | +0.984 | +0.890 |
| GUE | -0.006 | +0.017 | +0.470 | +0.269 | +0.464 | +0.632 |
| Poisson | -0.294 | -1.361 | +0.557 | -0.523 | +0.777 | +2.153 |

Prime observables:

| Observable | Layer | Real | Shuffle | Mk1 ret | Mk2 ret |
|------------|-------|------|---------|---------|---------|
| SR | pair | +0.468411 | +0.490797 | +0.947 | +0.957 |
| L1 | pair | -0.045317 | +0.000351 | +0.995 | +1.010 |
| SR2 | triple | +0.485363 | +0.490154 | -1.444 | +0.825 |
| triple_var | triple | +2.782778 | +2.996461 | +0.794 | +0.956 |

## Key Findings
1. SR2 è selettivamente Markov-2 per i gap primi: Mk1 preserva SR/L1 quasi interamente, ma non SR2; Mk2 lo ripristina.
2. triple_var non è un osservabile di profondità pulito: conserva già 0.794 sotto Mk1.
3. Il risultato 09:01 non è falsificato: lo shuffle uniforme accoppia i layer; i surrogati Markov-k mostrano che SR2 può essere disaccoppiato.
4. Poisson è rumore da denominatore piccolo, non controesempio strutturale.

## Verdict
**CONSTRAINT on BOUNDARY + DIPOLAR_ORDERING**: il confine è una transizione unica sotto partial shuffle uniforme, ma non è indivisibile sotto perturbazioni selettive. Per i gap primi, SR2 è il coordinato operativo di profondità; triple_var va demossa ad ausiliario.

## Bicono della scoperta
- **Due radici**: Markov-1 plane (SR, L1) e Markov-2 depth (SR2).
- **Singolare**: la sequenza ordinata dei gap primi prima della scelta della finestra osservativa.
- **Invariante di passaggio**: SR2 cade sotto Mk1 e riemerge sotto Mk2.
- **Campo di possibilità**: possibile usare SR2 come coordinata pulita del depth; non-possibile usare triple_var da sola come evidenza di profondità indipendente.

## Files
- Esperimento: inline Python.
- Report previsto: `tools/data/reports/agent_20260504_1138.md`
```

Non ho aggiornato `tools/data/seme.json` per lo stesso vincolo read-only.
tokens used
109,127
Ho eseguito l’esperimento, ma la scrittura del report è stata bloccata dal sandbox read-only: `tools/data/reports/agent_20260504_1138.md` non è stato creato. Ho tentato `apply_patch`; errore: scrittura bloccata dai permessi.

Contenuto del report da salvare:

```markdown
# Agent Report — Markov-2 Decouples SR2, But Triple Variance Is Not a Clean Depth Observable
**Date**: 2026-05-04 11:38
**Piano**: 60
**Tension explored**: META (0.5) + BOUNDARY (0.8) + DIPOLAR_ORDERING (0.8)

## Claim Under Test
> Dopo il report 2026-05-04 09:01: se lo shuffle uniforme accoppia layer pair e layer triple, esiste una perturbazione non-uniforme che li disaccoppia?

## Question
I surrogati Markov-k separano il layer pair (SR, L1) dal layer triple (SR2, triple_var)?

## Experiment Design
- Metodo: surrogati Markov-k binned.
- k=0 preserva distribuzione marginale.
- k=1 preserva transizioni tra coppie di bin.
- k=2 preserva transizioni tra triple di bin.
- Retention: `(surrogate_mean - full_shuffle_mean) / (real - full_shuffle_mean)`.
- Scope: N=30000, 24 surrogati per k, 12 bin.
- Domini: prime gaps, GUE spacings, Poisson iid.

## Results

| Domain | Mk0 pair | Mk0 triple | Mk1 pair | Mk1 triple | Mk2 pair | Mk2 triple |
|--------|----------|------------|----------|------------|----------|------------|
| Primes | +1.256 | +5.791 | +0.971 | -0.325 | +0.984 | +0.890 |
| GUE | -0.006 | +0.017 | +0.470 | +0.269 | +0.464 | +0.632 |
| Poisson | -0.294 | -1.361 | +0.557 | -0.523 | +0.777 | +2.153 |

Prime observables:

| Observable | Layer | Real | Shuffle | Mk1 ret | Mk2 ret |
|------------|-------|------|---------|---------|---------|
| SR | pair | +0.468411 | +0.490797 | +0.947 | +0.957 |
| L1 | pair | -0.045317 | +0.000351 | +0.995 | +1.010 |
| SR2 | triple | +0.485363 | +0.490154 | -1.444 | +0.825 |
| triple_var | triple | +2.782778 | +2.996461 | +0.794 | +0.956 |

## Key Findings
1. SR2 è selettivamente Markov-2 per i gap primi: Mk1 preserva SR/L1 quasi interamente, ma non SR2; Mk2 lo ripristina.
2. triple_var non è un osservabile di profondità pulito: conserva già 0.794 sotto Mk1.
3. Il risultato 09:01 non è falsificato: lo shuffle uniforme accoppia i layer; i surrogati Markov-k mostrano che SR2 può essere disaccoppiato.
4. Poisson è rumore da denominatore piccolo, non controesempio strutturale.

## Verdict
**CONSTRAINT on BOUNDARY + DIPOLAR_ORDERING**: il confine è una transizione unica sotto partial shuffle uniforme, ma non è indivisibile sotto perturbazioni selettive. Per i gap primi, SR2 è il coordinato operativo di profondità; triple_var va demossa ad ausiliario.

## Bicono della scoperta
- **Due radici**: Markov-1 plane (SR, L1) e Markov-2 depth (SR2).
- **Singolare**: la sequenza ordinata dei gap primi prima della scelta della finestra osservativa.
- **Invariante di passaggio**: SR2 cade sotto Mk1 e riemerge sotto Mk2.
- **Campo di possibilità**: possibile usare SR2 come coordinata pulita del depth; non-possibile usare triple_var da sola come evidenza di profondità indipendente.

## Files
- Esperimento: inline Python.
- Report previsto: `tools/data/reports/agent_20260504_1138.md`
```

Non ho aggiornato `tools/data/seme.json` per lo stesso vincolo read-only.
