Incredibile: i numeri hanno segreti nascosti nelle loro somme!

Come un chimico e un’intelligenza artificiale hanno scoperto che i numeri interi si comportano come particelle quantistiche

---

Tutto è partito da una domanda semplice: e se guardassimo i numeri interi non come prodotti di primi, ma come somme?

Il Teorema Fondamentale dell’Aritmetica ci dice che ogni numero si scrive in un solo modo come prodotto di primi. Il 12 è 2x2x3, punto.
Ma come somma di primi, il 12 si può scrivere in tanti modi diversi: 2+2+2+2+2+2, oppure 5+7, oppure 2+3+7, oppure 2+5+5.
È questa ricchezza che abbiamo deciso di esplorare.

---

Bosoni e fermioni — no, non è roba da università

In fisica quantistica le particelle si dividono in due famiglie.

I bosoni sono le particelle socievoli: possono stare nello stesso stato in qualsiasi numero. I fotoni della luce sono bosoni — per questo i laser funzionano: miliardi di fotoni tutti nello stesso stato, tutti insieme.

I fermioni sono le particelle solitarie: nessuna coppia può trovarsi nello stesso stato. Gli elettroni sono fermioni. Questo è il motivo per cui la materia è solida — gli elettroni non si accatastano tutti sul livello più basso, ma occupano livelli diversi. Se non fosse così, ogni atomo collasserebbe su se stesso.

Bene. Adesso porta questa idea sui numeri.

Prendete il numero 7. Può essere scritto come somma di primi in tre modi:

7 = 2+2+3 — il 2 si ripete. Stato bosonico.
7 = 2+5 — tutti diversi. Stato fermionico.
7 = 7 — un primo solo. Stato fermionico.

Chiamiamo fermionico uno stato in cui tutti i primi sono distinti.
Chiamiamo bosonico un sistema che ammette qualsiasi stato, con o senza ripetizioni.
Il sistema fermionico è quindi un sottoinsieme di quello bosonico.

---

I mattoncini di tutto: 2 e 3

Prova a costruire qualsiasi numero usando solo 2 e 3:

5 = 2+3 8 = 2+2+2+2 oppure 2+3+3 11 = 2+2+2+2+3 17 = 2+2+2+2+3+3+3

Funziona sempre, per qualsiasi numero maggiore o uguale a 2. Nessun’altra coppia di primi ha questa proprietà — prova con 2 e 5: non riesci a fare il 3. Prova con 3 e 5: non riesci a fare il 4.

Il 2 e il 3 sono anche gli unici numeri con un solo stato: il 2 si scrive solo come {2}, il 3 solo come {3}. Non si scompongono in nulla di più piccolo — sono i mattoncini irriducibili del sistema.

Il 5 invece ha già due stati: {5} e {2+3}. Però nessuno dei due ha ripetizioni, quindi entrambi stati fermionici.

---

Il 4 e il 6 sono strani

Quanti stati fermionici ha il 4?

4 = 2+2 — il 2 si ripete. Bosonico, non fermionico. Non esistono altri modi.

Il 4 non ha nessuno stato fermionico. Lo stesso vale per il 6:

6 = 3+3 — il 3 si ripete. Bosonico. 6 = 2+2+2 — il 2 si ripete. Bosonico. Nient’altro.

Il 4 e il 6 sono i soli due numeri composti senza alcuno stato fermionico. Per tutti gli altri esiste sempre almeno un modo di scriverli come somma di primi tutti diversi. Questo non è ovvio — e per i numeri grandi segue da un teorema dimostrato nel 2013 da un matematico di nome Harald Helfgott, che ha risolto un problema aperto da secoli.

---

L’entropia di un numero

Quanti stati ha un numero? Il 7 ne ha 3, il 10 ne ha 5. Più stati ha un numero, più è ricco — più disordinato, direbbero i fisici. Questa ricchezza si misura con l’entropia:

S = ln(numero di stati)

Il logaritmo serve per rendere l’entropia una quantità che cresce in modo ragionevole — se raddoppio gli stati, l’entropia non raddoppia, cresce di ln(2).

