Coarse-graining amorphous plasticity: impact of rejuvenation and disorder

Sintesi Tecnica: Coarse-graining della Plasticità Amorfa: Impatto di Rejuvenation e Disordine

Enunciato del Problema
Il comportamento meccanico dei solidi amorfi, caratterizzato da fenomeni quali l'invecchiamento (aging), il rejuvenation, il banding di taglio (shear banding) e l'effetto Bauschinger, deriva da riarrangiamenti strutturali localizzati (Trasformazioni di Taglio) che interagiscono attraverso un mezzo elastico. Sebbene le simulazioni atomistiche catturino queste caratteristiche, esse sono computazionalmente costose. Al contrario, i modelli Elasto-Plastici (EPM) mesoscopici offrono un framework computazionalmente efficiente, ma spesso si basano su ipotesi eccessivamente semplificate che possono fallire nel riprodurre quantitativamente i risultati atomistici. Il problema centrale affrontato è la sfida del "coarse-graining": determinare come trasferire i parametri dalle scale atomistiche ai modelli a reticolo mesoscopici per riprodurre la fisica della plasticità amorfa, concentrandosi specificamente sul ruolo del sottostante paesaggio di energia potenziale (PEL) disordinato, l'evoluzione del disordine sotto carico e il rejuvenation dello stato vetroso.

Metodologia
Lo studio impiega un approccio comparativo tra simulazioni atomistiche Athermal Quasistatic (AQS) e un modello EPM scalare mesoscopico semplificato.

1. Simulazioni Atomistiche:

   • Sono stati preparati tre distinti campioni di vetro binario Lennard-Jones 2D con storie termiche contrastanti:
     • MCT: Raffreddato istantaneamente da una temperatura vicina alla Transizione di Mode Coupling ($T_p \approx 0.38$).
     • ESL: Raffreddato istantaneamente da una temperatura più alta ($T_p = 0.351$), rappresentando un vetro poco ricotto (soft).
     • GQ: Raffreddato gradualmente da $T_p = 0.351$, rappresentando un vetro ben ricotto (hard).
   • I campioni sono stati sottoposti a carico di taglio AQS fino a una deformazione globale di $\gamma = 5$.
   • Le proprietà locali sono state caratterizzate utilizzando il "metodo della matrice congelata" (frozen matrix method), che isola regioni circolari per misurare gli stress di snervamento locali ($\sigma_c$), le fluttuazioni di stress interno ($\delta\sigma$) e i cali di stress ($\Delta\sigma$) attraverso 18 direzioni di taglio.

2. Modello Mesoscopico:

   • È stato implementato un modello a reticolo scalare su una griglia quadrata con condizioni al contorno bi-periodiche.
   • Dinamica: Il modello presenta una dinamica di soglia locale in cui una cella cede se lo stress locale supera una soglia $\sigma_c$. Gli eventi plastici innescano una ridistribuzione dello stress elastico a lungo raggio (kernel di Eshelby) e un rinnovo della soglia locale.
   • Parametri Chiave: Lo studio si concentra su due parametri adimensionali:
     • Forza del Disordine ($s$): Il rapporto tra le fluttuazioni dello stress di snervamento e lo scatto di stress indotto da un evento plastico.
     • Parametro di Età/Rejuvenation ($r$): Definito come $r = (\sigma_{c}^{iniziale} - \sigma_{c}^{rinnovato}) / \sigma_{c}^{rinnovato}$. Quantifica la differenza tra lo stato iniziale del vetro e lo stato del materiale dopo un riarrangiamento plastico.
   • Strategia di Up-scaling: I parametri per il modello mesoscopico (distribuzioni degli stress di snervamento, incrementi di scivolamento e distribuzioni di rinnovo) sono stati estratti direttamente dai dati atomistici. Lo studio ha testato due ipotesi specifiche:
     • Scenario Disordinato: La distribuzione di rinnovo degli stress di snervamento è identica alla distribuzione iniziale ($r=0$).
     • Scenario di Rejuvenation: La distribuzione di rinnovo corrisponde al vetro MCT (uno stato più "fresco"), implicando $r > 0$.
     • Ipotesi di Correlazione: È stata introdotta una correlazione diretta tra lo stress di snervamento locale e l'ampiezza dell'incremento di scivolamento successivo ($\varepsilon_0 \propto \sigma_c^2$), riflettendo la struttura gerarchica del PEL.

Contributi Chiave e Risultati

