How glass breaks -- Damage explains the difference between surface and fracture energies in amorphous silica

Sintesi Tecnica: Come si rompe il vetro – Il danno spiega la differenza tra energia superficiale ed energia di frattura nella silice amorfa

Problema
Nella meccanica della frattura, esiste una discrepanza fondamentale tra l'energia libera superficiale ($\gamma$), che quantifica l'energia necessaria per creare nuove superfici, e l'energia di frattura ($G_c$), che rappresenta l'energia totale dissipata durante la propagazione di una cricca. Nei materiali fragili ideali, $G_c$ dovrebbe essere uguale a $2\gamma$. Tuttavia, le misurazioni sperimentali sui vetri di silicato mostrano costantemente che $G_c$ supera $2\gamma$ di un fattore di circa 5. Questo eccesso di energia viene spesso attribuito alla dissipazione plastica all'apice della cricca, ma tale spiegazione rimane controversa per la silice amorfa. Le prove di plasticità sono difficili da ottenere sperimentalmente su scala atomica, e i modelli di plasticità standard non sono applicabili ai silicati a causa della loro mancanza di incrudimento. Di conseguenza, l'origine fisica di questa discrepanza energetica nei silicati amorfi fragili è rimasta irrisolta.

Metodologia
Gli autori hanno impiegato un approccio di simulazione multiscala combinando la dinamica molecolare (MD) su scala atomica con la modellazione a campo di fase (phase-field) continua per investigare la frattura nella silice amorfa.

• Simulazioni Atomiche: Le simulazioni MD sono state condotte utilizzando il potenziale BKS per modellare le interazioni atomiche e la troncatura di Wolf per le forze coulombiche. I campioni sono stati generati tramite posizionamento casuale degli atomi, equilibrati e raffreddati a $10^{-5}$ K. La propagazione della frattura è stata simulata in condizioni quasi-statiche e atermiche utilizzando una condizione al contorno di spostamento a campo K per guidare la crescita della cricca.
• Coarse-Graining (Raggruppamento): Per colmare il divario tra le scale atomica e continua, è stata applicata una tecnica di convoluzione gaussiana basata sulla fisica. Questo metodo ha generato campi omogenei per spostamento, deformazione, sforzo ed energia potenziale, conservando massa ed energia.
• Quantificazione del Danno e dell'Energia:
  • Danno: Definito meccanicamente come la perdita di rigidezza, calcolato confrontando l'energia elastica reale con l'energia di deformazione non danneggiata.
  • Energia Libera Superficiale (FSE): Isolata calcolando la differenza tra l'energia potenziale totale e l'energia di deformazione elastica ($\psi_{FSE} = \psi_{pot} - \psi_{\epsilon}$), separando efficacemente l'energia di formazione della superficie dalla deformazione elastica.
  • Calibrazione del Campo di Fase: È stato utilizzato uno schema di Aggiornamento di Elementi Finiti (FEMU) per adattare il campo di danno ottenuto dalle simulazioni MD a un modello a campo di fase continuo. Ciò ha permesso di estrarre i parametri del materiale, specificamente la tenacità alla frattura ($g_c$) e la scala di lunghezza interna ($l_c$), direttamente dagli istantanee atomiche.
• Analisi Strutturale: I cambiamenti strutturali locali sono stati quantificati tracciando l'evoluzione delle strutture ad anelli (utilizzando l'algoritmo di Le Roux & Jund) e il numero di coordinazione degli atomi di silicio.

Risultati Chiave

1. Danno Diffuso vs. Plasticità: Lo studio dimosta che la dissipazione di energia attorno all'apice della cricca è dovuta principalmente a un danno diffuso piuttosto che alla plasticità. La zona di danno si estende per circa 16–23 Å dal percorso della cricca. L'analisi della risposta sforzo-deformazione e del bilancio energetico conferma che la dissipazione plastica e l'elasticità non lineare sono trascurabili; il materiale perde rigidezza (danno) invece di mantenerla (plasticità).
2. Separazione dei Contributi Energetici: Gli autori hanno separato con successo l'energia di frattura totale in due componenti distinte:
   • Energia Libera Superficiale: Localizzata strettamente sulla superficie, associata alla rottura dei legami e al rilassamento atomico.
   • Energia di Danno: Estesa più profondamente nel materiale (~20 Å), associata ai riarrangiamenti strutturali.
     Il rapporto tra l'energia libera superficiale e l'energia di frattura totale ($2\gamma/g_c$) è stato calcolato essere circa il 23%, in linea con le osservazioni sperimentali dove l'energia superficiale costituisce solo circa il 20% dell'energia totale.
3. Firme Strutturali:
   • Danno Diffuso: Correla fortemente con i cambiamenti nella struttura ad anelli del reticolo del vetro, che presenta un'estensione spaziale più ampia.
   • Energia Libera Superficiale: Correla con i cambiamenti nel numero di coordinazione degli atomi di silicio, che sono più localizzati vicino alla superficie.
4. Validazione: I valori di tenacità alla frattura derivati dall'adattamento del campo di fase, dal bilancio energetico globale e dal metodo dell'integrale J sono stati coerenti. La scala di lunghezza interna calcolata ($l_c \approx 12\text{--}20$ Å) si allinea con la larghezza osservata della zona di danno.

Significato e Rivendicazioni
L'articolo sostiene di fornire un quadro pienamente coerente per comprendere la frattura nella silice amorfa, risolvendo la storica discrepanza tra energie superficiali e di frattura. Il contributo primario è l'identificazione del danno diffuso — distinto dalla plasticità — come il meccanismo dominante per la dissipazione di energia nei vetri di silicato fragili. Collegando questo danno a specifici riarrangiamenti strutturali (rottura della struttura ad anelli) e quantificandone l'estensione spaziale, lo studio offre una spiegazione strutturale della natura "fragile" ma capace di dissipazione energetica della silice.

Gli autori affermano che il loro approccio valida il concetto di danno diffuso utilizzato nelle tecniche di frattura a campo di fase, fornendo una dimostrazione diretta della sua base fisica su scala atomica. Inoltre, i risultati suggeriscono che le visioni tradizionali della frattura fragile potrebbero trascurare la sottile dissipazione meccanica inerente ai materiali amorfi. Il documento conclude con moderazione, osservando che, sebbene queste scoperte chiariscano il meccanismo nei silicati, sono necessarie ulteriori indagini per determinare se tali concetti siano applicabili ad altri sistemi amorfi, come i vetri metallici, o a scenari di carico dinamico e di fatica.