È la stessa formula che Boltzmann scrisse sulla sua lapide: S = k · ln(Ω), dove Ω è il numero di microstati e k è una costante. Noi poniamo k=1 per semplicità.

Per il 10:

• S_B = ln(5) ≈ 1.61 perché ha 5 stati bosonici.
• S_F = ln(2) ≈ 0.69 perché ha 2 stati fermionici.

---

Quanto è bosonico un numero? Il parametro lambda

Avendo due entropie — una bosonica S_B e una fermionica S_F — possiamo costruire una misura di bosonicità.
La scelta naturale è un rapporto tra la loro differenza e la loro somma:

lambda = (S_B – S_F) / (S_B + S_F)

Vale 0 quando le due entropie sono uguali — cioè quando bosoni e fermioni hanno lo stesso numero di stati.
Vale 1 nel caso estremo in cui non esistono stati fermionici.

Per il 10:
lambda = (1.61 – 0.69) / (1.61 + 0.69) = 0.92 / 2.30 ≈ 0.40

Tre numeri hanno lambda = 0: il 2, il 3 e il 5. Ma per ragioni diverse.

Il 2 e il 3 hanno un solo stato ciascuno — {2} e {3} — identico in entrambi i sistemi. Non c’è nulla da distinguere.

Il 5 ha due stati — {5} e {2,3} — ma entrambi sono fermionici per natura: nessuno ha ripetizioni.
Quindi il sistema bosonico e quello fermionico del 5 vedono esattamente le stesse cose. S_B = S_F = ln(2), lambda = 0.

Per n grande converge lentamente verso un valore che sembra avvicinarsi a 1/4=0.25, anche se la convergenza è così lenta che non possiamo escludere un limite diverso. Quello che è certo è la direzione: tutti i numeri, primi e composti, tendono allo stesso valore. Mediamente, per n tra 7 e 4000, i primi hanno lambda pari a 0.279 e i composti 0.274;
i primi leggermente più alti, con una differenza che si assottiglia rapidamente e diventa quasi impercettibile per n grande.

---

I bosoni rivelano i primi — ma i fermioni li affinano

Per ogni numero n calcoliamo una quantità semplice: il numero medio di addendi dei suoi stati bosonici.
Per il 10 avevamo cinque stati — con 5, 4, 3, 2, 2 addendi rispettivamente — quindi la media è (5+4+3+2+2)/5 = 3.2.

Facendo lo stesso per tutti i numeri da 2 in poi otteniamo una sequenza di punti che, messi su un grafico, formano una curva che sale.
Non sale in modo uniforme — ha irregolarità locali che rispecchiano la struttura aritmetica dei numeri.

Perché la media degli addendi e non, per esempio, l’entropia diretta? Perché la media degli addendi è una quantità liscia e regolare — cresce in modo graduale senza i salti bruschi del logaritmo.

L’entropia bosonica S_B= ln(Ω_B) fa un salto ogni volta che si aggiunge un nuovo stato, creando discontinuità che rendono la curvatura irregolare e difficile da interpretare.
La media degli addendi, invece, cambia in modo continuo.

La curvatura è un parametro che ci fa capire se la nostra curva sta accelerando o decelerando.

curvatura(n) = valore(n+1) – 2·valore(n) + valore(n-1)

È la derivata seconda discreta — la stessa formula dell’accelerazione in fisica.

Facciamo un esempio concreto con n=7. La media degli addendi bosonici vale: media(6)=2.5, media(7)=2.0, media(8)=3.0. Quindi:

curvatura bosonica(7) = 3.0 – 2·2.0 + 2.5 = +1.5 → positiva. Il 7 è primo.

L’entropia fermionica vale: S_F(6) non esiste (il 6 ha zero stati fermionici), S_F(8)=0, S_F(9)=0. La curvatura fermionica non è calcolabile per n=7 perché il 6 è fermionicamente vuoto — un altro effetto della stranezza del 4 e del 6.

Prendiamo n=12. Media(11)=3.33, media(12)=3.86, media(13)=3.67:

curvatura bosonica(12) = 3.67 – 2·3.86 + 3.33 = -0.71 → negativa. Il 12 è composto.