1. Comportamento Stress-Deformazione e Transienti:

   • Modelli semplici senza correlazione tra stress di snervamento e incrementi di scivolamento non sono riusciti a riprodurre il picco di stress transiente osservato nelle simulazioni atomistiche, producendo solo curve monotone.
   • L'introduzione di una correlazione tra lo stress di snervamento locale e l'ampiezza dell'incremento di scivolamento ha permesso al modello mesoscopico di riprodurre l'intera curva stress-deformazione, incluso il picco transiente e lo stress di flusso allo stato stazionario, sia per i vetri soft (ESL) che per quelli hard (GQ).
   • Si è scoperto che lo stress di flusso allo stato stazionario è controllato principalmente dal parametro di forza del disordine $s$ (legato all'ampiezza media dello scivolamento), piuttosto che dai dettagli specifici dello scenario di rejuvenation.

2. Rejuvenation e Evoluzione dello Stress di Snervamento Locale:

   • Vetro Soft (ESL): Le simulazioni atomistiche hanno mostrato un leggero indurimento (aumento dello stress di snervamento medio locale) seguito da un plateau. Il modello mesoscopico con rejuvenation ($r \approx 0.14$) ha ottenuto un accordo quasi quantitativo con questa evoluzione. Senza rejuvenation, il modello sovrastimava lo stress di snervamento.
   • Vetro Hard (GQ): Le simulazioni atomistiche hanno mostrato un indurimento iniziale seguito da un regime di addolcimento (softening) associato al banding di taglio. Il modello mesoscopico senza rejuvenation non è riuscito a riprodurre la fase di addolcimento, mostrando solo un continuo indurimento. L'inclusione dello scenario di rejuvenation ($r \approx 0.65$) è stata cruciale per recuperare qualitativamente la diminuzione dello stress di snervamento medio, sebbene l'accordo quantitativo rimanesse imperfetto.

3. Localizzazione e Shear Banding:

   • Per il vetro hard (GQ), le simulazioni atomistiche hanno mostito un chiaro banding di taglio.
   • Il modello mesoscopico con rejuvenation e correlazione ha riprodotto con successo l'evoluzione dell'indice di localizzazione: un rapido aumento durante la nucleazione della banda seguito da una graduale diminuzione man mano che la banda si allarga.
   • I modelli senza rejuvenation hanno mostrato una localizzazione significativa, ma hanno sottostimato l'intensità e la persistenza della banda di taglio.

4. Fluttuazioni di Deformazione e Stress:

   • Sebbene il modello mesoscopico abbia riprodotto la tendenza dell'evoluzione della varianza della deformazione, ha sistematicamente sovrastimato l'entità delle fluttuazioni di un fattore di circa 1,5 rispetto ai risultati atomistici.
   • Le fluttuazioni del campo di stress interno sono state sovrastimate ancora più significativamente (da un fattore di 2 a 3). Gli autori attribuiscono questa discrepanza alla natura scalare del modello (che ignora il disordine orientazionale delle direzioni di scivolamento) e alla negligenza dell'eterogeneità elastica e degli effetti non affini nella formulazione mesoscopica.

Significatività e Rivendicazioni
L'articolo sostiene che un modello scalare coarse-grained può riprodurre le caratteristiche salienti della plasticità amorfa, a condizione che vengano incorporate specifiche informazioni strutturali dalla scala atomistica. La principale significatività risiede nell'identificazione di due ingredienti critici per l'accordo quantitativo:

1. Correlazione: La necessità di correlare i cali di stress locali con gli stress di snervamento locali, il che riflette la natura gerarchica del Paesaggio di Energia Potenziale (pozzi più profondi corrispondono a scivolamenti maggiori).
   2.Rejuvenation: La validazione dello scenario di rejuvenation, dove i riarrangiamenti plastici trasformano la struttura vetrosa locale in uno stato che somiglia a un vetro raffreddato vicino alla Transizione di Mode Coupling. Questo meccanismo è essenziale per riprodurre il comportamento di addolcimento e il banding di taglio osservati nei vetri ben ricotti.

Gli autori riconoscono con modestia che il modello è "troppo semplice e caricaturale" per un completo accordo quantitativo, specialmente per quanto riguarda le ampiezze delle fluttuazioni. Concludono che, sebbene il modello scalare catturi la fisica essenziale del flusso e della localizzazione, il lavoro futuro dovrà affrontare l'eterogeneità elastica, il disordine tensoriale e gli effetti di memoria per risolvere completamente le discrepanze nella statistica delle fluttuazioni.