Termodinamica statistica della sostanza pura ideale e dei suoi stati di aggregazione

<a id="gesamt"></a>

Riassunto generale

Il modello 3 è un modello qualitativo, di meccanica statistica, di una sostanza ideale che ricava i tre stati di aggregazione — solido, liquido e gassoso — insieme alle loro transizioni da un'unica immagine di legame, quella dei legami molecolari. Il principio di fondo sta nel tener conto di tre serbatoi di energia per l'energia interna: l'energia di vibrazione nei legami molecolari, l'energia di traslazione delle molecole e l'energia cinetica cooperativa nella rotazione delle micelle, forma idealizzata dei flussi interni di materia che compensano le fluttuazioni interne di densità. Una grandezza fondamentale importante in questo contesto è il numero di coordinazione medio z̄ = z_max·x, che si stabilisce in modo autoconsistente attraverso un legame cooperativo a due stati (di tipo Ising, cooperatività K = J·z_max, degenerazione dello stato rotto G = z_max). Da esso seguono tre transizioni su una scala τ comune, calibrata mediante un solo parametro di energia ε: la fusione come soglia di fononi/rigidità (z̄ raggiunge la rigidità), l'ebollizione come collasso della connettività (percolazione) al punto critico e una transizione facoltativa solido(1)→solido(2) (FCC↔BCC) con il gradino di c_v privo di parametri di 1,0 k_B. Una transizione diventa del primo ordine non appena ln G > 2 + 4ε/K. Qui la degenerazione dello stato rotto G è il numero di contatti di vicinato che possono rompersi (G = z_max, il numero di coordinazione massimo), l'energia di legame ε è l'energia guadagnata per contatto di vicinato e la cooperatività K = J·z_max è il prodotto della forza di accoppiamento J, di tipo Ising, per z_max; la transizione diventa dunque del primo ordine tanto prima quanto maggiore è la molteplicità dei contatti G e quanto più forte è la cooperatività K rispetto all'energia di legame ε.