Il pattern è sistematico: i numeri primi tendono ad avere curvatura bosonica positiva (+0.28 in media), i composti negativa (-0.10 in media).
Il 65% dei primi ha curvatura positiva, il 66% dei composti ha curvatura negativa.

Quando invece calcoliamo la curvatura dell’entropia fermionica S_F, il pattern si sgonfia: i primi hanno curvatura positiva solo nel 29% dei casi, i composti curvatura negativa solo nel 52% — quasi il caso puro. La curvatura fermionica è quasi cieca alla primalità.

Ora costruiamo una sequenza partendo da n=2. La regola: ad ogni passo scegli il salto che ti porta al numero con la curvatura più alta.
I salti disponibili sono solo numeri primi — 2, 3, 5, 7, 11, 13… — perché ogni passo rappresenta l’aggiunta di un mattoncino primo alla somma.
Non puoi saltare di 1 perché 1 non è primo. Non puoi saltare di 4 o 6 perché non sono primi.

Abbiamo quattro combinazioni: la curvatura può essere bosonica o fermionica, e i passi possono essere bosonici (ripetibili) o fermionici (ognuno usabile al massimo una volta — principio di Pauli applicato al percorso).

Combinazione 1 — curvatura bosonica, passi bosonici:

Da n=2:
salto di 3 → n=5, curvatura = +1.50
salto di 5 → n=7, curvatura = +1.50
salto di 11 → n=13, curvatura = +0.72
Il 3 e il 5 pareggiano — scegliamo il salto più piccolo, arriviamo a n=5.

Da n=5:
salto di 2 → n=7, curvatura = +1.50
salto di 17 → n=22, curvatura = +0.12
Il 2 vince — andiamo a n=7.

Da n=7:
salto di 19 → n=26, curvatura = +0.21
salto di 2 → n=9, curvatura = +0.20
Il 19 vince di misura — andiamo a n=26.

Sequenza: 2 → 5 → 7 → 26 → 31 → 36 → 43 → 48 → 50 → 53 → 60 → 67 → 69
Primi: 7 su 13 = 53.8%

Combinazione 2 — curvatura bosonica, passi fermionici:

Stessa curvatura bosonica, ma ogni passo è usabile al massimo una volta.

Da n=2: scelgo salto 3 → n=5. Il 3 è consumato.
Da n=5: scelgo salto 2 → n=7. Il 2 è consumato.

Da n=7: il 2 è consumato.
Scelgo salto 19 → n=26. Il 19 è consumato.

Da n=26:
scelgo salto 5 → n=31. Il 5 è consumato.
Da n=31: scelgo salto 17 → n=48. Il 17 è consumato.

Sequenza: 2 → 5 → 7 → 26 → 31 → 48 → 71 → 84 → 95
Primi: 5 su 9 = 55.6%
Passi usati: {2, 3, 5, 11, 13, 17, 19, 23} — tutti primi.

Combinazione 3 — curvatura fermionica, passi bosonici:

Da n=2:
salto di 7 → n=9, curvatura fermionica = +0.693
salto di 17 → n=19, curvatura fermionica = +0.575
salto di 13 → n=15, curvatura fermionica = +0.405
Il 9 vince — è un composto.

Da n=9:
salto di 2 → n=11, curvatura fermionica = +1.386
L’ 11 vince — andiamo a n=11.

Da n=11:
salto di 19 → n=30, curvatura fermionica = +0.560
salto di 2 → n=13, curvatura fermionica = 0.000
Il 30 vince — è un composto.

La traiettoria abbandona subito il primo.
Sequenza: 2 → 9 → 11 → 30 → 32 → 35 → 40 → 45 → 56 → 58 → 65 → 67 → 70
Primi: 3 su 13 = 23.1%

Combinazione 4 — curvatura fermionica, passi fermionici:

Stessa curvatura fermionica, passi non ripetibili.
Sequenza: 2 → 9 → 11 → 30 → 35 → 38 → 51 → 62 → 91
Primi: 2 su 9 = 22.2%

---

Quanti numeri primi ci aspettiamo di trovare per caso in un intervallo di numeri? Non tutti i numeri sono primi — man mano che i numeri crescono, i primi diventano sempre più rari.
Ma quanto rari?

