On the limits of the energetic coupling between field dislocation mechanics and phase field crystal

Sintesi Tecnica: Limiti dell'accoppiamento energetico tra la meccanica delle dislocazioni di campo e il Phase Field Crystal

Enunciato del problema
Il documento affronta la sfida teorica e computazionale di riconciliare due distinti framework di modellazione per la plasticità mediata da dislocazioni: la Meccanica delle Dislocazioni di Campo (FDM) e il modello Phase Field Crystal (PFC).

• FDM risolve con successo il problema del valore iniziale (IBVP) per corpi continui con campi di dislocazione, catturando correttamente l'equilibrio dello stress elastico e la dinamica del trasporto. Tuttavia, manca di una cristallografia intrinseca; la localizzazione del nucleo della dislocazione deve essere prescritta tramite termini di energia non convessi dipendenti da variabili di distorsione plastica dipendenti dal percorso, il che introduce inconsistenze termodinamiche.
• PFC cattura la sottostante struttura cristallografica su scala atomica e mantiene naturalmente nuclei di dislocazione compatti attraverso un funzionale di energia libera di Swift-Hohenberg non convesso. Tuttavia, fallisce nel catturare il corretto comportamento elastico poiché le sue equazioni governanti sono derivate esclusivamente da considerazioni dissipative, portando a un rilassamento meccanico diffusivo innaturale anziché alla propagazione di onde elastiche.

Un framework accoppiato (PFCFDM), proposto da Acharya e Viñals (2020), ha tentato di unificare questi modelli introducendo un termine di penalità energetica. Tale termine penalizza la distanza $L^2$ tra la distorsione elastica ($U^e$) proveniente da FDM e la distorsione configurazionale ($Q$) derivata dal parametro d'ordine PFC. L'obiettivo è dotare il PFC di una corretta elasticità e l'FDM di una cristallografia intrinseca senza introdurre variabili dipendenti dal percorso nel funzionale di energia libera. Questo articolo investiga se tale accoppiamento raggiunga i risultati previsti.

Metodologia
Gli autori impiegano una combinazione di analisi variazionale e simulazioni numeriche per decostruire l'accoppiamento PFCFDM:

1. Analisi Variazionale: Gli autori derivano le derivate funzionali del funzionale di energia libera accoppiato rispetto al parametro d'ordine $\psi$ e alla distorsione elastica $U^e$. Analizzano la struttura dei termini di forza risultanti nell'equazione di evoluzione del phase-field e nel tensore degli sforzi modificato.
2. Decomposizione Matematica: Utilizzando la decomposizione di Stokes-Helmholtz, gli autori separano i campi di distorsione in parti compatibili (senza rotore) e incompatibili (senza divergenza). Esaminano come il termine di penalità interagisca con queste componenti, concentrandosi specificamente sull'incompatibilità che trasporta la topologia della dislocazione (vettore di Burgers).
3. Simulazioni Numeriche: Il modello accoppiato è implementato utilizzando metodi spettrali per un reticolo esagonale 2D. Le simulazioni coprono tre scenari:
   • Un dipolo di dislocazioni di bordo immobile per esaminare gli stati di equilibrio tra $U^e$ e $Q$.
   • L'evoluzione del nucleo in assenza di carichi esterni per testare la capacità di prevenire la diffusione del nucleo.
   • Scorrimento (glide) e annichilazione di dislocazioni sotto sforzo di taglio macroscopico, confrontando due paradigmi di guida: guidato da FDM (risolvendo l'equazione di trasporto per la densità di dislocazione) e guidato da PFC (impostando il tasso di distorsione plastica uguale al flusso di dislocazione PFC).

Contributi Chiave e Risultati Teorici
Il documento dimostra che l'accoppiamento energetico $L^2$ soffre di limitazioni strutturali fondamentali:

• Insensibilità alle distorsioni incompatibili: La derivata variazionale del termine di penalità rispetto al parametro d'ordine $\psi$ dipende dalla divergenza del mismatch ($\nabla \cdot (Q - U^e)$). Di conseguenza, la forza di accoppiamento agisce solo sulle parti compatibili dei campi di distorsione. Essa è cieca rispetto alle parti incompatibili (senza divergenza), che contengono tutta l'informazione relativa alla topologia della dislocazione (rotore della distorsione). Pertanto, l'accoppiamento non può imporre l'uguaglianza delle densità polari di dislocazione ($\alpha = \nabla \times U^e$ e $\tilde{\alpha} = \nabla \times Q$).
• Trasmissione Diffusiva della Meccanica: L'accoppiamento trasmette le condizioni al contorno meccaniche dal dominio FDM al dominio PFC tramite diffusione di fase piuttosto che rilassamento elastico. Il termine di penalità agisce come un operatore di diffusione sulle fasi delle ampiezze complesse. Vicino ai nuclei delle dislocazioni, dove le ampiezze svaniscono, questo coefficiente di diffusione diverge, portando a comportamenti non fisici e impedendo la trasmissione di informazioni elastiche su scale temporali rilevanti.
• Diffusione del Nucleo e Pinning: L'accoppiamento non riesce a prevenire l'innaturale diffusione dei nuclei delle dislocazioni in FDM. Inoltre, il termine di penalità agisce come una "molla di richiamo" che blocca la dislocazione PFC in posizione. Gli autori identificano una contraddizione nel parametro di accoppiamento $c_{pen}$: esso deve essere sufficientemente grande da superare il rilassamento di Swift-Hohenberg (richiedendo $c_{pen} \gg 1$) ma abbastanza piccolo da evitare il bloccaggio della dislocazione (richiedendo $c_{pen} \ll 1$). Ciò crea una finestra operativa vuota in cui il modello non funziona né come rilassatore veloce né come motore di dislocazioni mobili.
• Inconsistenza Termodinamica nell'Approccio Guidato da PFC: Mentre l'approccio guidato da FDM soffre di due densità di dislocazione distinte e irconciliabili, l'approccio guidato da PFC (imponendo $\dot{U}^p = J^\psi$) non garantisce una dissipazione plastica non negativa. I risultati numerici mostrano che la dissipazione plastica può diventare negativa, indicando che il moto indotto da PFC non è puramente guidato dallo stress e viola la coerenza termodinamica in certi regimi.

Risultati Numerici

• Equilibrio: Nelle simulazioni statiche, il sistema raggiunge un minimo locale in cui le divergenze di $U^e$ e $Q$ coincidono ($\nabla \cdot U^e = \nabla \cdot Q$), ma le parti incompatibili (e quindi i nuclei delle dislocazioni) rimangono disallineate.
• Evoluzione del Nucleo: In assenza di carico esterno, la densità FDM $\alpha$ si diffonde diffusivamente indipendentemente dalla forza di accoppiamento. La densità PFC $\tilde{\alpha}$ rimane compatta, ma le due densità evolvono indipendentemente ($\alpha \neq \tilde{\alpha}$).
• Annichilazione: Sotto sforzo applicato, il moto delle dislocazioni è severamente ostacolato. Nel caso guidato da FDM, i nuclei si diffondono e non riescono ad annichilirsi efficientemente. Nel caso guidato da PFC, il moto è o bloccato dalla "molla" di penalità o avviene a velocità inferiori rispetto alla pura evoluzione diffusiva di PFC, senza alcun recupero del rilassamento elastico.

Significatività e Assertions
Il documento conclude che l'accoppiamento energetico $L^2$ proposto non integra con successo la meccanica delle dislocazioni con la cristallografia.

• Il termine di accoppiamento agisce come una forza guidata dalla divergenza che fa corrispondere solo le parti compatibili dei campi di distorsione, lasciando l'informazione topologica (parti incompatibili) non accoppiata.
• L'informazione meccanica viene trasmessa diffusivamente anziché elasticamente, portando a un comportamento innaturale del nucleo e all'incapacità di soddisfare l'equilibrio elastico.
• Il modello soffre di un trade-off fondamentale nel parametro di accoppiamento che impedisce la simultanea localizzazione del nucleo e la mobilità della dislocazione.
• Entrambe le implementazioni, sia guidata da FDM che da PFC, presentano difetti critici: la prima fallisce nel riconciliare le densità di dislocazione, mentre la seconda viola il secondo principio della termodinamica (dissipazione plastica negativa).

Gli autori affermano che, sebbene l'accoppiamento sia termodinamicamente coerente nella sua formulazione, la sua struttura variazionale limita intrinsecamente la sua capacità di riconciliare la meccanica del continuo di FDM con la risoluzione cristallografica di PFC. Lo studio definisce i domini di applicabilità di questo specifico accoppiamento, dimostrando efficacemente che non può servire come ponte robusto tra le due teorie nella sua forma attuale.