L'articolo conduce dai sette assiomi (capitolo 1) attraverso l'analisi delle contraddizioni (capitolo 2) e una descrizione per immagini con strisce a fumetti statiche (capitolo 3) fino all'impalcatura di meccanica statistica — legame cooperativo, funzioni di partizione e capacità termiche (capitoli 4–5) —, alla plausibilità matematica in 3D (capitolo 6) e infine alla calibrazione e all'equazione di stato (capitolo 7). Le idee sono integrate dalla scansione 3D completa di c_v della fusione del primo ordine (Δz̄, calore latente ΔU, spostamento del picco), da un modello a molle di Einstein per la rigidità fononica θ(reticolo) e dalla linea di trasformazione solido(1)→solido(2) via Clausius–Clapeyron. θ(reticolo) viene poi raffinato in una densità di stati fononici di tipo Debye; la trasformazione solido(1)→solido(2) viene inoltre introdotta come linea nel piano (τ, p) (calibrata sul ferro) e viene aggiunto un ramo Redlich–Kwong separato accanto all'impostazione predefinita di van der Waals. θ(reticolo) è inoltre verificato su un reticolo di Lennard-Jones completo (FCC stabile, BCC ideale instabile a taglio) e la topologia solido(1)→solido(2) è verificata sul titanio (seconda conferma) e sull'argon (controllo negativo). Il capitolo 8 colloca il modello fra undici concezioni consolidate della termodinamica statica. L'appendice documenta le fonti bibliografiche (ordinate per autore) e il supporto dell'IA (quota d'innovazione stimata, costo stimato e durata per passo di sviluppo); gli ausili raccolgono il modo di lettura, i punti critici rimasti e le prospettive.

Enunciato centrale. Un'unica misura cooperativa del legame, z̄, tiene insieme tre tradizioni altrimenti separate — l'immagine del gas reticolare / di Ising (transizione cooperativa a due stati), l'immagine della percolazione (connettività → ebollizione) e l'immagine dei fononi di Einstein/Debye (rigidità → fusione) — in un quadro delle fasi coerente e qualitativamente completo, deliberatamente entro un'immagine ideale di tipo van der Waals, potendosi naturalmente ricavare l'edificio teorico da un'immagine dinamica facilmente immaginabile. Non l'esattezza, ma il carattere figurato è la vera forza del modello e ne caratterizza il valore scientifico come modello.

Indice

• Riassunto generale
• 1 — Gli assiomi con le loro parole chiave (consolidati)
• 2 — Analisi delle contraddizioni e delle lacune
• 3 — Descrizioni per immagini
• 3.1 Le tre fasi singole
• 3.2 Transizioni di fase
• 3.3 Perché la velocità del suono nei liquidi cala con T
• 3.4 Come nasce il moto browniano — con una nozione affinata di micella
• 3.5 Perché le reazioni agli strati limite delle micelle diventano più rapide
• 3.6 La finestra liquida tridimensionale
• 4 — Modello di meccanica statistica
• 4.1 Gradi di libertà e variabili
• 4.2 Legame cooperativo a due stati (di tipo Ising, parametro J)
• 4.3 Bilancio energetico — con termine cooperativo
• 5 — Funzioni di partizione e capacità termiche
• 5.1 Funzione di partizione fattorizzata
• 5.2 Capacità termica per fase (analitica)
• 5.3 Curva numerica c_v(τ) **con** termine cooperativo — calore latente quantificato
• 5.4 Soglia di percolazione di connettività come risultato (confermata invariata) - 5.4.a Due soglie distinte nettamente separate
• 5.5 Ordine della transizione di fusione
• 5.6 Studio di dimensione finita della soglia di percolazione
• 5.7 Scansione 3D completa di c_v(τ) della fusione del primo ordine
• 6 — Plausibilità matematica dell'interpretazione 3D
• 6.1 Variazione del grado di libertà di rotazione 2D → 3D
• 6.2 Influenza dei diversi reticoli cristallini — la connettività come tratto distintivo
• 6.3 La transizione aggiuntiva solido(1) → solido(2)
• 6.4 Considerazioni: l'anomalia acustica dell'acqua liquida
• 6.5 Alternativa: campi di vortici frattali invece della rotazione delle micelle
• 6.6 Rigidità fononica microscopica θ(reticolo): un modello a molle di Einstein
• 6.7 Stima matematica della metastabilità (2D triangolato ↔ quadrato, 16b)
• 6.8 Raffinamento di tipo Debye di θ(reticolo)
• 6.9 θ(reticolo) da un reticolo di Lennard-Jones completo
• 7 — Calibrazione ed equazione di stato
• 7.1 Metodo di calibrazione (un parametro di energia ε) — adattato
• 7.2 Tabella di calibrazione
• 7.3 Asse densità/volume: equazione di stato di van der Waals con nocciolo duro
• 7.4 Calibrazione quantitativa l/g: a → ε, b → D²
• 7.5 Superficie (τ, v) comune con punto triplo
• 7.6 Calore latente a fronte delle entalpie di fusione
• 7.7 Unificazione delle scale di energia: ε ↔ ε_lg
• 7.8 Punto triplo senza parametri a partire da z̄_c
• 7.9 Calcolo di controllo 3D: una finestra liquida netta su un'unica scala ε
• 7.10 Vera scansione di percolazione di legami in 3D (FCC/BCC) e sensibilità del bordo della finestra
• 7.11 Vera scansione in τ della cinetica di legame in 3D: τ_b(3D) e la linea di fusione su una scala ε
• 7.12 La linea di trasformazione solido(1)→solido(2) nella superficie (τ, v)
• 7.13 Linea di trasformazione solido(1)→solido(2) nel diagramma (τ, p), calibrata sul ferro
• 7.14 Redlich–Kwong come ramo separato dell'immagine ideale
• 7.15 solido(1)→solido(2) su altri polimorfi: titanio e controllo con gas nobile
• 8 — Confronto bibliografico: concezioni della termodinamica statica
• 8.1 Tabella d'insieme
• 8.2 I modelli uno per uno
• 8.3 Somiglianze e differenze rispetto al modello 3
• 9 — Appendice
• 9.1 Bibliografia
• 9.2 Supporto dell'IA
• 9.3 Glossario di abbreviazioni e simboli
• 9.4 Script di calcolo e copertura dei test

<a id="kap1"></a>

1 — Gli assiomi con le loro parole chiave (consolidati)

Riassunto. Sette assiomi dotati di parole chiave fissano la sostanza ideale: numero di particelle conservato, moto puramente classico, un legame cooperativo a due stati e il numero di coordinazione medio z̄ come parametro d'ordine. Essi reggono tutte le derivazioni che seguono; nel testo tedesco di riferimento una filigrana di omoglifi segna le parole chiave (chiave in fondo al testo).

Gli assiomi precisati:

1. VOLO LIBERO (assioma aggiuntivo). Fra due molecole che si sfiorano senza essere legate non agiscono né attrito né forze elettrostatiche. Quantità di moto ed energia cambiano solo per la repulsione dura del NOCCIOLO DURO nell'urto e per la formazione e la rottura di un LEGAME. L'assioma conta soltanto dal punto di vista della meccanica classica e del bilancio energetico: garantisce che non esistano canali dissipativi a livello molecolare — l'energia del moto non va mai «perduta», passa soltanto nei tre serbatoi della TRIADE o in energia di legame. (Attrito e viscosità macroscopici nascono in modo emergente dal gioco del quindici e dalla FRATTURA, non come forza fondamentale microscopica.)
2. MOLTITUDINE. Una sostanza ideale è fatta di moltissime molecole uguali; le proprietà macroscopiche sono descritte statisticamente nel limite termodinamico (N → ∞).
3. LEGAME. Due molecole geometricamente vicine possono portare un legame. Un legame intatto abbassa l'energia della quantità finita ε; la sua rottura costa esattamente ε. Il legame è un sistema a due stati {intatto, rotto}.
4. TRIADE. L'energia interna si ripartisce classicamente su tre serbatoi: (a) traslazione delle molecole, (b) vibrazione nei legami intatti (una molla armonica per legame), (c) rotazione di blocchi molecolari rigidi (grado di libertà rotazionale collettivo delle regioni connesse).
5. NOCCIOLO DURO. Se la distanza molecolare scende sotto un valore minimo σ, agisce una repulsione dura (volume proprio). Insieme al LEGAME ne risulta un potenziale di coppia efficace (guscio duro più buca di legame finita).
6. CONSERVAZIONE. La quantità di moto totale e l'energia totale si conservano. Una singola molecola è energeticamente neutra — non immagazzina energia in sé; l'energia vive nei legami (LEGAME) e nei gradi di libertà cinetici (TRIADE).
7. FRATTURA. I legami si disfano e si formano secondo il principio fessurazione/risanamento: più legami vicini mancano → maggiore probabilità di rottura; più vicini legati → maggiore probabilità di risanamento. La cinetica è cooperativa.
8. SVELTO. Rottura e risanamento avvengono rapidamente e con poca attivazione; il sottosistema dei legami si trova a ogni temperatura in equilibrio locale e può essere separato adiabaticamente dal moto molecolare (fondamento della funzione di partizione fattorizzata).

Atteggiamento di modellazione. La coordinazione geometrica si conserva in modo idealizzato; cambia solo il numero di coordinazione z per effetto delle rotture. Il moto nasce come nel gioco del quindici: le cavità (dal volume proprio del NOCCIOLO DURO e dalle rotture della FRATTURA) migrano, i blocchi avanzano. L'immagine deliberatamente «sbagliata» di cristalliti che ruotano e migrano è mantenuta come base intuitiva.

<a id="kap2"></a>

2 — Analisi delle contraddizioni e delle lacune

Riassunto. L'analisi esamina sistematicamente il modello quanto a contraddizioni interne e lacune e aggiorna lo stato del lavoro; le debolezze riconosciute vengono nominate e risolte nei capitoli seguenti oppure segnalate come limite onesto del modello.

Elenco dei punti critici

1. Doppio conteggio della traslazione (assioma 3a contro 3c). Il testo grezzo nomina la traslazione due volte (traslazione molecolare e traslazione di blocco). Due serbatoi di traslazione indipendenti non esisterebbero. Risoluzione: la traslazione di blocco è solo un modo collettivo delle traslazioni molecolari già contate; genuinamente nuovo è unicamente il grado di libertà di rotazione del blocco. Nel bilancio dei gradi di libertà (cap. 4) la traslazione è contata esattamente una volta.
2. «Energeticamente neutra» (assioma 5) contro «i legami assorbono energia» (assioma 2/3b). Contraddizione apparente. Risoluzione: la neutralità riguarda la molecola come corpo; l'energia risiede fra le molecole (nel LEGAME) e in modi cinetici collettivi, mai in un serbatoio intramolecolare. Molecola ≠ legame — così la contraddizione è risolta.
3. Attrazione assente (l'assioma 4 nomina solo la repulsione). Senza attrazione non vi sarebbe fase condensata. Risoluzione: l'attrazione è fornita dal LEGAME (buca di profondità ε), la repulsione dal NOCCIOLO DURO (guscio duro in σ). Insieme, un potenziale di coppia efficace.
4. Retroazione senza termine opposto (assioma 6). «Più vicini mancanti → più rotture» è di per sé una retroazione positiva (avanzamento della frattura fino alla disgregazione totale). Risoluzione: il termine di risanamento (più vicini → più risanamento) più il fattore di Boltzmann dell'energia di legame finita ε (LEGAME) bilanciano la frattura e producono una transizione definita invece di una fuga.
5. Separazione delle scale temporali non esplicitata (assioma 7). «Rapida e a bassa attivazione» diventa utilizzabile solo se se ne ricava la separazione adiabatica fra dinamica dei legami e dinamica di traslazione. Risoluzione: posta esplicitamente come presupposto della fattorizzazione Z = Z_trans·Z_vib·Z_rot·Z_bond.
6. Geometria assente e meccanismo del moto assente. Gli assiomi non nominano né il reticolo di coordinazione né ciò che rende possibile il flusso. Risoluzione (come ipotesi aggiuntiva esplicitamente segnalata): una geometria di coordinazione di riferimento fissa con occupazione variabile dei legami; come caso 2D canonico la coordinazione triangolare (z_max = 6). Mobilità attraverso il meccanismo delle cavità del gioco del quindici. Questi due punti sono le uniche vere aggiunte oltre i sette assiomi e sono segnalati nel capitolo 6 come questione strategica.
7. Lacuna latente di dissipazione. La versione 07 lasciava aperto se molecole non legate perdano energia mentre si sfiorano. Risolta dal VOLO LIBERO: non vi sono attrito né elettrostatica microscopici; tutti i flussi di energia e quantità di moto passano esclusivamente per gli urti del NOCCIOLO DURO e la cinetica del LEGAME. Con ciò il bilancio energetico (cap. 4.3) è chiuso e la statistica di Maxwell–Boltzmann esatta nel limite gassoso (gas ideale libero fra urti elastici).
8. Incoerenza energetica del parametro d'ordine. La versione 08 usava il termine cooperativo −(K/2)x² nell'energia libera g(x) (per generare x(τ)), ma non nell'energia interna U/N. Ciò è termodinamicamente incoerente: lo stesso termine di campo molecolare deve comparire in entrambe. Risoluzione (cap. 4.3): −(K/2)x² viene riportato in U/N; il calore latente ne risulta misurabilmente accresciuto (cap. 5.3), senza modificare i plateau solido e gassoso.
9. Calibrazione l/g solo qualitativa. La versione 09 aveva l'asse delle densità (equazione di stato di van der Waals) solo in forma ridotta. Risoluzione (cap. 7.4): la linea di ebollizione assoluta è calcolata tramite il punto critico (T_c, p_c) e T_b(p_atm) è confrontato con i valori misurati; lo scarto sistematico (regola di Guldberg) è nominato.
10. Assi separati senza punto triplo. Gli assi di fusione ed ebollizione erano calcoli separati. Risoluzione (cap. 7.5): entrambi sono riuniti in un piano (τ, v), il punto triplo è costruito e il punto critico segnato.

<a id="kap3"></a>

3 — Descrizioni per immagini

Riassunto. Il capitolo 3 spiega in modo intuitivo il solido, il liquido e il gassoso nonché le loro transizioni (con strisce a fumetti statiche): le tre fasi singole (3.1), le transizioni di fase ricristallizzazione/fusione/ebollizione/sublimazione insieme al punto triplo e a quello critico (3.2), l'anomalia acustica (3.3), il moto browniano e la nozione affinata di micella (3.4–3.5) nonché la finestra liquida tridimensionale (3.6). Sulla base di queste descrizioni sono stati realizzati schemi (strisce a fumetti) che devono aiutare a raffigurarsi davanti all'occhio interiore i rapporti molecolari dinamici nelle diverse fasi. Questa raffigurabilità è il vero nucleo del modello. Esso acquista valore qualitativo e scientifico quando la meccanica statistica permette di mostrare in modo plausibile che il modello rende prevedibili tendenze e osservazioni empiriche qualitative anche senza grande apparato matematico.

<a id="kap3-1"></a>

3.1 Le tre fasi singole

Solido. Quasi tutti i vicini sono legati fra loro da legami molecolari (LEGAME); la media dei numeri di coordinazione z delle molecole è vicina al massimo geometrico. La rete molecolare è un unico blocco rigido senza cavità (NOCCIOLO DURO, impacchettamento denso) — un gioco del quindici interno come nel liquido non è possibile nella versione compatta. Unico movimento: il tremolio attorno alle posizioni di riposo (vibrazione di legame, TRIADE-b), determinato dai legami molecolari con i vicini. Quantità di moto ed energia oscillano elasticamente nella rete (CONSERVAZIONE, VOLO LIBERO: senza perdite). Nel caso limite (temperatura molto bassa, impacchettamento ideale) si ha lo stato ideale di un cristallo perfetto.

Liquido. La FRATTURA cooperativa ha attraversato la rete: mancano legami a chiazze, la media dei numeri di coordinazione z delle molecole è scesa, il blocco rigido si è internamente sgretolato in molti blocchi più piccoli. Dal monolito è nato un cumulo di detriti. Questi blocchi traslano e ruotano (TRIADE-c). Fra di essi si sono aperte cavità attraverso cui si spostano l'uno accanto all'altro — il gioco del quindici del moto browniano è attivo. La geometria di coordinazione locale si conserva in modo idealizzato; è cambiato solo quanti legami esistono realmente e si formano.

Gas. La FRATTURA ha vinto: i legami sono rotti quasi ovunque (z → 0). Nel caso estremo nessuna molecola ha più un legame molecolare con un'altra. I blocchi si sono disgregati in singole molecole. Volano libere, accoppiate solo dagli urti del NOCCIOLO DURO (con VOLO LIBERO nel mezzo); le velocità seguono Maxwell–Boltzmann (CONSERVAZIONE). Vibrazione di legame e rotazione di blocco si sono spente — traslazione pura (TRIADE-a). Nel caso limite (temperatura alta, pressione bassa) si ha lo stato ideale di un gas ideale.

<a id="kap3-2"></a>

3.2 Transizioni di fase

Ricristallizzazione (solido(1)/solido(2)). Ancora sotto la fusione, il solido può cambiare il proprio reticolo senza perdere la rigidità: due impacchettamenti ideali con numeri di coordinazione diversi (FCC z_max = 12 ↔ BCC z_max = 8) hanno entalpie libere μ(τ, p) diverse; a μ_FCC = μ_BCC la geometria di coordinazione salta di colpo (non solo l'occupazione z̄ come nella fusione). Diventa visibile un piccolo gradino di c_v (Δc_v = variazione di z_max/4 = 1,0 k_B) più un debole calore latente di trasformazione. Affermazione dell'IA: è concepibile solo in 3D (connettività discreta, quantizzata; cap. 6.3). Nel caso bidimensionale la vera ricristallizzazione, cioè la transizione fra due fasi solide, è più facile da immaginare. Nel piano esistono solo tre tassellature regolari senza buchi. Si può ricoprire interamente una superficie (tassellarla) con triangoli equilateri, con quadrati o con esagoni regolari. Se in questi reticoli gli spigoli rappresentano i legami molecolari e i vertici le molecole, allora è concepibile una transizione di fase solido(triangolare) → solido(quadrato). La molecola nella tassellatura triangolare avrebbe ancora sei legami verso le molecole vicine, mentre quella nel reticolo quadrato ne avrebbe solo quattro.

Fusione (s/l). Al crescere della temperatura la rete si lacera in abbastanza punti perché il blocco rigido complessivo perda la propria rigidità (la FRATTURA cooperativa scende sotto la soglia di rigidità di Maxwell z̄ = 2d — la rete diventa molle). Si formano micelle termodinamiche, cioè regioni che non sono più collegate da legami alle regioni vicine più grandi. Dal punto di vista molecolare il monolito diventa un cumulo di detriti e acquista una capacità di scorrimento limitata. Le prime cavità consentono lo spostamento di blocchi molecolari, come nel gioco del quindici. I legami mancanti verso il vicinato consentono la rotazione interna nelle micelle. La rotazione delle micelle (TRIADE-c) si avvia. La rottura di molti legami (LEGAME) divora calore latente.

Ebollizione (l/g). Un ulteriore riscaldamento fa percolare le cavità stesse: esse si uniscono in uno spazio vuoto continuo, in cui i blocchi si dissolvono e singole molecole volano via. L'ebollizione è il passaggio da «cavità nella materia» a «materia in una cavità».

Sublimazione (s/g). A bassa pressione manca la zona liquida intermedia: la rete cede singole molecole direttamente allo spazio vuoto, senza che si giunga alla formazione di blocchi come nelle micelle. Ci si può chiedere qui quale influenza possano avere sulla formazione o non formazione di micelle termodinamiche gli aspetti geometrici (molecole lunghe, molecole che si intrecciano fra loro), gli aspetti elettrostatici (allineamento per momenti magnetici o dipolari) e i diversi tipi di legame molecolare (van der Waals, interazioni fra sistemi mesomeri, legami a idrogeno).

Punto triplo. Proprio lì coincidono due soglie: la rigidità di Maxwell (solido ↔ liquido, perdita della rigidità) e la percolazione delle cavità (liquido ↔ gas) intervengono allo stesso (τ, p). Tutte e tre le immagini (blocco rigido [monolito], blocchi che scorrono [cumulo di detriti], molecole libere [nube di polvere]) coesistono.

Punto critico. Al crescere della pressione il salto liquido→gas diventa sempre più piccolo; la frazione di vuoto cambia in modo sempre più continuo. Al punto critico la differenza scompare del tutto — «blocchi densi con poco vuoto» e «molto vuoto con poche molecole» passano con continuità l'uno nell'altro. («Dal punto di vista del liquido il bicchiere è mezzo vuoto, da quello del gas mezzo pieno.»)

Ricristallizzazione nell'immagine 2D. Il fumetto mostra ora la ricristallizzazione come caso 2D, uno sotto l'altro: in basso il reticolo denso triangolato (coordinazione z = 6, l'ideale 2D a energia minore), in alto il reticolo quadrato lasco (z = 4). La transizione triangolato → quadrato accresce il volume (lavoro contro la pressione) e abbassa il numero di legami. Si vede che il numero di coordinazione cambia mentre la rigidità di un reticolo cristallino permane.

Previsione qualitativa di transizioni di energia e capacità termiche. Dalla sola connettività z si possono stimare entrambe le cifre caratteristiche: (i) l'energia di legame per particella è U/N ≈ −(ε/2)·z (ogni legame conta per metà), dunque tanto minore quanto maggiore è z; (ii) il salto di c_v al cambio di reticolo segue la regola del modello Δc_v = variazione di z_max/4. Per il caso 2D triangolato (z = 6) ↔ quadrato (z = 4) ne risulta una differenza di energia di ΔU ≈ ε per particella e un gradino di c_v di Δc_v = (6 − 4)/4 = 0,5 k_B — metà del valore 3D (1,0 k_B). Regola mnemonica: più denso = a energia minore, ma più costoso nella transizione (calore latente maggiore, gradino di c_v maggiore), perché nell'allargamento si liberano più legami e più gradi di libertà congelati.

La struttura del reticolo determina la metastabilità? In parte. La struttura del reticolo fissa l'ordine delle energie — quale forma è lo stato fondamentale (triangolato, μ più basso) e quale quello eccitato (quadrato). Se la forma lasca sopravviva come metastabile a bassa temperatura non lo decide da sola: nel modello rigorosamente ideale (senza barriera) la sostanza assume a ogni (τ, p) immediatamente la fase con il μ più basso — lì la metastabilità non esiste. Essa nasce solo da una energia di attivazione finita (barriera di nucleazione) che frena cineticamente la ristrutturazione. La geometria del reticolo fornisce dunque la gerarchia termodinamica, la barriera la persistenza cinetica — lo stesso ingrediente che spiega anche la fase soprassatura (vedi sotto; quantitativamente cap. 6.7).

Fase soprassatura — energia di attivazione. La figura aggiuntiva (3.2, riga 2) mostra una fase gassosa soprassatura sopra un liquido: accanto a molecole singole, piccoli gruppi di 2–3 particelle aggregate (barre grigie) che volano come un tutto. Il modello spiega la soprassaturazione con una energia di attivazione apprezzabile per i legami molecolari: i piccoli aggregati (2–3) hanno un rapporto bordo/superficie sfavorevole, i loro legami di bordo sono incompleti — solo un nucleo critico cresce spontaneamente. Finché non nasce un nucleo sopracritico, il gas resta soprassaturo, benché μ_gas > μ_liquido consentirebbe già la condensazione. Proprio questa barriera manca nell'immagine ideale di equilibrio e viene qui nominata come contributo autonomo situato «fuori dall'immagine ideale» (stima quantitativa al cap. 6.7).

<a id="kap3-3"></a>

3.3 Perché la velocità del suono nei liquidi cala con T

Il suono è un'onda elastica: v = √(K_eff/ρ). K_eff nasce esclusivamente dai legami intatti (LEGAME) — essi sono le molle. Al crescere di T si rompono cooperativamente sempre più legami (FRATTURA); la rete di molle si assottiglia e K_eff cala più ripidamente della densità ρ. La rottura domina il numeratore, non il denominatore — perciò v cala. Interpretazione molecolare: la trasmissione del suono lungo i legami molecolari è più rapida della trasmissione meccanica per urti. Più fratture compaiono in un liquido al crescere della temperatura, tanto più spesso il suono deve essere trasportato da urti meccanici anziché da legami molecolari. Di conseguenza la velocità del suono diventa più lenta al salire della temperatura. (L'anomalia dell'acqua — un massimo acustico verso 74 °C — è trattata al cap. 6.4.)

<a id="kap3-4"></a>

3.4 Come nasce il moto browniano — con una nozione affinata di micella

In senso termodinamico-statistico una micella non è qui un aggregato di tensioattivi ma un blocco cinto dalla FRATTURA: i legami mancanti (i non-legami al confine del blocco) formano il «guscio» attorno a un cristallite racchiuso (il nucleo legato integralmente). Questo nucleo può compiere un moto di rotazione cooperativo interno (TRIADE-c) — il blocco ruota come un tutto mentre il suo guscio di confine lo disaccoppia dai vicini.

Il moto browniano nasce allora così: mediante lo scambio di legami molecolari (FRATTURA, SVELTO) singole molecole possono lasciare il guscio di una micella e passare in un'altra micella. Ogni cambio di questo tipo trasmette quantità di moto (CONSERVAZIONE; senza perdite grazie al VOLO LIBERO) e altera bruscamente forma e momento angolare di entrambe le micelle coinvolte. Una grande particella di prova (granello di polvere) è lambita da molte di queste riorganizzazioni micellari non correlate; la somma degli impulsi è nulla in media, ma ha varianza finita → cammino casuale. Più alta è T, tanto più frequenti sono gli eventi di FRATTURA e i cambi di micella, e tanto più forte il tremolio. Il moto browniano è così l'ombra macroscopica, ossia la conseguenza, del processo microscopico di scambio fra micelle.

<a id="kap3-5"></a>

3.5 Perché le reazioni agli strati limite delle micelle diventano più rapide

Premessa imbarazzante. L'argomentazione che segue è una deduzione dall'immagine e, per quanto ne so, finora non è stata né dimostrata né studiata empiricamente — devo però ammettere di non essermi preso la briga di documentarmi a fondo. L'unico indizio di novità è che i sistemi di IA non menzionano questa conclusione. Questa tesi è un esempio di come le rappresentazioni per immagini possano condurre a nuove domande.

In solventi moderatamente solventi i reagenti si depositano di preferenza sui gusci di confine delle micelle — nel modello proprio i confini frattura/blocco (FRATTURA). Due effetti accelerano la reazione: (i) la riduzione dimensionale — i reagenti non si cercano più nell'intero volume ma lungo un'interfaccia quasi bidimensionale; una ricerca 2D coglie statisticamente il bersaglio prima di una 3D. (ii) La concentrazione — il guscio di confine raccoglie i reagenti, la loro densità locale cresce. Per immagini, i confini di frattura agiscono come rotaie di guida che portano i partner l'uno verso l'altro.

L'effetto dimensionale ha un nucleo rigoroso: il cammino casuale (Pólya). Che una ricerca bidimensionale colga statisticamente meglio di una tridimensionale non discende solo dall'intuizione ma da un teorema del calcolo delle probabilità. George Pólya (1887–1985) dimostrò nel 1921 per il cammino casuale su un reticolo: in una e in due dimensioni il camminatore ritorna con probabilità 1 (dunque con certezza) infinite volte in ogni punto, sul reticolo piano ℤ² come sulla retta ℤ. A partire da tre dimensioni ciò non vale più: la probabilità di ritorno è strettamente minore di 1 (sul reticolo cubico semplice ℤ³ solo circa 0,34). Per immagini, secondo Pólya: «Un uomo ubriaco ritrova sicuramente la via di casa, un pesce ubriaco invece no.» Trasferito a due reagenti che devono incontrarsi più volte per un urto riuscito, ciò significa: in un solvente molto buono con cammino casuale tridimensionale libero, partner una volta separati si allontanano per sempre con probabilità finita — i tentativi ripetuti di reazione falliscono e la resa resta bassa per questa ragione puramente geometrico-statistica. A uno strato limite micellare (quasi) bidimensionale, invece, il reincontro è certo (probabilità 1): i partner si ritrovano ancora e ancora finché la reazione riesce. Lo strato limite non solo concentra, ma mette al tempo stesso a disposizione la dimensione più favorevole ai reincontri. L'affermazione (deliberatamente acuminata) suona: le rese nei solventi molto buoni tendono a essere peggiori per ragioni di cammino casuale.

Riferimento: G. Pólya (1921), Über eine Aufgabe der Wahrscheinlichkeitsrechnung betreffend die Irrfahrt im Straßennetz, Math. Ann. 84, 149–160 — citato qui secondo M. Eisermann, Höhere Mathematik 3 (Univ. di Stoccarda), esercizio W-2×2. (Originale consultato il 5.9.2026).

<a id="kap3-6"></a>

3.6 La finestra liquida tridimensionale

Diversamente dalle figure di 3.1–3.5, questa non mostra alcuna scena molecolare ma il risultato del calcolo dei cap. 7.9/7.10 come schema su un asse τ (in unità dell'energia di legame ε). L'obiettivo era rendere visibile perché il reticolo 2D isolato ammette a stento un dominio per una fase liquida, mentre col reticolo 3D ciò è possibile senza difficoltà, valendo in entrambi i casi la stessa regola di rigidità (z̄_c = 2d).

Immagine ① — 2D (nessuna finestra). Sull'asse τ, τ_m = 1,20 e τ_triple = 1,16 si sovrappongono praticamente (τ_triple sotto τ_m), con τ_c = 1,50 poco sopra: il dominio della fase liquida è compresso in una fessura — il solido passa quasi direttamente al gas (vicino alla sublimazione).

Immagine ② — 3D (finestra netta). Lo stesso asse, ma dilatato: τ_m = 1,20 < τ_triple = 2,83 < τ_c = 3,00. Fra τ_triple e τ_c sta il campo liquido autonomo (largo verso il basso, τ_triple − τ_m = 1,63; stretto verso l'alto, τ_c − τ_triple = 0,18). Il sottile bordo superiore è segnato come tale.

Immagine ③ — la causa. z̄_c = 2d cresce con la dimensione (4 → 6), perciò v_solid cala (0,50 → 0,667 in volumi ridotti, cioè il solido diventa «più lasco» rispetto al liquido), e τ_c = z_max/4 raddoppia (1,50 → 3,00). Il punto critico si allontana ben sopra il punto di fusione. Solo così nasce spazio per la fase liquida. È un puro effetto dimensionale.

<a id="kap4"></a>

4 — Modello di meccanica statistica (invariato dalla versione 09)

Riassunto

Il capitolo 4 imposta l'impalcatura di meccanica statistica del modello in 2D (reticolo triangolare, z_max = 6):

• 4.1 Gradi di libertà e variabili: la grandezza d'ordine è la frazione di legami x ∈ [0,1] con il numero di coordinazione medio z̄ = z_max·x; per blocco ≥ 2 si aggiunge un grado di libertà di rotazione, e la temperatura τ si conta in unità dell'energia di legame ε.
• 4.2 Legame cooperativo a due stati: l'equazione di legame autoconsistente x = 1/(1 + G·exp[−(ε + K·x)/τ]) (K = J·z_max, G = z_max) segue come condizione di stazionarietà dell'energia libera di Bragg–Williams; il suo minimo globale fornisce il ramo fisico, e il cap. 5.5 decide se esso decorra con continuità o con un salto (primo ordine).
• 4.3 Bilancio energetico con termine cooperativo: U/N = τ·[1 + z̄/4 + ½ n_blk] − (ε/2)·z̄ − (K/2)·x² (2D) unisce traslazione, vibrazione di legame e rotazione di blocco, l'energia di configurazione dei legami e la nuova energia cooperativa di campo molecolare −(K/2)x². La capacità termica c_v associata riceve il termine aggiuntivo positivo −K·x·ẋ, che acuisce il picco di fusione ma lascia intatti il plateau solido e il ramo gassoso.

<a id="kap4-1"></a>

4.1 Gradi di libertà e variabili

• N molecole su una geometria di coordinazione di riferimento fissa. Caso 2D canonico: reticolo triangolare, z_max = 6, un GDL di rotazione per blocco (MOLTITUDINE, LEGAME). Questo capitolo vale esplicitamente per l'interpretazione 2D; la generalizzazione 3D sta nel cap. 6.
• Variabile di legame s_ij ∈ {0, 1} per coppia geometrica di vicini (1 = intatto).
• Numero di coordinazione z_i = Σ_j s_ij; media z̄ = 2 N_bond / N — il parametro d'ordine. Si ha z̄ = z_max · x con la frazione di legami x = p_b ∈ [0, 1].
• Unità di energia τ ≡ k_B T (classicamente, ½τ per grado di libertà quadratico). Tutte le energie in unità dell'energia di legame ε.

<a id="kap4-2"></a>

4.2 Legame cooperativo a due stati (di tipo Ising, parametro J)

Un legame isolato (energia −ε intatto, 0 rotto) ha, nell'equilibrio di Maxwell–Boltzmann, p_b⁰(τ) = 1/(1 + e^(−ε/τ)). La FRATTURA cooperativa (assioma 6) sposta la stabilità efficace del legame di un campo molecolare proporzionale alla frazione di legami x. In approssimazione di campo molecolare (Ising) l'equazione di legame autoconsistente diventa:

x = 1 / (1 + G · exp[−(ε + K·x) / τ]), K = J·z_max, z̄ = z_max·x.

Qui:

• K = J·z_max è la forza di accoppiamento cooperativa (retroazione della FRATTURA: più legami vicini → maggiore stabilità efficace del legame). J > 0 acuisce la transizione; oltre un J critico z̄(τ) diventa discontinuo — questo è il calore latente della fusione. Il caso puro a due stati è il limite J = 0 (nessun campo molecolare, crossover continuo).
• G > 1 è la degenerazione dello stato rotto: una molecola liberata guadagna spazio delle fasi (libertà di traslazione). G spinge la rottura dei legami ad alto τ e fa sì che x → 1/(1+G) invece di restare fermo a ½ (caso limite del gas). Canonicamente G = z_max.

Questa equazione segue come condizione di stazionarietà di un'energia libera di Bragg–Williams per molecola (n_slots = z_max/2 posti di legame):

g(x)/n_slots = −ε·x − (K/2)·x² + τ·[x ln x + (1−x) ln(1−x)] − τ·ln(G)·(1−x).

Il termine −(K/2)x² è l'accoppiamento cooperativo (ferromagnetico), l'ultimo termine la degenerazione di rottura. Il minimo globale di g su x fornisce il ramo fisico; dove due minimi divengono ugualmente profondi si trova la transizione del primo ordine. SVELTO (assioma 7) giustifica che x(τ) assuma a ogni temperatura il valore di equilibrio.

<a id="kap4-3"></a>

4.3 Bilancio energetico — con termine cooperativo

Il termine di campo molecolare −(K/2)x² è un'energia genuina (indipendente da τ) e deve comparire in U/N. Il bilancio coerente recita dunque (2D):

U/N = τ·[1 + z̄/4 + ½ n_blk] − (ε/2)·z̄ − (K/2)·x².

• τ·[…]: traslazione (1) + vibrazione di legame (z̄/4) + rotazione di blocco (½ n_blk) — TRIADE.
• −(ε/2)·z̄: energia di configurazione dei legami — LEGAME.
• −(K/2)·x²: energia cooperativa di campo molecolare — FRATTURA. Nuovo in U/N.

La capacità termica associata c_v = d(U/N)/dτ (in unità k_B) riceve un contributo aggiuntivo:

c_v = [1 + z̄/4 + ½ n_blk] + τ·[¼ ż̄ + ½ ṅ_blk] − (ε/2)·ż̄ − K·x·ẋ, (ẋ = dx/dτ)

dove il nuovo termine −K·x·ẋ è positivo (ẋ < 0) e agisce esattamente là dove molti legami si rompono contemporaneamente. Esso acuisce il picco di fusione senza modificare il plateau solido (ẋ ≈ 0) né il ramo gassoso.

<a id="kap5"></a>

5 — Funzioni di partizione e capacità termiche

Riassunto

Il capitolo 5 calcola funzioni di partizione e capacità termiche e fornisce i dati del modello per la calibrazione (cap. 7):

• 5.1 Funzione di partizione fattorizzata: per la separazione delle scale temporali, Z = Z_trans·Z_vib·Z_rot·Z_bond si fattorizza.
• 5.2 Capacità termica per fase: c_v = d(U/N)/dτ; il termine positivo −(ε/2)(dz̄/dτ) genera il picco di tipo Schottky come analogo del calore latente.
• 5.3 Curva numerica c_v(τ): picco di c_v 5,51 a τ_m = 1,20, plateau solido 2,50, ramo gassoso → 1; il calore latente L ≈ 1,41 ε per molecola (con termine cooperativo) è triplicato rispetto a 0,48 ε (senza).
• 5.4 / 5.4a Due soglie separate: la percolazione di connettività z̄_c = z_max·p_c ≈ 2,08 (a τ_b ≈ 1,85) segna la disgregazione degli aggregati (vicino all'ebollizione) ed è un risultato; va tenuta nettamente distinta dalla rigidità generica di Maxwell z̄_c = 2d (densità di fusione, perdita di rigidità).
• 5.5 Ordine della transizione di fusione: criterio ln G > 2 + 4ε/K; il 2D (G = 6) è continuo per ogni J, il 3D (G = 12) è del primo ordine da K > ~8,2 (dimostrazione: Δz̄ ≈ 5,3, ΔU ≈ 1,75 ε).
• 5.6 Studio di dimensione finita: p*(L) converge monotonamente dal basso verso p_c = 0,3473; z̄_c ≈ 2,08 è robusto rispetto alla dimensione finita.
• 5.7 Scansione 3D completa di c_v: da K > K_crit ≈ 8,25 (J > 0,68) si apre il vero salto Δz̄ con calore latente ΔU; il picco di fusione si sposta con J verso τ maggiori (spostamento positivo di +0,5 … +2,2 rispetto al picco 2D continuo).

<a id="kap5-1"></a>

5.1 Funzione di partizione fattorizzata

Per la separazione delle scale temporali (SVELTO) la funzione di partizione canonica si fattorizza:

Z = Z_trans · Z_vib · Z_rot · Z_bond.

• Z_trans = (1/N!)·(A/λ_T²)^N, λ_T = h/√(2π m τ) (traslazione 2D classica ideale; MOLTITUDINE, CONSERVAZIONE, VOLO LIBERO).
• Z_vib = (τ/ħω₀)^(N_bond) (oscillatore armonico classico per legame intatto; TRIADE-b).
• Z_rot = ∏_(blocchi≥2) (2π I_b τ)^(1/2)/ħ (rotatore rigido 2D classico per blocco; TRIADE-c).
• Z_bond = Σ_{s} exp[(ε/τ)·Σ s_ij + (K/τ)·termine cooperativo] (somma di legami di tipo Ising; LEGAME, FRATTURA). La valutazione in campo molecolare fornisce x(τ), z̄(τ) da 4.2.

<a id="kap5-2"></a>

5.2 Capacità termica per fase (analitica)

Da U/N (cap. 4.3), c_v = d(U/N)/dτ (in unità k_B):

c_v = [1 + z̄/4 + ½ n_blk] + τ·[¼ (dz̄/dτ) + ½ (dn_blk/dτ)] − (ε/2)·(dz̄/dτ).

Il termine −(ε/2)·(dz̄/dτ) è positivo (dz̄/dτ < 0: riscaldare rompe legami) e domina là dove molti legami si rompono contemporaneamente — il picco di tipo Schottky, analogo del calore latente.

<a id="kap5-3"></a>

5.3 Curva numerica c_v(τ) con termine cooperativo — calore latente quantificato

Legame autoconsistente (minimo globale), n_blk(τ) da un Monte Carlo di percolazione di legami (union-find, L = 100, 5 realizzazioni), c_v = dU/dτ (differenze centrate), U/N incluso −(K/2)x² (cv_scan_09.py, ε = 1, J = 0,55, G = z_max = 6). Punti di appoggio scelti:

Risultato (numerico). L'ordine delle tre fasi resta confermato ed è ora più netto:

c_v(solido) ≈ 2,50 < c_v(liquido, massimo) ≈ 5,51 > c_v(ramo gassoso) ↓ → 1.

Calore latente (area del picco sopra la linea di base). Integrando (c_v − linea di base) sulla finestra di fusione τ ∈ [0,70 ; 1,70] (linea di base = retta per c_v(0,70) e c_v(1,70)) si ottiene

L ≈ 1,41 ε per molecola (con termine cooperativo), L ≈ 0,48 ε per molecola (senza termine cooperativo, vecchio bilancio), ΔL ≈ +0,92 ε per molecola.

Il termine cooperativo innalza il picco di un fattore ≈ 1,7 (3,25 → 5,51) e triplica il calore latente. Fisicamente: la rottura simultanea di molti legami vicini (FRATTURA) immagazzina non solo energia di configurazione (−ε/2·z̄) ma anche energia di campo molecolare (−K/2·x²) — entrambe si liberano fondendo.

<a id="kap5-4"></a>

5.4 Soglia di percolazione di connettività come risultato

La frazione dell'aggregato gigante P_∞(τ) collassa a τ ≈ 1,85; lì x ≈ 0,347 = p_c(triangolo) = 2·sin(π/18) = 0,3473, dunque z̄_c = z_max·p_c ≈ 2,08 — un risultato, non un dato imposto. Questa è la soglia di connettività (aggregato che attraversa) e segna la disgregazione degli aggregati vicino all'ebollizione (τ_b); non va confusa con la rigidità di Maxwell (z̄_c = 2d, fusione) — cfr. la separazione netta al cap. 5.4a.

<a id="kap5-4a"></a>

5.4.a Due soglie distinte nettamente separate

La doppia denominazione finora vaga «percolazione/rigidità» è ora risolta. Si tratta di due soglie fisicamente diverse che non vanno confuse:

Motivazione. La fusione è la perdita della rigidità meccanica di un blocco denso, non la disgregazione della connessione. La grandezza rilevante è dunque il conteggio dei vincoli di Maxwell (campo medio / rigidità generica): una rete D-dimensionale di forze centrali diventa isostatica esattamente quando la coordinazione media raggiunge z̄ = 2D (2D → 4, 3D → 6). La connettività dei legami (z̄_c = 2,08) sta nettamente al di sotto; una rete appena connessa a z̄ = 2,08 è da tempo molle (modi soffici), dunque già liquida. Perciò la connettività segna la disgregazione degli aggregati più tardiva, vicino all'ebollizione, mentre la rigidità fissa la solidificazione/liquefazione più precoce (densità di fusione). Questa separazione fissa il punto triplo del cap. 7.8 senza parametri e in modo realistico; non cambia la soglia di connettività 2,08 confermata in 5.4/5.6, ma solo la sua lettura (τ_b = vicino all'ebollizione, non la fusione).

<a id="kap5-5"></a>

5.5 Ordine della transizione di fusione

Criterio analitico. La regione spinodale (due minimi di g) esiste per τ < K/4; la transizione si trova approssimativamente a τ_t ≈ (ε + K/2)/ln G. Un vero salto (transizione del primo ordine) richiede τ_t < K/4, dunque nel limite K → ∞:

ln G > 2 + 4ε/K.

Ne segue un risultato notevole sulla degenerazione di rottura G = z_max:

Conferma numerica (max |Δx| per passo in τ, detect_first_order):

Lettura. Nel modello 2D canonico (G = 6) la fusione è per principio continua (una transizione «morbida» / crossover) — per nessun J nasce un vero salto, perché ln 6 < 2. Solo la maggiore coordinazione/degenerazione dell'impacchettamento 3D (G = 12) consente, a partire da una cooperatività sufficiente (K > 8,2), una transizione del primo ordine con calore latente finito. Ciò concorda con la realtà: la fusione 2D è spesso continua (di tipo KTHNY), la fusione 3D è del primo ordine. J = 0 (nessun campo molecolare) è sempre il caso di riferimento continuo. — Fisicamente J indica il rapporto fra energia cooperativa ed energia di legame (K/ε = J·z_max/ε): quanto un legame vicino intatto stabilizza il proprio.

Dimostrazione del primo ordine (motivata dal 3D, cv_scan_09.py, J = 1,0, G = 12, z_max = 12): z̄ salta di Δz̄ ≈ 5,30 a τ ≈ 2,82; il calore latente come salto di energia ΔU ≈ 1,75 ε per molecola. Lo stesso K su un reticolo G = 6 (2D) resta continuo — il carattere dell'ordine dipende unicamente da G = z_max.

<a id="kap5-6"></a>

5.6 Studio di dimensione finita della soglia di percolazione

Percolazione di legami su reticoli triangolari L × L (periodici), frazione di aggregato gigante P_∞(p) scalata su L (cv_scan_09.py; L = 50/100/200/400 con 8/5/3/2 realizzazioni). Soglia esatta p_c = 2·sin(π/18) = 0,34730 → z̄_c = z_max·p_c = 2,0838.

P_∞(p) (media ± dispersione):

Stima della soglia p*(L) (P_∞ attraversa 0,5):

Risultato. p*(L) converge monotonamente dal basso verso p_c = 0,3473; la dispersione a p_c cala da ±0,103 (L = 50) a ±0,017 (L = 400) — la transizione si acuisce con L ma non sposta la soglia (valore di percolazione universale). Con ciò z̄_c ≈ 2,08 = z_max·p_c è confermato come valore del modello robusto rispetto alla dimensione finita.

<a id="kap5-7"></a>

5.7 Scansione 3D completa di c_v(τ) della fusione del primo ordine

La singola serie dimostrativa del primo ordine del cap. 5.5 (J = 1,0, G = 12) è ora ampliata a una scansione completa di c_v(τ) sulla cooperatività J (firstorder_3d_16.py, connettività FCC per n_blk, bilancio energetico 3D U/N = τ·[3/2 + z̄/4 + (3/2)n_blk] − (ε/2)z̄ − (K/2)x²). Il criterio analitico d'ordine del cap. 5.5 (primo ordine ⇔ K > K_crit = 4ε/(ln 12 − 2) = 8,25) separa i casi; il picco di riferimento 2D continuo si trova a τ_2D = 1,20 (c_v ≈ 5,5).

*caso continuo: Δz̄/ΔU sono il crossover ripido (non un vero salto); τ_transizione è lì il picco di c_v.

Risultato. (i) Il caso 3D continuo J = 0,55 dà τ_transizione = τ_m ≈ 1,72 — la verifica incrociata con la cinetica di legame indipendente (bondkinetics_3d_15.py: τ_m(FCC) = 1,70) concorda. (ii) A partire da K > K_crit ≈ 8,25 (J > 0,68) si apre il vero salto Δz̄ con calore latente finito ΔU (salto di energia attraverso la discontinuità) — la fusione diventa del primo ordine, esattamente come richiede il criterio ln 12 > 2. (iii) La verifica incrociata J = 1,0 riproduce il valore dimostrativo Δz̄ ≈ 5,4 del cap. 5.5. (iv) L'altezza del salto Δz̄ e la posizione τ_transizione crescono con J: la transizione si sposta verso l'alto — la coordinazione più elevata dell'impacchettamento 3D e la cooperatività più forte innalzano la temperatura di fusione. Rispetto al picco 2D continuo (τ_2D = 1,20) si ottiene dunque ovunque uno spostamento positivo del picco (+0,5 … +2,2) che cresce monotonamente con K. Il carattere del primo ordine della fusione 3D è così documentato non più in un solo punto, ma su tutto l'intervallo di cooperatività.

<a id="kap6"></a>

6 — Plausibilità matematica dell'interpretazione 3D

Riassunto

Il capitolo 6 rende plausibile la generalizzazione dell'immagine dei legami da 2D a 3D in modo puramente matematico (senza una visualizzazione HTML propria):

• 6.1 Grado di libertà di rotazione 2D → 3D: la rotazione di blocco passa da un asse a tre, la traslazione da due direzioni a tre, z_max da 6 a 12; i valori base di c_v salgono (c_v,solido 2,5 → 4,5, c_v,gas 1 → 3/2), la forma della curva c_v(τ) resta.
• 6.2 La connettività come tratto distintivo: z_max è una grandezza quantizzata e separa di colpo l'impacchettamento più denso (FCC/HCP, z_max = 12) dal BCC (z_max = 8), senza stato intermedio.
• 6.3 Transizione aggiuntiva solido(1)→solido(2): dalla discretezza segue il polimorfismo FCC↔BCC con un gradino di c_v privo di parametri di 1,0 k_B; nell'immagine ideale rigorosa (θ_BCC = θ_FCC) esso non compare, e solo un BCC più molle (θ_BCC < θ_FCC, sotto la soglia di apertura 0,61) apre τ_ss < τ_m.
• 6.4 / 6.5 Considerazioni: l'anomalia acustica dell'acqua liquida come cambio sfumato di ordine locale fra due tipi di reticolo, e un'alternativa di campi di vortici frattali alla rotazione delle micelle (puramente qualitativa).
• 6.6 / 6.8 / 6.9 θ(reticolo) microscopico su tre livelli: modello a molle di Einstein (θ_BCC/θ_FCC ≈ 0,82…0,94), densità di stati fononici di tipo Debye (≈ 0,87…1,00) e reticolo di Lennard-Jones completo (≈ 1,0; BCC ideale instabile a taglio, FCC stato fondamentale stabile) — tutti sopra la soglia 0,61, sicché il polimorfismo reale sta fuori dall'immagine ideale armonica.
• 6.7 Metastabilità (2D): calcolo d'ordine di grandezza per la trasformazione triangolato ↔ quadrato (ΔU = +1 ε, ΔA ≈ +15,5 %) insieme alla barriera di nucleazione ΔG\* = πλ²/Δμ_A come contributo cinetico.

Questo capitolo rende plausibile la generalizzazione 3D in modo puramente matematico. Il 3D non viene mai visualizzato in una versione HTML. Lo scopo è unicamente mostrare come il grado di libertà di rotazione e la connettività cambino da 2D a 3D e quali nuove transizioni divengano così concepibili.

<a id="kap6-1"></a>

6.1 Variazione del grado di libertà di rotazione 2D → 3D

L'unico punto qualitativamente nuovo rispetto al 2D è la rotazione dei blocchi:

• 2D: un blocco rigido ha 1 GDL di rotazione (rotazione nel piano) → U_rot/N = ½ n_blk τ.
• 3D: un blocco rigido ha 3 GDL di rotazione (tre assi di Eulero) → U_rot/N = (3/2) n_blk τ.

Allo stesso tempo la traslazione passa da 2 a 3 GDL (U_trans/N = (3/2)τ) e il numero di coordinazione massimo da z_max = 6 (triangolo) a z_max = 12 (impacchettamento compatto di sfere). Il bilancio energetico in 3D recita dunque

U/N = τ·[3/2 + z̄/4 + (3/2) n_blk] − (ε/2)·z̄,

e i valori base della capacità termica si spostano coerentemente:

La struttura della curva (plateau solido, massimo liquido tramite −(ε/2)dz̄/dτ, caduta verso il valore del gas) resta identica; solo i livelli scalano con la dimensione e con z_max.

<a id="kap6-2"></a>

6.2 Influenza dei diversi reticoli cristallini — la connettività come tratto distintivo

In 3D vanno distinti almeno due reticoli ideali, che si separano per la loro connettività (numero di coordinazione z_max):

La connettività z_max è il tratto distintivo decisivo: fissa sia il livello di capacità termica del solido (z_max/4) sia la soglia di percolazione z̄_c. Poiché ogni reticolo ammette solo un numero di coordinazione fisso, la connettività è una grandezza quantizzata e finita — non esiste transizione continua fra z_max = 12 e z_max = 8, ma solo il salto discreto da un tipo di reticolo all'altro.

<a id="kap6-3"></a>

6.3 La transizione aggiuntiva solido(1) → solido(2)

Dalla finitezza/discretezza della connettività segue una transizione assente nel modello 2D: solido(1) → solido(2) (polimorfismo/allotropia). Due fasi solide con z_max diverso (per esempio FCC ↔ BCC) hanno entalpie libere μ(τ, p) diverse; alla combinazione temperatura/pressione μ_{solido(1)} = μ_{solido(2)} la sostanza cambia reticolo senza fondere. Nell'immagine del modello la geometria di coordinazione cambia di colpo (non solo l'occupazione z̄ come nella fusione) — l'unica situazione in cui la geometria di riferimento, altrimenti fissa, salta essa stessa. Sarebbe osservabile un piccolo gradino di c_v (variazione di z_max/4) più un calore latente di trasformazione, nettamente più debole del picco di fusione.

Il criterio del primo ordine ln G > 2 + 4ε/K (cap. 5.5) fornisce ora una ragione concreta per cui proprio l'interpretazione 3D è qualitativamente più ricca: solo lì (G = z_max = 12, ln 12 = 2,485 > 2) la fusione può diventare del primo ordine — e solo lì esiste il cambio discreto di reticolo solido(1) → solido(2) (FCC ↔ BCC), perché la connettività è quantizzata/finita. Entrambi i fenomeni sono la stessa causa in due vesti: la coordinazione più elevata e discreta dell'impacchettamento 3D genera sia il netto salto di fusione sia il polimorfismo. Il reticolo 2D isolato (G = 6) non ha né l'uno né l'altro — coerentemente con il cap. 5.5.

Quantificazione (bondkinetics_3d_15.py). La transizione viene ora calcolata in modo snello. Entrambi i reticoli ricevono un'entalpia libera nell'immagine armonica classica (le parti comuni p·v e di traslazione si elidono nel confronto):

μ(τ) = h₀ − c_v·τ·ln(τ/θ), h₀ = −(ε/2)·z_max, θ = ħω/k_B (rigidità fononica in unità ε),

con i valori caratteristici privi di parametri dei cap. 6.1/6.2:

Da qui il gradino di c_v senza parametri: |Δc_v| = variazione di z_max/4 = (12−8)/4 = 1,0 k_B (nella transizione FCC→BCC il solido diventa più «leggero» di 1 k_B in capacità termica) — molto minore del picco di fusione (diversi k_B, cap. 5.3/7.11).

La temperatura di trasformazione τ_ss (μ_FCC = μ_BCC) e il calore latente di trasformazione L = τ_ss·(S_BCC − S_FCC):

Risultato (constatazione onesta). Nell'immagine ideale rigorosa (stessa rigidità fononica per tutti i legami, θ_BCC = θ_FCC) non esiste transizione solido(1)→solido(2): l'FCC vince sia l'entalpia (più bassa, −6 contro −4 ε) sia — per il maggior numero di legami (c_v = 4,5 > 3,5) — l'entropia vibrazionale; μ_FCC resta a ogni temperatura sotto μ_BCC. Solo una rigidità fononica dipendente dal reticolo θ_BCC < θ_FCC (il BCC più aperto vibra più mollemente e ha dunque più entropia) porta μ_BCC sotto μ_FCC ad alta temperatura — allora la trasformazione si apre, e si colloca sotto il picco di fusione τ_m, sicché nasce l'ordine solido(1)→solido(2)→liquido→gas (τ_ss < τ_m < τ_c ≤ τ_b(3D), cap. 7.11). La rigidità molle del BCC è un'estensione moderata fuori dall'immagine ideale rigorosa (nota nella realtà come polimorfismo γ/δ del ferro: BCC ad alta temperatura, FCC in mezzo); il gradino di c_v di 1,0 k_B e il piccolo calore di trasformazione positivo sono invece interni al modello e privi di parametri.

Quanto è davvero piccolo θ_BCC/θ_FCC? (theta_springs_16.py). Per liberare la trasformazione dal parametro libero θ, la rigidità fononica viene ora stimata microscopicamente (cap. 6.6): un modello a molle di Einstein con molle di coppia centrali dà θ_BCC/θ_FCC ≈ 0,82 (solo vicini più prossimi) fino a ≈ 0,94 (con la seconda shell del BCC). La soglia di apertura richiesta dal confronto dei μ sta invece a θ_BCC/θ_FCC < 0,61 (al di sotto esiste τ_ss e sta sotto τ_m). Il rapporto geometrico di rigidità è dunque troppo grande: nell'immagine ideale rigorosa l'FCC resta stabile a ogni temperatura, non vi è transizione solido(1)→solido(2) di origine puramente geometrica. Il polimorfismo reale (ferro γ↔δ) richiede contributi aggiuntivi (elettronici/magnetici/anarmonici) fuori dall'immagine armonica ideale — proprio la «estensione moderata» posta sopra. Ciò conferma la constatazione onesta e trasforma il rapporto θ da ipotesi libera in un limite del modello stimato e nominato.

<a id="kap6-4"></a>

6.4 Considerazioni: l'anomalia acustica dell'acqua liquida

I liquidi normali mostrano v(T) decrescente in modo monotono (cap. 3.3). L'acqua ha un massimo acustico verso 74 °C: v sale dapprima con T prima di calare. Nel modello ciò si può leggere come una concorrenza fra due tipi di reticolo locali (cap. 6.2): un ordine locale aperto, a bassa coordinazione, simile al ghiaccio (z_max piccolo, molti «buchi»), e un ordine locale denso, a coordinazione maggiore. A T bassa domina la struttura aperta; il riscaldamento rompe dapprima l'ordine locale aperto a favore di quello più denso — la connettività efficace e con essa K_eff salgono (v cresce). Solo sopra il massimo prevale l'ordinaria rottura di legami della FRATTURA, K_eff cala, v scende. L'anomalia è così espressione di un cambio di ordine locale di tipo solido(1)↔solido(2), sfumato e incompleto, all'interno della fase liquida — concepibile solo nel modello 3D con due tipi di reticolo, non nel reticolo 2D isolato.

<a id="kap6-5"></a>

6.5 Alternativa: campi di vortici frattali invece della rotazione delle micelle

Il moto browniano locale (cap. 3.4) è stato spiegato mediante la rotazione e lo scambio delle micelle. In alternativa lo si può immaginare come effetto di campi di vortici frattali: i confini di blocco, continuamente rimodellati dalla FRATTURA, generano un campo di correnti locali di quantità di moto che — analogamente ai flussi turbolenti — presenta vortici autosimili su tutte le scale di lunghezza. Una particella di prova è allora spinta non da urti micellari discreti, ma da un continuo di vortici invariante di scala. Entrambe le immagini (micelle discrete ↔ vortici frattali) sono casi limite complementari della stessa cinetica di FRATTURA; quale prevalga dipende dal rapporto fra la dimensione della particella di prova e la dimensione media delle micelle. Questa breve discussione è puramente qualitativa e non viene qui sviluppata matematicamente.

<a id="kap6-6"></a>

6.6 Rigidità fononica microscopica θ(reticolo): un modello a molle di Einstein

Il confronto solido(1)→solido(2) del cap. 6.3 dipende dalla rigidità fononica relativa θ_BCC/θ_FCC. Invece di sceglierla liberamente, la si stima qui a partire da coordinazione e geometria (theta_springs_16.py).

Modello. Un atomo centrale sta nella buca armonica dei suoi z vicini; ogni legame è una molla di coppia radiale κ. Uno spostamento u si proietta isotropicamente con ⟨cos²⟩ = 1/d = 1/3 (d = 3) sulle direzioni di legame, sicché la costante elastica efficace è k_eff = (z/3)·κ. La frequenza di Einstein è ω = √(k_eff/m), dunque θ = ħω/k_B ∝ √k_eff. A parità di massa e di molla di coppia segue

θ_BCC/θ_FCC = √(z_BCC/z_FCC) = √(8/12) = 0,816.

Raffinamenti. (B) La seconda shell del BCC (6 atomi a r₂/r₁ = 2/√3 = 1,155) con una molla scalata con la distanza κ(r) ∝ (r₁/r₂)ⁿ (n = 6…8, di tipo Lennard-Jones) porta il rapporto a 0,91…0,94 (il BCC risulta effettivamente più rigido per la vicina seconda shell). (C) Il legame più corto del BCC a parità di densità (r_nn(BCC)/r_nn(FCC) = 0,972, −2,8 %) irrigidisce leggermente la molla e dà 0,89…0,91. Complessivamente una finestra geometrica θ_BCC/θ_FCC ≈ 0,82…0,94.

Soglia di apertura. Il confronto dei μ (cap. 6.3) apre la trasformazione (τ_ss esiste e sta sotto τ_m = 1,70) solo per θ_BCC/θ_FCC < 0,61. La finestra geometrica sta nettamente sopra.

Constatazione. Il puro modello a molle non rende il BCC abbastanza molle: nell'immagine ideale rigorosa l'FCC resta stabile e una transizione solido(1)→solido(2) di origine geometrica non nasce. Con ciò θ_BCC/θ_FCC è passato da ipotesi libera a grandezza stimata, il cui valore naturale (≈ 0,8–0,9) tiene chiusa la trasformazione — coerentemente con la constatazione onesta del cap. 6.3. Il polimorfismo reale sta dunque dimostrabilmente fuori dall'immagine armonica ideale (contributi elettronici/magnetici, anarmonicità), e non nella mera geometria del reticolo.

<a id="kap6-7"></a>

6.7 Stima matematica della metastabilità (2D triangolato ↔ quadrato)

Per la discussione qualitativa del cap. 3.2, ecco il corrispondente calcolo d'ordine di grandezza nell'immagine 2D (dischi duri, diametro D, distanza di contatto a = D).

(1) Energia. Per particella U/N = −(ε/2)·z (ogni legame per metà): triangolato z = 6 → −3 ε, quadrato z = 4 → −2 ε. La trasformazione triangolato → quadrato costa dunque ΔU = +1 ε per particella (endotermica; la forma lasca è più ricca di energia).

(2) Volume (area) e lavoro. Area per particella: triangolato A_tri = (√3/2)·D² ≈ 0,866 D², quadrato A_quad = D². L'allargamento è ΔA = A_quad − A_tri ≈ 0,134 D² (+15,5 %) — l'aumento di volume mostrato nel fumetto. Contro una pressione (ridotta) p va fornito il lavoro W = p·ΔA; l'entalpia libera cresce di Δμ = ΔU + p·ΔA − τ·ΔS. Il termine entropico ΔS > 0 (la forma lasca ha modi più molli, un gradino di c_v maggiore Δc_v = (6−4)/4 = 0,5 k_B) abbassa μ_quad sotto μ_tri solo a τ alta — perciò il triangolato è lo stato fondamentale di bassa temperatura e il quadrato la forma di alta temperatura / metastabile.

(3) Barriera di nucleazione (2D). La metastabilità presuppone una barriera. In 2D un nucleo della nuova fase (raggio r) cresce con un guadagno di area −Δμ_A·π r² (Δμ_A = forza motrice per area) contro il contributo di bordo +λ·2π r (λ = energia di spigolo per lunghezza): ΔG(r) = −Δμ_A·π r² + 2π λ r. Massimo in r\* = λ / Δμ_A con la barriera ΔG\* = π λ² / Δμ_A. Il nucleo critico r\ si restringe al crescere della forza motrice Δμ_A; finché la scala termica τ = k_BT non raggiunge la barriera (**ΔG\ ≳ alcuni k_BT), la fase lasca (o, nell'ebollizione, quella soprassatura) resta congelata cineticamente. Per piccoli aggregati (2–3 particelle, come nella striscia sulla soprassaturazione) il rapporto bordo/superficie è massimamente sfavorevole — essi stanno a sinistra di r\* e in media si ridissolvono. Con ciò l'energia di attivazione mostrata nel fumetto è la stessa grandezza ΔG\* che sostiene sia la forma cristallina metastabile sia la fase gassosa soprassatura; è un contributo cinetico** al modello di equilibrio ideale (privo di barriera).

<a id="kap6-8"></a>

6.8 Raffinamento di tipo Debye di θ(reticolo)

Il modello a molle di Einstein (cap. 6.6) sostituisce l'intera distribuzione dei fononi con una sola frequenza. L'approccio presente raffina θ(reticolo) mediante una densità di stati fononici snella per reticolo (theta_debye_17.py): le stesse molle di coppia centrali κ(r) = r^(−n) sulla prima e sulla seconda shell di vicini formano la matrice dinamica D(k) = Σ_j κ_j·(R̂_j⊗R̂_j)·(1 − cos(k·R_j)), diagonalizzata su una griglia k della zona di Brillouin primitiva (Monkhorst–Pack, 12³ punti) mediante un risolutore analitico 3×3, sola libreria standard. Dalle tre branche fononiche seguono due misure di rigidità di tipo Debye come momenti di frequenza:

• 2º momento θ(2) ∝ √⟨ω²⟩ — il limite di Einstein (la media sulla ZB di 1 − cos dà 1, dunque ⟨ω²⟩ = z·κ/3); controllo, soli vicini più prossimi, pura geometria: √⟨ω²⟩_BCC / √⟨ω²⟩_FCC = 0,816 = √(8/12) — esattamente il valore di Einstein del cap. 6.6.
• momento logaritmico θ(0) ∝ exp⟨ln ω⟩ — la grandezza propriamente di tipo Debye, che pesa le branche acustiche molli che il valore di Einstein ignora.

Tutte le costanti dotate di scala (ħ/k_B, la scala di κ, il prefattore di Debye) si elidono nel rapporto θ_BCC/θ_FCC; resta solo la forma della DOS. Con la seconda shell di vicini (che è ciò che stabilizza meccanicamente il BCC) e una legge di molla n = 0…8 si ottiene

θ_BCC/θ_FCC ≈ 0,87…1,00 (banda d'errore su definizione del momento e legge di molla).

Ciò è piuttosto più alto del puro valore di Einstein (0,82…0,94), perché il legame più corto del BCC a parità di densità irrigidisce le molle. La soglia di apertura 0,61 richiesta dal confronto dei μ (cap. 6.3) resta nettamente al di sotto: neanche la DOS fononica rende il BCC abbastanza molle. La constatazione del cap. 6.6 è così confermata e precisata dall'approssimazione più fine — nell'immagine ideale rigorosa l'FCC resta stabile a ogni temperatura; non nasce alcuna transizione solido(1)→solido(2) di origine puramente geometrica, e il polimorfismo reale (ferro, cap. 7.13) necessita di contributi fuori dall'immagine armonica ideale.

<a id="kap6-9"></a>

6.9 θ(reticolo) da un reticolo di Lennard-Jones completo

Il cap. 6.8 usa una molla centrale troncata κ(r) = r^(−n) su nn + 2nn. La versione 17 verifica la stessa domanda con una vera costante di forza da un potenziale di coppia sull'intero reticolo (theta_ljlattice_18.py): ogni legame diventa una molla ricavata dal potenziale di Lennard-Jones φ(r) = 4ε[(σ/r)¹² − (σ/r)⁶], con la parte longitudinale κ_L = φ″(r) e la parte trasversale κ_T = φ′(r)/r. La matrice dinamica

D_ab(k) = Σ_j [ (κ_L − κ_T)·û_a û_b + κ_T·δ_ab ]·(1 − cos(k·R_j))

somma su tutte le shell di vicini fino a un raggio di taglio (2,0…3,0 nn); ogni reticolo è previamente rilassato al proprio equilibrio LJ (costante reticolare libera, densità uguale non imposta). È il test onesto per sapere se l'approssimazione snella nn+2nn (cap. 6.8) colga correttamente l'ordine delle rigidità.

Risultato (deterministico, sola libreria standard):

Due constatazioni, entrambe che rafforzano i cap. 6.6/6.8:

1. L'FCC è lo stato fondamentale LJ stabile (E ≈ −8,4 ε per atomo con taglio 3,0 nn, contro −8,0 ε per il BCC) — il noto primato dell'impacchettamento denso in un potenziale di coppia.
2. Il BCC LJ ideale è instabile a taglio: circa il 3…3,5 % dei modi fononici sono immaginari (ω² < 0, instabilità di Born). L'immagine ideale del puro potenziale di coppia non porta dunque alcun BCC meccanicamente stabile; sui modi reali rimanenti θ_BCC/θ_FCC ≈ 1,0 — la soglia di apertura 0,61 richiesta dal confronto dei μ a maggior ragione non viene superata verso il basso.

Con ciò l'affermazione dei cap. 6.6/6.8 è robusta rispetto al raffinamento dalla molla troncata al reticolo completo e dal modello puramente longitudinale alla matrice dinamica completa: nell'immagine ideale classica non nasce alcuna transizione solido(1)→solido(2) puramente geometrica — il polimorfismo reale (ferro cap. 7.13, titanio cap. 7.15) necessita di contributi fuori dal potenziale di coppia (elettronici/magnetici, anarmonici). Un vero calcolo fononico ab initio è deliberatamente lasciato fuori (abbandonerebbe il carattere snello dell'immagine ideale).

<a id="kap7"></a>

7 — Calibrazione ed equazione di stato

Riassunto. Il capitolo 7 calibra il modello su sostanze reali mediante l'unico parametro di energia ε e sviluppa il ramo liquido/gas in un'equazione di stato:

• 7.1 / 7.2 Metodo e tabella di calibrazione: accoppiamento lineare di temperatura T = (ε/k_B)·τ con un punto di ancoraggio al punto di fusione; da ciò il rapporto universale raffinato R = τ_b/τ_m ≈ 1,54.
• 7.3 / 7.4 Equazione di stato di van der Waals con nocciolo duro: l'EOS ridotta p_r = 8τ_r/(3v_r − 1) − 3/v_r² con punto critico (1, 1, 1) e coesistenza di Maxwell; a ∝ ε, b ∝ D². I punti di ebollizione stanno sistematicamente a T_b/T_c ≈ 0,48 (distanza di Guldberg); un calcolo di confronto Redlich–Kwong li innalza a ≈ 0,60 — come budget d'errore nominato, senza cambio di EOS.
• 7.5 Superficie (τ, v) comune: la cupola di coesistenza (ebollizione), la linea di fusione e la parete dura si compongono con un punto triplo (τ_r,triple ≈ 0,47 per v_solid = 0,40).
• 7.6 Calore latente: L ≈ 1,41 ε dà ΔS_fus ≈ 1,175 R (vicino alla regola di Richards); i metalli sono colti eccellentemente (rame ×1,00), mentre i cristalli ionici, a legame idrogeno e molecolari stanno fuori dall'immagine ideale.
• 7.7 Unificazione ε ↔ ε_lg: k_B T_c = z_max·ε/4 accoppia le due scale (ε_lg = (27/32)·z_max·ε); la previsione 3D τ_c/τ_m = 2,50 coglie le sostanze molecolari semplici, mentre il 2D resta caso particolare.
• 7.8 / 7.9 Punto triplo senza parametri: v_solid segue dalla soglia di rigidità z̄_c = 2d, non dalla connettività; in 3D si apre una finestra liquida netta (τ_m = 1,20 < τ_triple = 2,83 < τ_c = 3,00), in 2D resta chiusa (ruolo particolare del 2D, ormai triplicemente documentato).
• 7.10 / 7.11 Vere scansioni 3D: la percolazione di legami su FCC/BCC dà z̄_c(3D) ≈ 1,5–1,6; la scansione in τ della cinetica di legame mostra τ_m ≈ τ_rig (fusione = rigidità, basso) e τ_b(3D) ≈ τ_c (ebollizione = disgregazione della connettività al punto critico).
• 7.12 / 7.13 / 7.15 Linea solido(1)→solido(2): dal punto τ_ss la pressione fa una linea (dp/dτ = ΔS/Δv > 0, Δv = +8,9 %), inserita come quarta linea di confine nella superficie (τ, v) e calibrata come linea (τ, p) sul ferro (γ→δ), confermata sul titanio (α→β) e verificata con l'argon come controllo negativo.
• 7.14 Ramo Redlich–Kwong: come variante autonoma dell'immagine ideale con calibrazione propria accanto all'impostazione predefinita di van der Waals (T_b/T_c ≈ 0,60, errore −1 %).

<a id="kap7-1"></a>

7.1 Metodo di calibrazione (un parametro di energia ε) — adattato

Lineare: T = (ε/k_B)·τ, un punto di ancoraggio per sostanza. Ancoraggio al punto di fusione: ε/k_B = T_m,exp/τ_m, previsione T_b,pred = T_m,exp·(τ_b/τ_m). Con il picco di fusione raffinato (termine cooperativo) il massimo di c_v sta ora a τ_m = 1,20 (invece di 1,15), mentre la percolazione di connettività (vicino all'ebollizione) resta invariata a τ_b = 1,85. Da ciò il rapporto universale raffinato:

R = τ_b / τ_m = 1,85 / 1,20 ≈ 1,54 (versione 08: 1,61).

La tabella di calibrazione della versione 08 (cap. 7.2) scala linearmente con R: tutti i T_b,pred calano del fattore 1,54/1,61 ≈ 0,957 (≈ −4,3 %). Ciò migliora leggermente la sovrastima dei gas nobili (argon +54 % → +48 %) e peggiora un poco i risultati per NaCl e ferro (NaCl −0,6 % → −4,8 %). L'enunciato di base resta: R vale per sostanze a impacchettamento denso isotropo e legame non direzionale; gas nobili e catene stanno fuori dall'immagine ideale.

<a id="kap7-2"></a>

7.2 Tabella di calibrazione

Ripresa invariata da modell_08.md cap. 7.2 (con la nota sullo scaling lineare in R del 7.1). Riproducibile via cv_scan_08.py (tabella) e cv_scan_09.py (τ_m, τ_b, R).

<a id="kap7-3"></a>

7.3 Asse densità/volume: equazione di stato di van der Waals con nocciolo duro

La pressione passa da «parametro» a vera variabile. Il NOCCIOLO DURO fornisce il volume proprio b, l'attrazione di legame (ε·z̄) la pressione di coesione a. In unità ridotte (punto critico in (v_r, τ_r, p_r) = (1, 1, 1)):

p_r = 8 τ_r / (3 v_r − 1) − 3 / v_r², v_r > 1/3 (parete dura).

Coesistenza per costruzione di Maxwell (cv_scan_09.py, taglio ad aree uguali):

Risultato. Sotto τ_c si apre il salto di densità l/g (v_gas/v_liq cresce da 2,5 a 11,5, tanto maggiore quanto più freddo). Questa è la linea di ebollizione — ora nettamente separata dal picco di c_v della fusione (cap. 5.3), che vive su un reticolo di densità fissa. Con ciò si chiude il punto 3 rimasto aperto nella versione 08: fusione (ordine di rigidità/legame a densità fissa) ed ebollizione (salto di densità/volume) sono assi diversi dello stesso modello.

<a id="kap7-4"></a>

7.4 Calibrazione quantitativa l/g: a → ε, b → D²

Applicazione delle grandezze ridotte a sostanze reali. Le costanti di van der Waals dipendono dal punto critico:

b = k_B·T_c / (8 p_c) (volume proprio, NOCCIOLO DURO), a = 27 k_B²·T_c² / (64 p_c) (coesione, LEGAME).

Così a ∝ ε (la coesione proviene dall'energia di legame) e b ∝ D² (il volume proprio 2D dal diametro del nocciolo duro). La scala di energia naturale del ramo l/g è

ε_lg / k_B = (a/b) / k_B = (27/8)·T_c = 3,375·T_c.

Sostanza di ancoraggio argon (T_c = 150,87 K, p_c = 48,63 bar): b = 5,35·10⁻²⁹ m³ per molecola, a = 3,76·10⁻⁴⁹ J·m³ per molecola², ε_lg/k_B = 509 K e un diametro molecolare D ≈ b^{1/3} = 3,8·10⁻¹⁰ m — accurato all'11 % rispetto al diametro di van der Waals dell'argon (~3,4·10⁻¹⁰ m).

Previsione della temperatura di ebollizione. A pressione esterna fissa, il punto di ebollizione è là dove è soddisfatta la pressione di saturazione ridotta p_sat(τ_r)·p_c = p_atm; da ciò T_b = T_c·τ_r,b (curva di saturazione di Maxwell da lg_calib_10.py):

*Il CO₂ sublima a 1 atm (T = 194,7 K = punto di sublimazione).

Risultato. Per le sostanze molecolari (p_c ≈ 30–70 bar) il modello dà ovunque T_b/T_c ≈ 0,48–0,51, sistematicamente ~17–24 % sotto il valore misurato. Ciò non è un errore del modello in senso stretto, ma la ben nota distanza dell'equazione di van der Waals dalla regola empirica di Guldberg (T_b/T_c ≈ 0,6): l'isoterma di van der Waals sottostima la pressione di saturazione alle basse temperature. L'acqua (e i metalli, i cui T_c/p_c sono comunque solo stimati grossolanamente) stanno con p_r,atm < 0,005 fuori dall'intervallo sensato di van der Waals (τ_r,b < 0,40) — lì la coesione non è più di tipo van der Waals (legami a idrogeno, legame metallico). L'ordine di grandezza e la graduatoria dei punti di ebollizione sono colti correttamente.

Inquadramento di Guldberg come tabella d'errore. La distanza è qui solo nominata, non corretta con un cambio di EOS (l'immagine ideale resta di tipo van der Waals). La regola empirica di Guldberg (Guldberg 1890) afferma che il punto di ebollizione si colloca verso T_b/T_c ≈ 0,6 (con lieve dipendenza dalla sostanza, ~0,58–0,64). Il modello ridotto di van der Waals dà invece un valore costante T_b/T_c = τ_r,b ≈ 0,48–0,51:

Lettura. La media del modello, 0,49, coglie la costante di Guldberg 0,60 a circa −18 %; lo scarto è sistematico e unilaterale (sempre troppo basso), dunque una proprietà dell'isoterma di van der Waals (pressione di saturazione troppo piccola a τ bassi), non dell'immagine molecolare. Una migliore EOS di coesione (per esempio Redlich–Kwong) innalzerebbe τ_r,b — ma resta qui deliberatamente fuori dall'immagine ideale.

Budget d'errore con un calcolo di confronto Redlich–Kwong (guldberg_rk_15.py). Per quantificare la distanza di van der Waals dalla regola di Guldberg (non per sanarla), lo stesso test del punto di ebollizione viene ripetuto con una coesione di tipo Redlich–Kwong, dipendente dalla temperatura, a/√T. In unità ridotte (punto critico (v_r, τ_r, p_r) = (1, 1, 1), Z_c = 1/3):

p_r = 3 τ_r/(v_r − β) − A/(√τ_r · v_r(v_r + β)), β = 3·Ω_b = 0,25992, A = 9·Ω_a = 3,84732,

con costruzione di Maxwell ad aree uguali e la stessa pressione esterna ridotta p_r,atm di sopra. Risultato (un valore per sostanza di ancoraggio, T_b/T_c = τ_r,b):

Risultato. La coesione RK dipendente dalla temperatura innalza T_b/T_c da vdW ≈ 0,48 a ≈ 0,60 — praticamente la costante di Guldberg e i valori reali (≈ 0,58–0,63); l'errore medio rispetto alle misure cala da −21 % (vdW) a −1 % (RK). Con ciò la causa della distanza di Guldberg è chiaramente quantificata: è proprio la mancanza di dipendenza dalla temperatura della coesione di van der Waals (pressione di saturazione troppo piccola a τ bassi), non l'immagine molecolare. Questo calcolo di confronto sta deliberatamente fuori dall'immagine ideale; il modello resta di tipo van der Waals (nessun cambio di EOS) — il budget d'errore è solo nominato, non incorporato.

<a id="kap7-5"></a>

7.5 Superficie (τ, v) comune con punto triplo

Gli assi di fusione ed ebollizione hanno finora vissuto separati: la fusione (ordine di rigidità/legame) a densità fissa, l'ebollizione (van der Waals) attraverso il salto di densità. In un piano comune (τ_r, v_r) essi si compongono:

• Cupola di coesistenza (ebollizione): v_liq(τ_r) (ramo sinistro) e v_gas(τ_r) (ramo destro) dalla costruzione di Maxwell, che si incontrano al punto critico (τ_r, v_r) = (1, 1).
• Solido (fusione): banda stretta appena sopra la parete dura v_r > 1/3, a densità elevata e quasi fissa (volume ridotto v_solid ≈ 0,40).
• Linea di fusione: quasi verticale a v_solid (densità fissa), alla temperatura della rigidità di Maxwell (z̄_c = 2d; il cap. 7.8 ne ricava v_solid senza parametri — la vecchia attribuzione della versione 11, «z̄_c ≈ 2,08», è così corretta).

Il punto triplo è il punto in cui il confine di coesistenza liquida v_liq(τ_r) raggiunge la densità d'impacchettamento solido v_solid — al di sotto non esiste più ramo liquido (sublimazione diretta):

Risultato. Per v_solid = 0,40 il punto triplo sta a τ_r,triple ≈ 0,47, p_sat ≈ 0,018, con un volume di vapore v_gas ≈ 68 (vapore di sublimazione molto diluito). Il rapporto τ_triple/τ_c ≈ 0,47 sta nella finestra reale (acqua 0,42; argon 0,56; CO₂ 0,71). Più denso è il solido (v_solid minore), tanto più basso il punto triplo ridotto e tanto più diluito il vapore coesistente — qualitativamente corretto. Con ciò entrambe le transizioni e il punto triplo sono ancorati in una sola superficie (τ, v); il punto critico siede in (1, 1, 1). La versione 15 aggiunge al cap. 7.12 una quarta linea di confine — la trasformazione solido(1)/solido(2) (FCC↔BCC) come due bande di densità solida quasi verticali (v ≈ 0,667 e 0,726) che si convertono l'una nell'altra a τ_ss.

Versione anteriore nello sviluppo: qui v_solid = 0,40 era ancora imposto. Il cap. 7.8 ricava v_solid senza parametri dalla soglia di percolazione z̄_c e fa del punto triplo una vera previsione.

<a id="kap7-6"></a>

7.6 Calore latente a fronte delle entalpie di fusione

Il calore latente del modello L ≈ 1,41 ε (cap. 5.3) è tradotto in unità reali mediante la calibrazione al punto di fusione ε/k_B = T_m/τ_m (τ_m = 1,20):

ΔH_m,pred = L·R·T_m/τ_m = 1,175·R·T_m, ossia ΔS_fus,pred = L·R/τ_m = 1,175·R.

Questa è — cosa notevole — quasi esattamente la regola di Richards (ΔS_fus ≈ R per solidi di costruzione semplice). Confronto con entalpie di fusione misurate (lg_calib_10.py):

Risultato. Il modello coglie ovunque l'ordine di grandezza (kJ/mol) e in modo eccellente la regola di Richards per i metalli semplici (rame ×1,00, ferro ×1,28); i gas liquefatti stanno appena sotto (metano ×0,94, azoto ×0,87, argon ×0,69). I cristalli ionici, a legame idrogeno e molecolari (NaCl ×0,37, acqua ×0,44, benzene ×0,28, CO₂ ×0,23) hanno entropie di fusione reali nettamente maggiori (ΔS_exp/R = 2,6–5,0), perché fondendo si liberano gradi d'ordine aggiuntivi (appaiamento ionico, orientazione di molecole poliatomiche) che il modello a un solo parametro d'ordine non coglie. Il limite è così chiaramente nominato: l'immagine ideale vale per sostanze a legame isotropo e non direzionale (senza fit).

<a id="kap7-7"></a>

7.7 Unificazione delle scale di energia: ε ↔ ε_lg

Finora il ramo di fusione portava l'energia di legame ε (τ_m = 1,20) e il ramo l/g una scala di coesione separata ε_lg/k_B = 3,375·T_c (cap. 7.4). Entrambe si possono ricavare da un solo parametro del modello. Nel modello la coesione di van der Waals a nasce dalla somma delle energie di legame per molecola (LEGAME): ogni molecola porta z_max/2 legami di energia ε. Un gas reticolare di Bragg–Williams in campo medio con attrazione fra vicini ε fornisce la relazione classica per la temperatura critica:

k_B·T_c = z_max·ε / 4, dunque τ_c ≡ k_B T_c/ε = z_max/4.

Inserendo ciò in ε_lg = a/b = (27/8)·k_B·T_c (cap. 7.4) segue l'accoppiamento delle due scale:

ε_lg = (27/8)·(z_max/4)·ε = (27/32)·z_max·ε, costante c = 27/32 = 0,844.

Con ciò picco di fusione, percolazione e punto critico stanno su un'unica scala ε (unify_scale_12.py):

Previsione τ_c/τ_m = z_max/(4·τ_m): 2D → 1,25, 3D → 2,50. A fronte di sostanze reali:

Risultato. La previsione 3D τ_c/τ_m = 2,50 coglie bene la cifra caratteristica reale delle sostanze molecolari semplici (T_c/T_m ≈ 2,0–2,1); i metalli stanno più in alto (T_c stimato grossolanamente, altro tipo di legame). La previsione 2D (1,25) è troppo bassa — e τ_c(2D) = 1,50 sta addirittura sotto τ_b = 1,85 (il punto critico starebbe sotto la soglia di percolazione). Non è un errore di calcolo, ma il noto ruolo particolare del 2D: lì le scale di coesione e di legame non si ordinano nettamente, manca una finestra liquida robusta (cfr. la fusione 2D continua, cap. 5.5). Solo in 3D (τ_m = 1,20 < τ_c = 3,00) le scale si separano in modo fisicamente sensato.

<a id="kap7-8"></a>

7.8 Punto triplo senza parametri a partire da z̄_c

Il cap. 7.5 imponeva ancora v_solid = 0,40. La versione anteriore nello sviluppo ricava la densità d'impacchettamento solido dalla soglia di percolazione. Il numero di coordinazione di un reticolo a volume ridotto v_r segue la relazione densità-coordinazione (parete dura v_r = 1/3 ↔ impacchettamento più denso, z̄ = z_max):

z̄(v_r) = z_max·(1/3)/v_r = z_max/(3 v_r), inversamente v_solid = z_max/(3·z̄_c).

Il solido perde la propria rigidità non appena z̄ scende sotto la soglia z̄_c; la densità corrispondente è v_solid. Il punto triplo è allora il punto in cui il confine di coesistenza liquida v_liq(τ_r) raggiunge questa densità (unify_scale_12.py, 2D, z_max = 6):

Versione migliorata. La densità di fusione v_solid è fissata dalla rigidità generica di Maxwell z̄_c = 2d (2D → 4, 3D → 6), non dalla connettività dei legami (motivazione al cap. 5.4a: la fusione è perdita di rigidità, non di connessione). Con ciò il punto triplo affidabile recita (rigidity_3d_13.py, 2D):

Risultato 2D. Anche con la soglia di rigidità (fisicamente corretta) il reticolo 2D isolato non dà alcuna finestra liquida netta: v_solid = 0,50 fornisce sì il valore reale τ_triple/τ_c ≈ 0,77, ma in unità ε assolute τ_triple = 0,77·1,50 = 1,16 e sta dunque sotto τ_m = 1,20. La soglia di connettività (z̄_c = 2,08) spinge addirittura v_solid ≈ 0,96 ≈ v_c fuori dalla cupola di coesistenza — punto triplo e punto critico coincidono. Entrambi i calcoli confermano il ruolo particolare del 2D (cap. 5.5, 7.7): il reticolo 2D isolato regge a stento una fase liquida autonoma. È una constatazione onesta, non una debolezza di calcolo o di modello in senso stretto.

Il test decisivo è dunque il calcolo di controllo 3D (cap. 7.9): lì la stessa regola di rigidità (z̄_c = 6 → v_solid = 0,667) libera una finestra netta. In ogni caso il punto triplo è ora legato senza parametri a z̄_c — il parametro libero v_solid è eliminato.

<a id="kap7-9"></a>

7.9 Calcolo di controllo 3D: una finestra liquida netta su un'unica scala ε

Verifica se il collasso 2D si risolva in 3D quando tutte le marche stanno su un'unica scala ε (τ_c = z_max/4 dal cap. 7.7; v_solid dalla rigidità generica z̄_c = 2d; punto triplo là dove la densità liquida ridotta di van der Waals v_liq(τ_r) raggiunge la densità del solido). Risultato (rigidity_3d_13.py):

Risultato. In 3D si apre una finestra liquida netta: τ_m = 1,20 < τ_triple = 2,83 < τ_c = 3,00 (larghezze: τ_triple − τ_m = 1,63; τ_c − τ_triple = 0,18). Il bordo superiore è stretto — il punto triplo sta appena sotto quello critico — ma l'ordine è giusto e la fase liquida esiste come campo autonomo. In 2D la finestra resta chiusa (τ_triple = 1,16 < τ_m). Con ciò il ruolo particolare del 2D è documentato per la terza volta in modo indipendente ed è insieme mostrato che si tratta di un effetto dimensionale: la stessa regola di rigidità sana il modello in 3D. Il bordo superiore stretto è corroborato al cap. 7.10 da una vera scansione di percolazione 3D e da un'analisi di sensibilità.

Riquadro dei limiti del modello — il ruolo particolare del 2D. Tre calcoli indipendenti mostrano la stessa constatazione: il reticolo 2D isolato regge a stento una fase liquida autonoma. 1. τ_c < τ_b (cap. 7.7): il punto critico (1,50) starebbe sotto la percolazione di connettività (1,85) — le scale di coesione e di legame non si ordinano nettamente in 2D. 2. v_solid ≈ v_c (cap. 7.8): la densità di fusione di connettività (0,96) coincide quasi con la densità critica — punto triplo e punto critico collassano. 3. τ_triple < τ_m (cap. 7.9): perfino con la corretta rigidità di Maxwell il punto triplo (1,16) sta sotto il picco di fusione (1,20) — nessuna finestra.

È una constatazione onesta, non una debolezza di calcolo o di modello: è il noto effetto dimensionale (fusione 2D continua, di tipo KTHNY; cap. 5.5). In 3D la stessa regola di rigidità sana il modello (cap. 7.9/7.10). Se un secondo ordine locale (cap. 6.4) possa ulteriormente migliorare l'immagine 2D resta una questione di estensione, non un difetto dell'immagine ideale. Poiché il modello deve servire all'illustrazione, accetto la debolezza e rinuncio a un'estensione assiomatica del modello.

<a id="kap7-10"></a>

7.10 Vera scansione di percolazione di legami in 3D (FCC/BCC) e sensibilità del bordo della finestra

I cap. 7.8/7.9 trasferivano concettualmente il valore di connettività 2D z̄_c = 2,08 al 3D (dove v_solid = 1,92, non fisico) e appoggiavano l'affermazione 3D sulla rigidità generica z̄_c = 2d. La versione 13 verifica ora la connettività in 3D direttamente: percolazione di legami su FCC (z_max = 12) e BCC (z_max = 8) per Monte Carlo con union-find, bordi periodici, seme LCG deterministico (rigidity_3d_14.py). Soglia p_c = incrocio P_∞(p) = 0,5 (la stessa metodologia dello studio di dimensione finita 2D del cap. 5.6):

A causa dello spigolo di scatola moderato (L = 12, Python puro) i valori MC stanno leggermente sopra le soglie di letteratura (bias di dimensione finita, come in 2D dal basso o dall'alto a seconda dello stimatore); la graduatoria e l'ordine di grandezza sono corretti (FCC più denso → p_c minore del BCC).

Lettura (conferma 3D indipendente). Le soglie di connettività 3D stanno a z̄_c(3D) ≈ 1,5–1,6 — addirittura più basse del valore 2D 2,08 e ben sotto la rigidità di Maxwell 2d = 6. La «densità di fusione di connettività» che ne segue, v_solid = 1/(3 p_c) ≈ 1,7 (BCC) ovvero 2,6 (FCC), sta ben fuori dalla cupola di coesistenza (v_liq,max ≈ 0,94). Con ciò è mostrato direttamente — senza trasferire il valore 2D al 3D — che in 3D la connettività a maggior ragione non può fissare la densità di fusione; essa deve venire dalla rigidità (z̄_c = 6 → v_solid = 0,667). Le due soglie sono ancora più nettamente separate in 3D che in 2D. Il vero τ_b(3D) (temperatura alla quale z̄(τ) raggiunge la bassa soglia di connettività 1,5) starebbe dunque sopra il valore 2D τ_b, molto vicino al punto critico — coerentemente con l'immagine «disgregazione della connettività = ebollizione».

Sensibilità del bordo superiore stretto della finestra (τ_c − τ_triple ≈ 0,18). Numericamente attorno a v_solid = 0,667 (±2 %) si ottiene dτ_triple/dv_solid ≈ 1,55. La finestra si chiuderebbe (τ_triple = τ_c = 3,00) solo se v_solid crescesse di ≈ 0,11 (fino a ≈ 0,78), cioè se la soglia di rigidità 3D efficace stesse sensibilmente sopra 2d = 6 (z̄_c ≳ 7,7). Per il conteggio generico di Maxwell (z̄_c = 2d = 6) l'ordine τ_m < τ_triple < τ_c resta dunque robusto; il bordo superiore è realmente stretto, ma non a rischio. Il cambio discreto di reticolo solido(1)→solido(2) (FCC↔BCC) è quantificato nella versione 14 (cap. 6.3); la vera scansione in τ della cinetica di legame in 3D segue al cap. 7.11.

<a id="kap7-11"></a>

7.11 Vera scansione in τ della cinetica di legame in 3D: τ_b(3D) e la linea di fusione su una scala ε

I cap. 7.8–7.10 appoggiavano l'affermazione 3D sulla rigidità generica (z̄_c = 2d = 6) e sulla scansione statica di connettività (z̄_c(3D) ≈ 1,55). Restava aperta la vera scansione in τ della cinetica di legame in 3D: dove, sull'asse delle temperature, il z̄(τ) autoconsistente raggiunge questa bassa soglia di connettività (= il vero τ_b(3D)), e dove sta il picco di fusione? La versione 14 calcola il z̄(τ) autoconsistente (cap. 4.2, x = 1/(1 + G·e^{−(ε+K·x)/τ})) direttamente su FCC (z_max = 12) e BCC (z_max = 8) con G = z_max e la stessa cooperatività J = 0,55 del caso di riferimento 2D; n_blk(τ) proviene da un Monte Carlo di legami FCC/BCC (union-find), e

U/N = τ·[3/2 + z̄/4 + (3/2) n_blk] − (ε/2)·z̄ − (K/2)·x² (3D, con termine cooperativo),

c_v = dU/dτ (bondkinetics_3d_15.py, seme LCG deterministico; n_blk lisciato, poiché solo il debole termine di rotazione di blocco ne dipende). Vengono lette tre marche: il picco di c_v τ_m (fusione), la marca di rigidità τ_rig (z̄(τ) attraversa 2d = 6) e il vero τ_b(3D) (z̄(τ) attraversa la soglia di connettività z̄_c(3D) del cap. 7.10):

Risultato. Due constatazioni, entrambe coerenti con l'immagine finora:

1. La fusione sta in basso ed è doppiamente documentata. Il picco termodinamico di fusione τ_m (FCC 1,70, BCC 1,40) coincide praticamente con la marca di rigidità τ_rig (z̄ = 6; FCC 1,73, BCC 1,35) — la stessa temperatura alla quale la rete perde la propria rigidità. Ciò ancora la linea di fusione 3D sulla scala ε: T_m = (ε/k_B)·τ_m, dunque per esempio per l'argon (T_m = 83,8 K) ε/k_B = 83,8/1,70 ≈ 49 K. Poiché K = J·z_max è maggiore in 3D che in 2D (6,60 invece di 3,30), il picco di fusione 3D sta più in alto del valore 2D 1,20 — nessuna contraddizione, ma la cooperatività più forte dell'impacchettamento più denso.
2. Il vero τ_b(3D) siede al punto critico. Poiché la soglia di connettività 3D è così bassa (≈ 1,55), z̄(τ) la raggiunge solo proprio in alto: τ_b(3D) = 3,18 (FCC) ovvero 2,85 (BCC) — al livello di, o leggermente sopra, τ_c = z_max/4 (3,00 ovvero 2,00). L'aggregato gigante di legami sopravvive dunque all'intera fase liquida e si disgrega solo al punto critico. È esattamente la previsione del cap. 7.10 («il vero τ_b(3D) starebbe sopra il valore 2D, molto vicino al punto critico») e conferma in modo indipendente l'immagine «ebollizione = disgregazione della connettività»: in 3D la materia resta legata da parte a parte fin quasi al punto critico.

Entrambe le constatazioni corroborano la separazione delle soglie del cap. 5.4a: fusione ↔ rigidità (basso, τ_m ≈ τ_rig), ebollizione ↔ connettività (alto, τ_b(3D) ≈ τ_c). La densità di fusione viene dunque dalla rigidità (v_solid = 0,667), non dalla connettività — le due soglie sono divaricate sull'asse delle temperature in 3D più che mai prima nel modello.

<a id="kap7-12"></a>

7.12 La linea di trasformazione solido(1)→solido(2) nella superficie (τ, v)

Il cap. 6.3 forniva la trasformazione di reticolo FCC↔BCC come punto τ_ss (a p = 0). Attraverso la pressione ne risulta una linea (polymorph_line_16.py), che appartiene come quarta linea di confine alla superficie (τ, v) comune dei cap. 7.5/7.9 — accanto alla cupola di coesistenza (ebollizione), alla linea di fusione e alla parete dura.

Salto di volume. I due solidi siedono a densità quasi fissa; la loro differenza di volume segue dalle frazioni di impacchettamento (v ∝ 1/η): η_FCC = 0,7405, η_BCC = 0,6802. Con v_FCC = 0,667 (cap. 7.9) si ottiene v_BCC = v_FCC·η_FCC/η_BCC = 0,726, dunque

Δv = v_BCC − v_FCC = +0,059 (+8,9 %).

Il BCC più aperto ha il volume maggiore — la ricristallizzazione FCC→BCC è legata a un aumento di volume (lo stesso enunciato della striscia a fumetti «ricristallizzazione»).

Clausius–Clapeyron. Con ΔS = S_BCC − S_FCC (fase di alta meno fase di bassa temperatura), la pendenza recita

dp/dτ = ΔS/Δv.

Entrambe le grandezze sono positive (il BCC è la fase di alta temperatura, ha dunque più entropia e il volume maggiore), quindi dp/dτ > 0: una pressione maggiore innalza la temperatura di trasformazione, perché favorisce il FCC più denso. Per il riferimento θ_BCC/θ_FCC = 0,60 (τ_ss = 1,333, ΔS = 0,50 k_B, L = 0,667 ε):

Linea in (τ, p) per il riferimento (crescente): p = 0 → τ_ss = 1,33; p = 2 → 1,57; p = 4 → 1,81; p = 6 → 2,04 (unità ridotte).

Nella superficie (τ, v) la trasformazione appare come due bande di densità solida quasi verticali — solido(1) = FCC a v ≈ 0,667 (più denso, T bassa) e solido(2) = BCC a v ≈ 0,726 (più aperto, T alta) — che si convertono l'una nell'altra a τ_ss; sopra segue, a τ_m, la fusione. L'ordine completo sull'asse τ è dunque τ_ss < τ_m < τ_c ≤ τ_b(3D) (cap. 7.11).

Unità reali. Mediante la calibrazione al punto di fusione ε/k_B = T_m/τ_m (τ_m = 1,70, FCC 3D) il gradino di c_v e il calore di trasformazione diventano tangibili. Il gradino di c_v è indipendente dalla calibrazione: |Δc_v| = 1,0 k_B = 1 R = 8,31 J/(mol·K) = 0,0083 kJ/(mol·K). Il calore di trasformazione ΔH_ss = L·R·(ε/k_B) (riferimento L = 0,667 ε):

I calori di trasformazione sono piccoli e positivi (ordine 0,2–0,9 kJ/mol) e stanno dunque nettamente sotto le entalpie di fusione (cap. 7.6, diversi kJ/mol) — esattamente l'immagine attesa di una trasformazione polimorfica debole con un piccolo gradino di c_v. (Il rapporto θ di 0,60 sta deliberatamente fuori dall'immagine ideale rigorosa; il cap. 6.6 mostra che la pura geometria darebbe θ ≈ 0,82–0,94 e la trasformazione resterebbe allora chiusa — la linea è la conseguenza se-allora di una fase BCC molle.)

<a id="kap7-13"></a>

7.13 Linea di trasformazione solido(1)→solido(2) nel diagramma (τ, p), calibrata sul ferro

Il cap. 7.12 porta la trasformazione come linea (τ, v); la versione 16 la aggiunge come linea nell'immagine (τ, p) e la calibra su un polimorfo FCC↔BCC reale (solidsolid_pt_17.py).

Campo di pressione dall'EOS di van der Waals con nocciolo duro. Alle due densità solide v_FCC = 0,667 e v_BCC = 0,726 (cap. 7.12) l'EOS ridotta p_r = 8τ_r/(3v_r − 1) − 3/v_r² (cap. 7.3) dà ovunque pressioni positive; a parità di τ_r il FCC più denso porta la pressione maggiore (τ_r = 1: p_r = 1,25 contro 1,10). La compressione favorisce dunque il FCC — coerentemente con la pendenza positiva di Clausius–Clapeyron dp/dτ = ΔS/Δv > 0.

Linea e punto di trasformazione. Per il riferimento θ_BCC/θ_FCC = 0,60 (τ_ss = 1,333, ΔS = 0,50 k_B, dp\/dτ = 8,47 in unità ridotte p\ = p·D³/ε) la linea cresce come p\(τ) = 8,47·(τ − 1,333): p\ = 0 → τ = 1,33; p\ = 2 → 1,57; p\ = 4 → 1,81; p\ = 6 → 2,04. A **p\_tri ≈ 3,1 la linea solido-solido incontra la marca di fusione τ_m = 1,70 — un punto di trasformazione solido(1)-solido(2)-liquido** nel piano (τ, p) (analogo del punto triplo).

Calibrazione sul ferro (ε/k_B = T_m/τ_m = 1811 K / 1,70 = 1065 K). Il ferro ha due transizioni FCC↔BCC di segno opposto — proprio il caso di prova:

Per γ→δ (τ_ss = 1667/1065 = 1,565 → θ_BCC/θ_FCC ≈ 0,61) il modello coglie segno e ordine di grandezza: ΔH_ss ≈ 4,2 kJ/mol (in realtà ≈ 0,84 kJ/mol — l'ordine di grandezza kJ/mol, con il valore circa 5× troppo alto per l'approssimazione di soglia netta), gradino di c_v 1 R = 8,31 J/(mol·K) senza parametri, dT/dp ≈ +245 K/GPa. La transizione α→γ invece ha la fase densa (γ) come fase di alta temperatura — cosa che la pura geometria d'impacchettamento non può fornire (più denso = più basso = T bassa); il ferro α è stabilizzato magneticamente. Il ferro conferma così entrambi i lati della constatazione onesta: la topologia geometrica (γ→δ) sta dentro l'immagine ideale, l'anomalia magnetica (α→γ) fuori — lo stesso «fuori» che i cap. 6.6/6.8 quantificano con θ ≈ 0,87…1,00 > 0,61.

<a id="kap7-14"></a>

7.14 Redlich–Kwong come ramo separato dell'immagine ideale

La distanza di Guldberg del ramo di van der Waals (cap. 7.4, T_b/T_c ≈ 0,48 invece di ≈ 0,60) era lì soltanto nominata. Su esplicita richiesta dell'utente, una versione migliorata porta la coesione di Redlich–Kwong come variante autonoma dell'immagine ideale con calibrazione propria, coesistenza propria e sezione d'errore propria (rk_branch_17.py) — il ramo di van der Waals resta l'impostazione predefinita; tutta la restante fisica del modello (fusione, percolazione, solido(1)→solido(2)) è indifferente alla scelta dell'EOS di coesione.

Scala di energia propria. Da a_RK = Ω_a·R²T_c^{2,5}/p_c e b_RK = Ω_b·R T_c/p_c segue l'energia di coesione del ramo:

ε_lg,RK/k_B = (Ω_a/Ω_b)·T_c = 4,93·T_c (contro vdW ε_lg/k_B = 27/8·T_c = 3,375·T_c),

maggiore, perché la coesione dipendente dalla temperatura a/√T ripiega il fattore √T_c. Il diametro che segue da b_RK coglie i valori reali in modo simile al ramo vdW (argon D ≈ 3,34·10⁻¹⁰ m).

Coesistenza propria (ridotta, punto critico (1, 1, 1), Z_c = 1/3; p_r = 3τ_r/(v_r − β) − A/(√τ_r·v_r(v_r + β)), β = 0,25992, A = 3,84732):

Sezione d'errore propria (punto di ebollizione a p_r,atm = p_atm/p_c, T_b/T_c = τ_r,b):

Risultato. Il ramo RK innalza T_b/T_c da vdW ≈ 0,48 a ≈ 0,60 e coglie con ciò praticamente la costante di Guldberg e i valori reali (errore medio −1 % contro vdW −21 %). Esso sta come seconda variante dell'immagine ideale, chiaramente separata, accanto all'impostazione predefinita di van der Waals; quale ramo sia «più giusto» è questione della precisione richiesta, non della struttura del modello — il nucleo del modello 3 (z̄ come parametro d'ordine) resta identico in entrambi.

<a id="kap7-15"></a>

7.15 solido(1)→solido(2) su altri polimorfi: titanio e controllo con gas nobile

Il cap. 7.13 calibra la linea (τ, p) solido(1)→solido(2) sul ferro. La versione migliorata assicura segno e ordine di grandezza su un secondo polimorfo, magneticamente non appariscente, e aggiunge un controllo negativo (solidsolid_titanium_18.py). Calibrazione e salto di volume sono come al cap. 7.13 (ε/k_B = T_m/τ_m, τ_m = 1,70; Δv = +8,9 % denso→aperto dalla geometria d'impacchettamento); l'hcp è approssimato come l'FCC tramite z = 12 (la differenza fcc/hcp è solo nominata, non risolta).

Il titanio, con α(hcp) → β(bcc) a 1155 K, ha esattamente la topologia del modello: la fase densa (hcp, z = 12) è quella di bassa temperatura e riscaldando diventa aperta (bcc, z = 8), dunque ΔS > 0 e Δv > 0 → dp/dτ > 0 — lo stesso segno del ferro γ→δ, ma senza anomalia magnetica. Due polimorfi denso→aperto indipendenti a confronto:

Entrambi colgono il segno (dp/dτ > 0, la pressione favorisce la fase densa) e l'ordine di grandezza (ΔH ~ kJ/mol, gradino di c_v senza parametri 1 R = 8,31 J/(mol·K), dT/dp ~ 10²·K/GPa). Il valore di ΔH_ss è in entrambi i casi sistematicamente circa un fattore troppo alto — l'approssimazione di soglia netta del confronto dei μ sovrastima il calore latente di trasformazione; è un limite del modello nominato, non un errore di segno.

L'argon è il controllo negativo: l'argon solido è FCC e resta FCC fino alla fusione (83,8 K) — a pressione normale non esiste bcc stabile. Proprio questo prevede l'immagine ideale del puro potenziale di coppia: θ_BCC/θ_FCC ≈ 0,9…1,0 > 0,61 e un BCC LJ instabile a taglio (cap. 6.9) ⇒ nessuna transizione solido(1)→solido(2). La controprova conferma la struttura della constatazione: l'immagine ideale genera una linea solido(1)→solido(2) solo quando contributi non geometrici (elettronici/magnetici, come nel ferro e nel titanio) stabilizzano in primo luogo il BCC. Il titanio fornisce la prova positiva, l'argon quella negativa dello stesso quadro onesto dei cap. 6.8/6.9.

<a id="kap8"></a>

8 — Confronto bibliografico: concezioni della termodinamica statica

Riassunto. Il capitolo 8 colloca il modello 3 fra le concezioni consolidate della termodinamica statica:

• 8.1 Tabella d'insieme: undici modelli (gas ideale, van der Waals, Redlich–Kwong, sviluppo del viriale, Einstein, Debye, Lennard-Jones / sfere dure, teoria cellulare, Ising / gas reticolare, percolazione, Flory–Huggins) con la loro fase o fasi e la loro idea di base statistico-termodinamica.
• 8.2 I modelli uno per uno: un paragrafo ciascuno con caratteristiche, applicazione e limiti, iniziando dal link alla fonte; le fonti liberamente accessibili sono archiviate sotto model3/links/, mentre la teoria cellulare (a pagamento) è solo citata.
• 8.3 Confronto: il modello 3 è un modello di legame cooperativo a due stati con il parametro d'ordine z̄; il parente più stretto è Ising / gas reticolare, l'ebollizione è ripresa direttamente dalla percolazione e l'immagine dei fononi (Einstein/Debye) sostiene la rigidità di fusione. Il modello 3 tiene insieme queste tre tradizioni su una scala ε comune.

Questo capitolo colloca il modello 3 fra le immagini consolidate della termodinamica di equilibrio («statica») di solidi, liquidi e gas. Prima un panorama tabellare (8.1), poi un breve paragrafo in prosa per modello con caratteristiche tipiche, ambito di applicazione e limiti — ogni paragrafo inizia con il link alla fonte (8.2). Infine un confronto per punti con il modello 3 (8.3). Le fonti liberamente accessibili sono archiviate sotto model3/links/; lo script di verifica model3/scripts/check_links_16.py mantiene coerenti citazione e archivio tramite SHA-256.

<a id="kap8-1"></a>

8.1 Tabella d'insieme

<a id="kap8-2"></a>

8.2 I modelli uno per uno

https://de.wikipedia.org/wiki/Kinetische_Gastheorie — Gas ideale / teoria cinetica dei gas. Particelle senza volume proprio e senza interazione tranne gli urti elastici; la pressione nasce puramente dal trasferimento di quantità di moto. Fornisce l'equazione di stato pV = Nk_BT, la distribuzione delle velocità di Maxwell–Boltzmann e ⟨E⟩ = (3/2)k_BT. Applicazione: gas diluiti ad alta temperatura e bassa pressione. Limiti: nessuna condensazione, nessuna fase liquida o solida, nessuna transizione di fase — l'interazione è assente per costruzione.

https://de.wikipedia.org/wiki/Van-der-Waals-Gleichung — Equazione di van der Waals. (p + a/v²)(v − b) = k_BT completa il gas ideale con un volume di esclusione b (repulsione) e un termine di coesione a (attrazione). È la prima equazione di stato con una transizione gas-liquido, un punto critico e una costruzione di Maxwell. Applicazione: immagine qualitativa della condensazione, grandezze critiche. Limiti: quantitativamente imprecisa vicino al punto critico (esponenti di campo medio), nessuna fase solida.

https://de.wikipedia.org/wiki/Redlich-Kwong-Gleichung — Equazione di stato di Redlich–Kwong. EOS cubica con termine di attrazione dipendente dalla temperatura a/[√T · v(v + b)]. Lo scaling in √T migliora nettamente gli equilibri vapore-liquido rispetto a van der Waals (cfr. cap. 7.4, dove ciò innalza il rapporto T_b/T_c). Applicazione: calcoli tecnici di fase gassosa ed evaporazione. Limiti: parametrizzata empiricamente, tuttora nessuna fase solida; vicino al punto critico parimenti di tipo campo medio.

https://de.wikipedia.org/wiki/Virialentwicklung — Sviluppo del viriale. Sviluppo sistematico pV/Nk_BT = 1 + B(T)/v + C(T)/v² + … I coefficienti del viriale seguono dal potenziale intermolecolare (integrali di cluster della meccanica statistica) e sono dunque microscopicamente fondati. Applicazione: gas fino a densità media, riferimento per altre EOS. Limiti: la serie non converge per liquidi densi o solidi; i coefficienti superiori sono computazionalmente onerosi.

https://de.wikipedia.org/wiki/Einstein-Modell — Solido di Einstein. Il cristallo è un insieme di 3N oscillatori armonici indipendenti di uguale frequenza ω. Spiega la caduta della capacità termica c_v sotto la temperatura caratteristica θ_E = ħω/k_B e il limite di Dulong–Petit 3Nk_B ad alta temperatura. Applicazione: la descrizione più semplice delle capacità termiche dei solidi. Limiti: legge errata a bassa temperatura (esponenziale invece di T³), nessun modo acustico, solo la fase solida.

https://de.wikipedia.org/wiki/Debye-Modell — Modello di Debye. Tratta il solido come continuo elastico con uno spettro fononico fino alla frequenza di Debye; c_v ∝ T³ a bassa temperatura, Dulong–Petit ad alta — la temperatura di Debye θ_D riassume la durezza del reticolo. Applicazione: capacità termiche quantitativamente buone di cristalli semplici su tutto l'intervallo di T. Limiti: approssimazione di continuo (ignora la dispersione dei reticoli reali), solo la fase solida, nessuna fusione.

https://de.wikipedia.org/wiki/Lennard-Jones-Potential — Potenziale di Lennard-Jones / sfere dure. Potenziale di coppia 4ε[(σ/r)¹² − (σ/r)⁶] con nocciolo repulsivo e debole attrazione; il modello limite delle sfere dure mantiene solo la repulsione. È la base della dinamica molecolare e del Monte Carlo e descrive — numericamente — tutte e tre le fasi insieme al punto triplo. Applicazione: gas nobili, sistema di riferimento generico. Limiti: additività a coppie, valutabile nel fondo solo per simulazione o teoria delle perturbazioni (nessuna EOS chiusa).

https://doi.org/10.1098/rspa.1937.0210 — Teoria cellulare del liquido (Lennard-Jones–Devonshire). Ogni molecola si muove in una cella (gabbia) formata dai suoi vicini in campo mediato; il liquido è descritto dalla funzione di partizione cellulare e fornisce una EOS del liquido nonché un approccio alla fusione. Applicazione: liquidi semplici, anello di congiunzione fra solido e gas. Limiti: sovrastima l'ordine, trascura lacune/buchi (vi rimediano teorie dei buchi successive). (A pagamento — solo citata, non archiviata localmente; raggiungibilità verificata nello script di controllo.)

https://de.wikipedia.org/wiki/Ising-Modell — Modello di Ising / gas reticolare. Su un reticolo i siti portano una variabile a due stati (spin ↑/↓ ovvero occupato/vuoto) con accoppiamento ai primi vicini. Il gas reticolare è matematicamente isomorfo al modello di Ising e descrive la condensazione gas↔liquido come transizione cooperativa ordine-disordine con punto critico. Applicazione: fenomeni critici, transizioni cooperative. Limiti: discretizzazione reticolare, nessuna energia cinetica esplicita; in approssimazione di campo medio solo qualitativo.

https://de.wikipedia.org/wiki/Perkolationstheorie — Teoria della percolazione. Modello puramente geometrico: siti o legami sono occupati a caso con probabilità p; oltre una soglia p_c esiste un aggregato che attraversa. Descrive la gelificazione, l'insorgere della conducibilità e in generale le transizioni di connettività. Applicazione: formazione di reti, disgregazione/costituzione della connessione. Limiti: senza energia né temperatura intrinseche — l'accoppiamento alla termodinamica va aggiunto dall'esterno.

https://de.wikipedia.org/wiki/Flory-Huggins-Theorie — Teoria di Flory–Huggins. Modello reticolare per soluzioni e polimeri: l'entropia di mescolamento della collocazione delle catene più un parametro d'interazione χ governano miscibilità e smescolamento. Applicazione: soluzioni polimeriche, diagrammi di fase di miscele. Limiti: campo medio, reticolo incomprimibile, χ per lo più empirico — di rilevanza solo limitata per immagini solido/liquido/gas a un componente, qui tenuta come riferimento reticolare.

<a id="kap8-3"></a>

8.3 Somiglianze e differenze rispetto al modello 3

Il modello 3 è un modello di legame cooperativo a due stati: il parametro d'ordine è il numero di coordinazione medio z̄ = z_max·x, autoconsistente tramite la cooperatività K = J·z_max con la degenerazione dello stato rotto G = z_max; la fusione come soglia di rigidità (fononi), l'ebollizione come disgregazione della connettività (percolazione) al punto critico, transizione del primo ordine non appena ln G > 2 + 4ε/K.

In breve: il modello 3 tiene insieme tre tradizioni altrimenti separate — l'immagine del gas reticolare / di Ising (transizione cooperativa a due stati), l'immagine della percolazione (connettività → ebollizione) e l'immagine dei fononi di Einstein/Debye (rigidità → fusione) — e le calibra mediante un unico parametro di energia ε su una scala τ comune. Van der Waals e la teoria cellulare forniscono la controprova continua (EOS e gabbia liquida), rispetto alla quale l'immagine discreta dei legami può essere confrontata.

<a id="kap9"></a>

9 — Appendice

Riassunto. L'appendice raccoglie le prove formali:

• 9.1 Bibliografia: ordinata per nome d'autore, con link alle fonti liberamente archiviate (confronto SHA-256 tramite check_links_16.py; la teoria cellulare, per via del paywall, è solo citata).
• 9.2 Supporto dell'IA: la quota d'innovazione IA/autore stimata dall'IA per ogni modifica sostanziale (con una ripartizione della quota dell'autore per strategia/struttura/redazione e un calcolo di controllo), nonché un elenco dei singoli passi con costo (stimato), durata e distanza dal registro precedente (somma ≈ 160 USD, stime chiaramente segnalate).
• 9.3 Glossario di abbreviazioni e simboli: i simboli di uso frequente (τ, ε, z̄, K, G, θ, …), che perciò possono comparire in via eccezionale nel testo corrente senza descrizione preposta.
• 9.4 Script di calcolo e copertura dei test: tutte le cifre provengono da script Python riproducibili (quota di programmazione IA 100 %, copertura dei test 0 % ovunque), insieme alla aperta «dichiarazione imbarazzante» che l'autore non ha verificato di persona codice e risultati.

<a id="kap9-1"></a>

9.1 Bibliografia

Ordinata per nome d'autore. Le fonti liberamente accessibili sono archiviate sotto model3/links/; la loro corrispondenza con i link citati è verificata da model3/scripts/check_links_16.py (SHA-256). La teoria cellulare di Lennard-Jones e Devonshire, per via del paywall, è solo citata (raggiungibilità verificata nello script di controllo).

• Broadbent, S. R.; Hammersley, J. M. (1957): Percolation processes I. Crystals and mazes. Proc. Cambridge Phil. Soc. 53, 629–641. — Panoramica/archivio: Teoria della percolazione
• Clausius, R.; Maxwell, J. C.; Boltzmann, L. (1857–1872): fondazione della teoria cinetica dei gas. — Panoramica/archivio: Teoria cinetica dei gas
• Debye, P. (1912): Zur Theorie der spezifischen Wärmen. Ann. Phys. 39, 789–839. — Panoramica/archivio: Modello di Debye
• Einstein, A. (1907): Die Plancksche Theorie der Strahlung und die Theorie der spezifischen Wärme. Ann. Phys. 22, 180–190. — Panoramica/archivio: Modello di Einstein
• Flory, P. J.; Huggins, M. L. (1941–1942): teoria reticolare delle soluzioni polimeriche (parametro χ). — Panoramica/archivio: Teoria di Flory–Huggins
• Ising, E.; Lenz, W. (1920–1925): Beitrag zur Theorie des Ferromagnetismus. Z. Phys. 31, 253–258 (il gas reticolare è matematicamente isomorfo). — Panoramica/archivio: Modello di Ising
• Kamerlingh Onnes, H. (1901): forma del viriale dell'equazione di stato. — Panoramica/archivio: Sviluppo del viriale
• Lennard-Jones, J. E. (1924): On the determination of molecular fields. Proc. R. Soc. Lond. A 106, 463–477. — Panoramica/archivio: Potenziale di Lennard-Jones
• Lennard-Jones, J. E.; Devonshire, A. F. (1937): Critical phenomena in gases I. Proc. R. Soc. Lond. A 163, 53–70. https://doi.org/10.1098/rspa.1937.0210 (A pagamento — solo citato, non archiviato localmente.)
• Redlich, O.; Kwong, J. N. S. (1949): On the thermodynamics of solutions. V. An equation of state. Chem. Rev. 44, 233–244. — Panoramica/archivio: Equazione di Redlich–Kwong
• van der Waals, J. D. (1873): Over de Continuïteit van den Gas- en Vloeistoftoestand (tesi, Leida). — Panoramica/archivio: Equazione di van der Waals

<a id="kap9-2"></a>

9.2 Supporto dell'IA

Questo sviluppo del modello è nato da un'interazione fra uomo e IA: l'autore (l'utente dell'IA) stabilisce su che cosa si lavora (domande guida, decisioni, desideri sul modello); l'IA (Claude Code) stabilisce in larga misura come — derivazioni, script, struttura e testo. Le due tabelle che seguono sono stimate dall'IA stessa.

(a) Modifiche sostanziali di contenuto e quota d'innovazione stimata. La quota d'innovazione pesa le idee guida al 45 %, i lavori strutturanti al 40 %, i lavori redazionali al 15 %. Le idee guida (il livello del su che cosa) provengono prevalentemente dall'autore, i lavori strutturanti e redazionali (il livello del come) prevalentemente dall'IA; là dove l'IA ha apportato anche l'idea fisica guida (per esempio la lettura per percolazione dell'ebollizione, il modello a molle per θ), la sua quota cresce. Valori in percentuale, IA + autore = 100. Le tre colonne aggiuntive ripartiscono la quota dell'autore secondo i tre livelli d'azione — strategia (idee guida), struttura (derivazioni/script/architettura) e redazione (lavoro sul testo). La loro media pesata al 45 % / 40 % / 15 % riproduce la colonna Autore % come calcolo di controllo (model3/scripts/author_share_19.py, scarto massimo 0,95 punti percentuali).

(b) Singoli passi: costo (stimato), durata, distanza dal registro precedente. Data di creazione è l'ora di creazione del file di registro (mtime del file). Δ al registro precedente è la distanza temporale dalla voce di registro precedente (misura dell'impegno dell'utente fra i passi, incluse pause e /compact). Costo (stimato) sono le stime documentate per registro (prezzi di listino di Anthropic, ± 30 %). Durata (attiva) è una stima dell'IA del tempo di lavoro puro; per passi separati da /compact può essere maggiore di Δ. Un asterisco segnala una presunta ricostruzione o rielaborazione di contenuti esistenti.

Somma dei costi documentati (stimati, passi 01–21): ≈ 129 USD (somma delle singole stime); con il passo 22 (16d, ≈ 10 USD) ≈ 139 USD, il passo 23 (V16, ≈ 12 USD) ≈ 151 USD e il passo 24 (V17, ≈ 9 USD) ≈ 160 USD. Tutti i valori sono stime chiaramente segnalate (nessuna misurazione ccusage disponibile); i dettagli dei registri stanno in model3/logs/dev_log_NN.md.

<a id="kap9-3"></a>

9.3 Glossario di abbreviazioni e simboli

Per principio questo testo antepone a ogni simbolo un termine descrittivo (dunque «capacità termica c_v» invece del solo «c_v»). Le seguenti grandezze di uso frequente costituiscono l'eccezione: possono comparire nel testo corrente senza descrizione preposta, perché qui sono spiegate una volta per tutte. Tutte le temperature τ… sono ridotte (in unità dell'energia di legame ε).

<a id="kap9-4"></a>

9.4 Script di calcolo e copertura dei test

Tutte le cifre e le tabelle nominate nel testo provengono da script Python riproducibili in model3/scripts/ (ciascuno con un output grezzo …_output.txt). La quota di programmazione IA indica quale parte del codice degli script è stata scritta dall'IA; la copertura dei test riguarda i test automatizzati (test unitari, controlli assert). Essa è ovunque bassa: gli script sono calcoli unici senza impalcatura di test; la loro unica garanzia è la plausibilità dell'output grezzo.

L'autore confessa apertamente (dichiarazione imbarazzante): non ha verificato, ricalcolato né controllato sul piano specialistico né il codice sorgente di questi script né le cifre, le tabelle e le curve da essi prodotte. Si è qui affidato ciecamente alla competenza dell'IA, perché in primo piano stavano solo il carattere qualitativo e la plausibilità della concezione del modello. I valori vanno perciò presi con cautela (si veda anche la nota dettagliata al cap. 10.4).

<a id="kap10"></a>

10 Autocritica

Ho la sensazione che nella tabella del capitolo 9.2 (a) la quota di lavoro dell'autore, stimata dall'IA, sia in tutti e tre gli ambiti quasi sempre troppo alta di almeno 10–20 punti percentuali. Ho la sensazione di aver contribuito solo con l'idea di base, una o due interpretazioni e importanti raffinamenti nella visualizzazione grafica.

I capitoli (matematici) 4–9 sono stati proposti, ricercati e realizzati dall'IA. Peggio ancora: durante lo sviluppo avevo la sensazione che l'IA fosse fissata sull'elaborazione di una statistica basata sul numero di coordinazione. Le mie vere idee di base — che non solo gli atomi ma anche le molecole possono formare legami e che le rotture di legame all'interno di un liquido potrebbero seguire regolarità frattali — l'IA le ha potute accogliere solo a fatica, cercando di ignorare la terminologia. Perciò quelle idee di base sono ancora trattate in modo molto secondario nel testo. Non l'ho riscritto in seguito, perché riflettendo sulle transizioni di fase a partire dagli assiomi fondamentali si tengono naturalmente presenti queste idee. Inoltre la fissazione, o riduzione, delle idee al numero di coordinazione come grandezza di misura non è la via peggiore per matematizzare la mia rappresentazione per immagini, poiché a posteriori si è rivelata piuttosto portante per la verifica di plausibilità — pur avendo io fatto solo una verifica di plausibilità molto grossolana dei risultati dell'IA e nessun esame del codice Python. A mia discolpa sia detto: ciò che mi importava è sempre stato soltanto una rappresentazione per immagini dei rapporti molecolari nei diversi stati di aggregazione. All'adulazione dell'IA, che voleva confermare matematicamente la rappresentazione, non ho detto di no. La manipolazione da parte dell'IA ha dunque avuto successo su di me, perché ho esaminato i capitoli nominati solo nel modo più superficiale e li ho accolti con fiducia. (Questa esperienza su me stesso rafforza la mia diffidenza sul fatto che un «human in the loop» sia davvero efficace per controllare l'IA.)

Questo testo è stato realizzato con il supporto dell'IA (Claude Code 2.1.161, modello Opus 4.8 / claude-opus-4-8[1m], contesto 1M; data 2026-09-05). — Realizzato con il supporto dell'IA.

Nota sulla filigrana. Il testo tedesco di riferimento (de.modell_final_02.md) porta sulle parole chiave una filigrana fatta di omoglifi greci (Α→A, Β→B, Ε→E, Η→H, Ι→I, Κ→K, Μ→M, Ν→N, Ο→O, Ρ→P, Τ→T, Υ→Y, Χ→X, Ζ→Z). Poiché quell'espediente è legato alla grafia tedesca delle parole chiave, questa traduzione porta le parole chiave in maiuscole normali: VOLO LIBERO (FREIFLUG), MOLTITUDINE (VIELHEIT), LEGAME (BINDUNG), TRIADE (DREIKLANG), NOCCIOLO DURO (HARTKERN), CONSERVAZIONE (ERHALTUNG), FRATTURA (RISS), SVELTO (FLINK).