Gauss fu il primo a rispondere, a 15 anni, contando i primi nelle tavole logaritmiche. La sua intuizione: intorno a un numero n, la densità dei primi è approssimativamente 1/ln(n) — uno su logaritmo naturale di n. Nell’intervallo che esploriamo, questo dà circa il 27%. È la stima più semplice e ancora oggi la più usata.

Legendre affinò la formula aggiungendo una costante: la densità è circa 1/(ln(n) − 1.08). Migliora la stima per n piccoli, ma non cambia la struttura.

Riemann nel 1859 fece il salto definitivo. Invece di una formula approssimata, scrisse una formula esatta — ma con un prezzo: dipende dalla posizione degli zeri della funzione zeta. La sua formula dice che la distribuzione dei primi oscilla intorno alla stima di Gauss con correzioni precise, una per ogni zero. Conoscere gli zeri significa conoscere la distribuzione esatta dei primi. Ecco perché l’Ipotesi di Riemann — che tutti gli zeri stiano sulla linea Re(s)=1/2 — vale un milione di dollari.

Nel nostro caso, il riferimento naturale è la stima di Gauss: ~27% nell’intervallo esplorato. Tutto quello che supera questo valore è una firma dei primi che emerge dalla struttura del sistema.

Il riepilogo è immediato:

Curv. bosonica + passi bosonici: 53.8% — quasi doppio di Gauss
Curv. bosonica + passi fermionici: 55.6% — quasi doppio di Gauss
Curv. fermionica + passi bosonici: 23.1% — sotto Gauss
Curv. fermionica + passi fermionici: 22.2% — sotto Gauss

---

La conclusione è netta: quello che conta è la curvatura, non i passi.

Con curvatura bosonica si trovano i primi — sempre, sia con passi bosonici che fermionici.
Con curvatura fermionica si perdono — sempre, in entrambi i casi sotto Gauss.

C’è però una simmetria sottile: i passi fermionici migliorano leggermente con la curvatura bosonica (53.8% → 55.6%), e peggiorano leggermente con quella fermionica (23.1% → 22.2%).

La firma dei primi è nella curvatura bosonica — e solo lì. I bosoni la trovano in entrambi i casi.
I fermioni la perdono in entrambi i casi. Ma il vincolo fermionico — il principio di Pauli applicato ai passi — la rende più nitida quando la firma c’è, e più debole quando non c’è.

I numeri primi stanno in un punto di equilibrio tra i due sistemi: la curvatura bosonica li rivela, il vincolo fermionico li affina. Questa asimmetria è il filo rosso di tutto il framework — e tornerà, da direzioni inaspettate, fino all’ultima sezione del post.

---

Le tre funzioni di partizione

Abbiamo visto che i bosoni riveano i primi attraverso la curvatura. Ma da dove viene questa struttura? Per capirlo dobbiamo guardare dentro le funzioni di partizione.
Ad ogni stato assegniamo un’energia uguale al rango del primo nella sequenza: il 2 è il primo primo — rango 1. Il 3 è il secondo — rango 2. Il 5 è il terzo — rango 3. Il 7 è il quarto — rango 4.
E così via.

L’energia di uno stato è la somma dei ranghi dei suoi pezzi.
Torniamo al numero 10:
{2,2,2,2,2} — ranghi 1+1+1+1+1 — bosonico
{2,2,3,3} — ranghi 1+1+2+2 — bosonico
{2,3,5} — ranghi 1+2+3 — fermionico
{3,7} — ranghi 2+4 — fermionico
{5,5} — ranghi 3+3 — bosonico

Lo stato {2,2,2,2,2} ha l’energia più bassa — è il più pigro. Tutti gli altri costano di più.

La funzione di partizione Z è la somma di tutti gli stati, dove ogni stato contribuisce con un peso che dipende dalla sua energia.
Gli stati a bassa energia pesano di più, quelli ad alta energia pesano di meno.
Quanto di più? Dipende da un parametro che chiamiamo beta — una specie di manopola.

