Percolation Criticality of Amorphous-Amorphous Transitions in Compressed Glasses

Questo studio impiega simulazioni di dinamica molecolare su larga scala e la teoria della percolazione per rivelare che la transizione da bassa ad alta densità nel vetro di silice compresso è guidata dalla percolazione critica di cluster strutturali, esibendo esponenti critici che suggeriscono un meccanismo di percolazione della rigidità e sottolineando principi di trasformazione comuni tra materiali amorfi legati e non legati.

L'essenziale

Immagina un bicchiere d'acqua o un pezzo di vetro di una finestra. Li conosci come oggetti solidi e rigidi. Ma se ti rimpicciolissi fino alle dimensioni di un atomo per guardare all'interno, vedresti una rete caotica e aggrovigliata di minuscoli blocchi costruttivi. Nel vetro di silice (quello dei vetri delle finestre), questi blocchi hanno la forma di piramidi (tetraedri) fatte di silicio e ossigeno.

Questo articolo è come un film ad alta tecnologia che zooma su ciò che accade quando schiacci questo vetro con una pressione immensa — fino a 350.000 volte la pressione dell'atmosfera. Gli scienziati volevano capire come il vetro cambi la sua forma senza fondere o rompersi, un processo chiamato "transizione amorfo-amorfa".

Ecco la storia della loro scoperta, raccontata attraverso semplici analogie:

1. La folla a un concerto (La struttura)

Immagina il vetro come una grande folla di persone a un concerto.

• A pressione normale: Tutti sono in piedi in una formazione lenta e aperta. Nel vetro di silice, le "persone" sono atomi di silicio, e stanno tenendosi per mano con quattro vicini, formando forme piramidali perfette (tetraedri). Ci sono molti spazi vuoti tra di loro, il che rende la struttura "molle" e facile da comprimere.
• All'aumentare della pressione: Immagina che la sala del concerto inizi a restringersi. La folla viene schiacciata. Le persone non riescono più a mantenere le loro forme piramidali perfette. Iniziano a urtarsi l'un l'altra, cambiando il numero di vicini con cui si tengono per mano. Alcuni iniziano a tenersi per mano con 5 persone, poi con 6.

2. Il gioco della "Percolazione" (Il grande cambiamento)

Gli scienziati hanno usato un concetto chiamato percolazione. Immagina di versare acqua attraverso una spugna.

• Bassa pressione: La spugna ha dei buchi, ma sono tutti separati. Se versi l'acqua, questa rimane bloccata in piccole sacche. Non scorre tutto attraverso. Nel vetro, le forme a "piramide" sono isole isolate.
• Pressione critica: Man mano che schiacci più forte, i buchi iniziano a connettersi. Improvvisamente, si forma un percorso gigante e continuo dal top al fondo della spugna. L'acqua scorre!
• Nel vetro: Gli scienziati hanno scoperto che a specifici punti di pressione, le nuove forme (come i blocchi a 5 o 6 lati) si connettono improvvisamente per formare una catena gigante e continua che attraversa l'intero pezzo di vetro. Questa è la "transizione di percolazione". È il momento in cui il vetro si riorganizza fondamentalmente in uno stato più denso.

3. Due modi per guardare la folla

I ricercatori hanno osservato questa folla in due modi diversi, come se guardassero un film da due diverse angolazioni di telecamera:

• La vista "Legata" (La stretta di mano): Hanno guardato chi si sta tenendo direttamente per mano (legami chimici). Hanno visto che le forme piramidali stavano cambiando le loro strette di mano.
• La vista "Non Legata" (Lo spazio personale): Hanno ignorato le strette di mano e si sono limitati a guardare chi stava vicino a chi, indipendentemente dal fatto che si toccassero. È come guardare una folla dove le persone non si tengono per mano, ma stanno solo vicine tra loro.