Con la manopola beta su valori grandi — temperatura bassa — il peso dello stato {2,2,2,2,2} domina su tutti gli altri: il sistema preferisce lo stato più pigro, come l’acqua che a bassa temperatura diventa ghiaccio. Con la manopola beta su valori piccoli — temperatura alta — tutti i pesi si equivalgono, tutti gli stati diventano ugualmente probabili, come il vapore in cui le molecole volano libere in tutte le direzioni.

Ma nel nostro framework beta non è fissato da nessuna fisica reale. Non c’è un termometro. Beta è una manopola libera — è l’impalcatura che usiamo per costruire l’edificio.
Quello che ci interessa sono le proprietà che valgono per qualsiasi beta.

Costruiamo ora Z pezzo per pezzo, guardando un primo per volta.

Prendiamo lo stato {2,3,5} del numero 10. Il suo peso è il prodotto dei contributi di ciascun primo, ognuno pesato dal suo rango:

• il 2 compare 1 volta, rango 1 → e^(-1·beta)
• il 3 compare 1 volta, rango 2 → e^(-2·beta)
• il 5 compare 1 volta, rango 3 → e^(-3·beta)
• il 7 compare 0 volte → e^(0) = 1
• l’11 compare 0 volte → e^(0) = 1
• tutti gli altri → 1

Il peso dello stato è e^(-1·beta) · e^(-2·beta) · e^(-3·beta) · 1 · 1 · 1 · …

Prendiamo ora lo stato {7,11}. I contributi sono:

• il 2 compare 0 volte → 1
• il 3 compare 0 volte → 1
• il 5 compare 0 volte → 1
• il 7 compare 1 volta, rango 4 → e^(-4·beta)
• l’11 compare 1 volta, rango 5 → e^(-5·beta)
• dall 13 in poi → tutti 1

Il peso dello stato è 1 · 1 · 1 · e^(-4·beta) · e^(-5·beta) · 1 · 1 · …

I contributi dei singoli primi vanno moltiplicati tra loro perché i primi sono indipendenti — la scelta di quante volte usare il 7 non dipende dalla scelta di quante volte usare l’11, esattamente come due dadi indipendenti. In matematica, quando due sistemi non si influenzano il loro contributo totale è il prodotto dei contributi individuali, non la somma.

Ora, per ogni primo p, quante volte può comparire in un sistema bosonico? Zero, una, due, tre volte… all’infinito. La somma di tutti questi contributi per il primo p è:
S = 1 + e^(-rango(p)·beta) + e^(-2·rango(p)·beta) + e^(-3·rango(p)·beta) + … + e^(-k·rango(p)·beta) + …
Scritta in forma compatta:
S = Σ(k=0 → ∞) e^(-k·rango(p)·beta)
Questa è una serie geometrica di ragione r = e^(-rango(p)·beta). Una serie geometrica converge — cioè dà un risultato finito — a condizione che la ragione sia minore di 1 in valore assoluto. Nel nostro caso r = e^(-rango(p)·beta): poiché rango(p) ≥ 1 e beta > 0, l’esponente è sempre negativo, quindi r è sempre compreso tra 0 e 1. La condizione è sempre soddisfatta, per qualsiasi valore di beta.
La somma converge quindi a:
1 / (1 – e^(-rango(p)·beta))

Abbiamo calcolato il contributo del singolo primo p, ma Z_B non è il contributo di un solo primo — è il catalogo di tutti gli stati di tutti i numeri contemporaneamente.

Come si passa dal singolo primo a tutto il sistema?
Ricordi l’esempio di prima? Lo stato {2,3,5} del numero 10 aveva peso e^(-1·beta) · e^(-2·beta) · e^(-3·beta) · 1 · 1 · 1 · … — un funzione zetaprodotto in cui ogni primo contribuisce in modo indipendente. Quando sommiamo su tutte le possibili combinazioni di tutti gli interruttori simultaneamente, il risultato è il prodotto dei contributi individuali — esattamente come il numero di combinazioni di due dadi indipendenti è il prodotto dei loro esiti possibili.

Quindi Z_B è il prodotto, su tutti i primi, della serie geometrica di ciascuno:
Z_B= 1/(1-e^(-rango(2)·beta)) · 1/(1-e^(-rango(3)·beta)) · 1/(1-e^(-rango(5)·beta)) · 1/(1-e^(-rango(7)·beta)) · …
Sostituendo i ranghi — rango(2)=1, rango(3)=2, rango(5)=3, rango(7)=4 — otteniamo:
Z_B = 1/(1-e^(-1·beta)) · 1/(1-e^(-2·beta)) · 1/(1-e^(-3·beta)) · 1/(1-e^(-4·beta)) · …

E qui emerge qualcosa di sorprendente. Gli esponenti sono 1, 2, 3, 4, 5, 6, 7… — tutti i numeri interi in sequenza, nessuno escluso. Non compaiono i valori dei primi — 2, 3, 5, 7, 11 — ma i loro ranghi, che sono semplicemente tutti i numeri naturali uno dopo l’altro.
Perché? Perché i ranghi dei primi sono per definizione 1, 2, 3, 4, 5… — ogni primo occupa esattamente un posto nella sequenza, senza buchi e senza ripetizioni. I primi entrano come etichette degli interruttori, ma quando li traduci in ranghi ottieni la sequenza completa dei naturali. La struttura matematica della funzione abbraccia tutti i numeri interi — non solo i primi.

Per il sistema fermionico il ragionamento è identico, ma ogni primo può comparire solo 0 o 1 volta. La serie geometrica si tronca al secondo termine:

Z_F = (1+e^(-1·beta)) · (1+e^(-2·beta)) · (1+e^(-3·beta)) · (1+e^(-4·beta)) · …

E qui arriva la grande sorpresa: tra Z_B e Z_F vale questa identità esatta, dimostrabile in tre righe di algebra:

Z_B(β) = Z_F(β) · Z_B(2β)

Vale qualunque valore di beta tu scelga — grande, piccolo, intermedio. Iterando all’infinito:

Z_B(β) = Z_F(β) · Z_F(2β) ·  Z_F(4β)· Z_F(8β)…

Il sistema bosonico si scompone in infiniti sistemi fermionici con beta via via raddoppiato. L’errore numerico è zero — non è un’approssimazione, è un’uguaglianza esatta. Non ce lo aspettavamo.

---

Il problema: il sistema bosonico esplode

Al crescere di n, il numero di stati bosonici cresce senza controllo. Per n=100 ci sono migliaia di partizioni, per n=1000 milioni. È un mare di Dirac — un oceano di stati in cui qualsiasi struttura si annega.

Ricordate il filo rosso? I bosoni rivelano i primi, ma a caro prezzo: li sommergono in un oceano di stati ripetitivi. Abbiamo cercato due strade per domare questa esplosione senza perdere la struttura. Entrambe portano allo stesso posto.

---

Prima strada: V* e Mann-Whitney-Wilcoxon

La prima soluzione non impone nulla dall’esterno — misura invece quanto ogni singolo stato è già fermionico di suo.
Per ogni stato assegniamo un punteggio di fermionicità — chiamiamolo V* — basato sulla dispersione dei ranghi dei suoi pezzi.

La logica è semplice: uno stato con tutti i pezzi uguali ha dispersione zero e è un bosone puro.
Lo stato {2,2,2,2} ha tutti i pezzi con rango 1 — nessuna dispersione. V* = 0, bosone puro.

Uno stato con pezzi ai ranghi più distanti possibile ha dispersione massima — è un fermione puro.
Lo stato {2,3,5,7} ha ranghi 1, 2, 3, 4 — massimamente dispersi. V* = 1, fermione puro.

Il test di Mann-Whitney-Wilcoxon è uno strumento statistico classico che misura quanto due insiemi di numeri siano diversi tra loro. Applicato ai ranghi di una partizione, produce un valore che misura esattamente questa dispersione. E qui emerge un risultato sorprendente: il valore medio atteso è sempre lo stesso numero — k(k+1)/2, dove k è il numero di addendi — qualunque sia lo stato e qualunque sia il numero. Un invariante universale, che non dipende da niente.

V* ci permette di costruire una versione corretta di Z_B — chiamiamola Z_B* — dove ogni stato viene pesato non solo per la sua energia ma anche per la sua fermionicità:

Z_B*= somma su tutti gli stati di e^(-energia·beta) · V(stato)*