La Sorpresa: Entrambe le telecamere hanno mostrato esattamente la stessa storia! La vista della "stretta di mano" e la vista dello "spazio personale" hanno mostrato che il vetro si trasforma nella stessa sequenza: prima le forme sciolte si connettono, poi le forme dense prendono il sopravvento. Ciò suggerisce che le regole che governano il modo in cui il vetro cambia sono universali, sia che gli atomi si stiano "tenendo per mano" (come nella silice) sia che si stiano solo urtando tra loro (come nell'acqua congelata/ghiaccio).

4. Il "Numero Magico" e le regole del gioco

Gli scienziati volevano sapere se questa trasformazione segua un manuale di regole standard (come un gioco d'azzardo) o se abbia le sue regole speciali.

• I Tetraedri (Le forme a 4 lati): Quando le forme piramidali originali (che tengono 4 mani) si sono disgregate, lo hanno fatto esattamente come un gioco casuale di probabilità. Era un comportamento "standard".
• Le forme superiori (5, 6 o più mani): Quando le nuove forme, più dense, si sono formate e connesse, hanno infranto le regole standard. Hanno seguito un insieme di regole diverso e più complesso. Gli scattivisti chiamano questo "percolazione di rigidità". È come se la folla non si fosse solo connessa casualmente; si era connessa in un modo che rendeva l'intera struttura improvvisamente molto più rigida e soda.

5. La conclusione

L'articolo conclude che quando schiacci il vetro, non si limita solo a diventare più piccolo; subisce un evento drammatico, simile a un salto di fase, in cui la struttura interna si riorganizza in un nuovo "stato" più denso.

• La transizione avviene a pressioni "critiche" specifiche.
• Il modo in cui le nuove strutture si connettono è un misto di caso casuale (per le vecchie forme) e una regola più rigida e strutturata (per le nuove forme dense).
• Questo comportamento è simile nella silice e nel ghiaccio amorfo, suggerendo che la natura utilizzi "progetti" simili per riorganizzare diversi tipi di materiali vetrosi sotto pressione.

In breve, l'articolo mappa esattamente come lo "scheletro" microscopico del vetro si spezza, si sposta e si ricostruisce quando viene schiacciato, rivelando che la transizione da un vetro sciolto e flessibile a uno denso e rigido avviene attraverso un punto di svolta specifico e prevedibile.

Sommario tecnico

Sintesi Tecnica: Criticità della Percolazione nelle Transizioni Amorfo-Amorfe in Vetri Compressi

Definizione del Problema
I meccanismi che governano le trasformazioni strutturali nei vetri compressi, in particolare le transizioni da amorfo ad amorfo in materiali a base tetraedrica come la silice ($a$-SiO$_2$) e il ghiaccio amorfo ($a$-H$_2$O), rappresentano una sfida significativa. A differenza delle transizioni di fase cristalline nette, queste trasformazioni comportano un'evoluzione progressiva di motivi strutturali locali in distinti poliaomorfi. Sebbene studi precedenti abbiano identificato l'emergere di grandi cluster percolanti di poliedri ad alta coordinazione sotto pressione, la criticità di queste transizioni nei vetri è rimasta ignota. Nello specifico, non è chiaro se queste transizioni strutturali rientrino nella classe di universalità della percolazione standard (casuale) non correlata o se esibiscano un comportamento critico distinto, come la percolazione della rigidità. I tentativi precedenti di stimare gli esponenti critici sono stati limitati dalle piccole dimensioni dei sistemi inerenti agli approcci ab initio.

Metodologia
Per affrontare queste limitazioni, gli autori hanno impiegato simulazioni di dinamica molecolare (MD) classica su larga scala.

• Sistemi e Potenziali: Le simulazioni sono state condotte su vetri di silice utilizzando il potenziale a coppie SHIK nel pacchetto LAMMPS, e su ghiaccio amorfo utilizzando il campo di forza TIP4P in Gromacs. Le dimensioni del sistema variavano da $\sim$1.000 a $10^6$ atomi per la silice e fino a 49.000 molecole per il ghiaccio.
• Protocollo: I campioni sono stati equilibrati, raffreddati a temperatura ambiente (300 K per la silice, 124 K per il ghiaccio) e compressi quasi-staticamente fino a 35 GPa (silice) o 15 kbar (ghiaccio).
• Descrittori Strutturali: Sono stati utilizzati due approcci complementari:
  1. Modello Legato (Bonded): Basato sulla connettività covalente Si-O, analizzando i poliedri SiO$_Z$ (dove $Z$ è il numero di coordinazione) e le loro connettività (condivisione di spigoli, vertici o facce).
  2. Modello Non Legato (Non-bonded): Basato sulla disposizione degli atomi centrali (Si-Si per la silice, O-O per il ghiaccio) entro un raggio di cutoff, definendo i poliedri SiSi$_Z$ e OO$_Z$. Questo approccio tratta le interazioni come forze non legate, permettendo un confronto diretto tra vetri legati (silice) e non legati (ghiaccio).
• Analisi della Percolazione: Un codice Python dedicato (Nexus-CAT) implementando un algoritmo Union-Find è stato utilizzato per identificare i cluster. Le proprietà di percolazione, inclusi la lunghezza di correlazione ($\xi$), l'ordine parametrico ($P_\infty$), la dimensione media del cluster ($\langle S \rangle$) e la dimensione del cluster massimo ($S_{max}$), sono state calcolate.
• Esponenti Critici: Metodi di scaling delle dimensioni finite (utilizzando tecniche di data collapse) sono stati applicati per determinare gli esponenti critici ($\nu, \beta, \gamma$) e le dimensioni frattali ($D_f, D_f^*$) vicino alle soglie di percolazione.

Risultati Chiave

1. Sequenza di Transizioni di Percolazione: Sotto compressione, lo studio ha identificato una distinta sequenza di eventi di percolazione. In $a$-SiO$_2$, le reti tetraedriche a bassa densità (LD, $Z=4$) depercolano mentre emergono e percolano i cluster ad alta densità (HD, $Z=5, 6$). Nello specifico:

   • I cluster SiO$_5$ percolano a $\approx$12 GPa.
   • I cluster SiO$_6$ percolano a $\approx$23 GPa.
   • I cluster "tipo-stishovite" (SiO$_6$ con due connessioni di condivisione di spigoli) percolano a $\approx$30 GPa.
   • Sequenze simili sono state osservate in $a$-H$_2$O (cluster LD, HD e VHD), che avvengono a pressioni circa 20 volte inferiori rispetto alla silice, suggerendo un meccanismo universale scalato dalle interazioni specifiche della sostanza.

2. Complementarità dei Modelli: Gli approcci legato (SiO$_Z$) e non legato (SiSi$_Z$/OO$_Z$) hanno fornito intuizioni coerenti e complementari. Il modello non legato ha catturato con successo le trasformazioni strutturali e ha permesso un quadro unificato per confrontare la silice e il ghiaccio amorfo.

3. Esponenti Critici e Universalità:

   • La lunghezza di correlazione $\xi$ è scalata con la dimensione del sistema $L$ come $\xi \propto L^a$ con $a \approx 1$, confermando il comportamento critico.
   • Depercolazione dei Tetraedri: La depercolazione dei cluster tetraedrici LD ($Z=4$) ha prodotto esponenti critici (ad esempio, l'esponente di Fisher $\tau \approx 2.12-2.14$) molto vicini a quelli della percolazione casuale 3D standard.
   • Percolazione dei Cluster ad Alta Coordinazione: Al contrario, la percolazione dei cluster a coordinazione più elevata ($Z=5, 6$) e delle strutture VHD ha mostrato deviazioni dalla percolazione standard. Le dimensioni frattali ($D_f$) sono state sistematicamente inferiori al valore standard di 2.53, e gli esponenti suggeriscono una classe di universalità differente.

4. Dimensioni Frattali: La dimensione frattale $D_f^*$ dei cluster finiti (animali di reticolo) è stata calcolata per vari numeri di coordinazione, con un intervallo compreso tra $\approx$2.29 e 2.46, generalmente inferiore al valore standard della percolazione.

Significato e Rivendicazioni
Il documento afferma che l'introduzione di descrittori a lungo raggio e l'uso di simulazioni su larga scala permettono la prima stima robusta degli esponenti critici di percolazione nei vetri compressi. La principale significatività risiede nel fatto che le transizioni amorfo-amorfo non sono uniformi nel loro comportamento critico:

• Il collasso della rete tetraedrica nativa sembra seguire la percolazione casuale standard.
• L'emergere e la percolazione di reti poliedriche rigide ad alta coordinazione (associate alla densificazione) seguono probabilmente una classe di universalità diversa, indicativa di un meccanismo di percolazione della rigidità. Ciò è attribuito alle eterogeneità topologiche ed elastiche che derivano quando cluster percolanti di alta connettività emergono all'interno di un mezzo a bassa connettività.

Lo studio stabilisce un legame tra queste transizioni di percolazione e le proprietà meccaniche dei vetri, come il picco di compressibilità e l'inizio della plasticità. Inoltre, propone che il monitoraggio delle trasformazioni strutturali guidate dalla percolazione fornisca una base per analogie tra gli stati amorfi e i poliaomorfi cristallini (ad esempio, coesite IV, coesite V, stishovite), suggerendo che gli stati vetrosi possano corrispondere a megabasini di energia associati a queste strutture cristalline. Il lavoro sottolinea l'utilità dell'approccio non legato per unificare la comprensione delle trasformazioni indotte dalla pressione sia in sistemi vetrosi legati che non legati.