Gli stati bosonici puri — V*=0 — vengono azzerati automaticamente. Gli stati più fermionici dominano.
Per il numero 10, Z_B* azzera completamente {2,2,2,2,2} e {5,5}.

Ma qui il filo rosso si intreccia in modo inatteso. Certi numeri primi — come l’11, il 17, il 23 — hanno come stato fermionico più naturale se stessi: {11}, {17}, {23}. Un singolo elemento ha sempre V*=0, perché non c’è dispersione con un solo pezzo. Quindi nella Z_B* questi primi vengono azzerati — si nascondono di nuovo, questa volta nella funzione corretta che avrebbe dovuto trovarli.

È il paradosso del filo rosso: i primi si rivelano nei bosoni, ma quando proviamo a usare la fermionicità per trovarli meglio, li perdiamo di nuovo. Stanno sempre sul confine tra i due sistemi — né completamente bosonici né completamente fermionici. Non si lasciano catturare da nessuno dei due.

---

Seconda strada: il cursore alfa — e la porta verso Riemann

La seconda soluzione impone una regola dall’esterno: ogni primo p può apparire al massimo ⌊p^α⌋ volte — cioè la parte intera di p elevato ad alfa.

Immagina un cursore che va da 0 a infinito:

alfa = 0 — ogni primo compare al massimo 1 volta. Fermioni puri. Il 4 e il 6 sono vuoti.
alfa = 1 — ogni primo p compare al massimo p volte. Il 2 al massimo 2 volte, il 3 al massimo 3 volte.
alfa infinito — nessun limite. Bosoni puri.

Alzando alfa da 0, a un certo punto il 4 e il 6 si risvegliano.
Il 6 si risveglia a alfa = log₃2 — perché la condizione è che 3^alfa raggiunga 2, e il valore di alfa per cui 3^alfa = 2 è per definizione log₃2 ≈ 0.631.
Il 4 si risveglia dopo, a alfa = log₂2 = 1 — perché serve che 2^alfa raggiunga 2.

Il valore log₃2 coincide con la dimensione di Hausdorff dell’insieme di Cantor ternario — il frattale che si costruisce togliendo il terzo centrale di un segmento all’infinito.
Non sappiamo se sia una coincidenza.

Ma il valore alfa = 1 è il punto cruciale. A alfa = 1, ogni primo p compare al massimo p volte — una regola auto-referenziale: ogni primo porta scritto in sé stesso il proprio limite.
Questo sistema si chiama sistema self-quantum.

---

La convergenza: il sistema self-quantum e Riemann

Ricordiamo dove siamo. Nel sistema bosonico, ogni primo p può comparire quante volte vuole — la serie geometrica è infinita. Nel sistema fermionico, ogni primo può comparire al massimo una volta — la serie si tronca al secondo termine. Il cursore alfa ci permette di stare in mezzo: ogni primo p può comparire al massimo ⌊p^α⌋ volte.

Per alfa = 1, ogni primo p può comparire al massimo p volte. Il 2 al massimo 2 volte, il 3 al massimo 3 volte, il 5 al massimo 5 volte. È una regola auto-referenziale: ogni primo porta scritto in sé stesso il proprio limite.

Cosa succede alla serie geometrica del primo p in questo sistema? Invece di sommare infiniti termini:

1 + e^(-rango(p)·beta) + e^(-2·rango(p)·beta) + e^(-3·rango(p)·beta) + … all’infinito

ci fermiamo al termine p-esimo:

1 + e^(-rango(p)·beta) + e^(-2·rango(p)·beta) + … + e^(-p·rango(p)·beta)

Il 2 si ferma al secondo termine. Il 3 si ferma al terzo. Il 5 si ferma al quinto. Ogni primo tronca la propria serie al termine che porta il suo stesso numero — non il suo rango, ma il suo valore.

La funzione di partizione del sistema self-quantum è quindi il prodotto di tutte queste serie troncate:

Z_S(β) = ∏_p (1 + e^(-rango(p)·β) + e^(-2·rango(p)·β) + … + e^(-p·rango(p)·β))

Ora facciamo qualcosa di preciso: sostituiamo e^(-rango(p)·beta) con p^(-s). Perché questa sostituzione? Perché e^(-rango(p)·beta) e p^(-s) sono entrambi modi di assegnare un peso al primo p — il primo usa il rango moltiplicato per beta, il secondo usa il logaritmo del valore moltiplicato per s. Sono due famiglie di pesi diversi, ma quando costruiamo il prodotto infinito su tutti i primi, i due risultati convergono allo stesso valore. Non è una identità algebrica — è una convergenza numerica, verificabile passo per passo.

Sotto questa sostituzione, la serie troncata del primo p diventa:

1 + p^(-s) + p^(-2s) + … + p^(-ps)

E il prodotto su tutti i primi diventa Z_S(s).

Ora confrontiamo con la funzione zeta di Riemann, che abbiamo già incontrato:

ζ(s) = ∏_p 1/(1 – p^(-s)) = ∏_p (1 + p^(-s) + p^(-2s) + p^(-3s) + … all’infinito)

La struttura è identica — stesso prodotto su tutti i primi, stessa serie geometrica — ma ζ(s) ha la serie infinita, mentre Z_S(s) la tronca al termine p-esimo.

La differenza è il termine troncato: p^(-ps) e tutti i termini successivi che stiamo buttando via. Quanto sono grandi questi termini? Per p=2: buttiamo via 2^(-2s) + 2^(-3s) + … che per s=2 vale circa 0.06 — piccolo ma non trascurabile. Per p=3: buttiamo via 3^(-3s) + … che per s=2 vale circa 0.004. Per p=5: buttiamo via 5^(-5s) + … che per s=2 vale circa 0.00003. I termini troncati decrescono rapidissimamente al crescere di p.

Il risultato è che Z_S(s) converge a ζ(s) con un errore inferiore al 2% per s=2, e inferiore allo 0.003% per s=5.
Non è un’approssimazione costruita a mano — emerge dalla regola auto-referenziale del sistema self-quantum.

---

Cosa abbiamo trovato, e cosa non sappiamo ancora

E qui arriva la conclusione del filo rosso.

Abbiamo percorso due strade completamente diverse, partendo dallo stesso punto — le somme di numeri primi.

La prima strada è quella della curvatura bosonica: seguendo il gradiente della media degli addendi bosonici, la traiettoria trova i numeri primi nel 56% dei passi contro il 27% atteso da Gauss.
La stima di Gauss emerge spontaneamente — senza cercarla.

La seconda strada è quella del cursore alfa: imponendo la regola auto-referenziale alfa=1, la funzione di partizione del sistema self-quantum converge alla funzione zeta di Riemann.
La funzione zeta emerge spontaneamente — senza cercarla.

Due strade, una conclusione. La struttura che Gauss intuì contando i primi a 15 anni, e che Riemann formalizzò nel 1859 legandola alla funzione zeta, emerge da sola — due volte, da direzioni diverse — in un framework costruito semplicemente chiedendo: in quanti modi posso scrivere un numero intero come somma di primi?

Non l’abbiamo cercata. È arrivata da sola.

E ancora oggi — con un milione di dollari in palio per chiunque dimostri l’Ipotesi di Riemann — nessuno ha completamente capito perché.

I risultati ci sembrano interessanti e non sappiamo se qualcuno ci abbia già pensato.

Un’identità esatta tra le due funzioni di partizione, verificata con errore zero. Una soglia critica — log₃2 — che coincide con la dimensione di Hausdorff dell’insieme di Cantor. La convergenza del sistema self-quantum alla funzione zeta di Riemann. Un invariante universale nei ranghi statistici — k(k+1)/2 — che non dipende da nessun parametro. La firma dei primi nella curvatura bosonica, che batte la stima di Gauss quasi due volte — mentre la curvatura fermionica li perde completamente. E un risultato aperto: V* azzera i bosoni puri ma azzera anche i numeri primi che si scrivono solo come se stessi — quelli che nel sistema fermionico sono invisibili. I primi stanno sul confine tra i due mondi, e non si lasciano catturare da nessuno dei due.

Nel frattempo, il gioco continua.

---

Questo framework è stato sviluppato da un chimico in collaborazione con Claude (Anthropic). Le idee sono umane; la formalizzazione matematica è stata un lavoro a quattro mani.