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

Résumé Technique : Limites du couplage énergétique entre la mécanique des dislocations de champ et le cristal à champ de phase

Énoncé du problème
Le document traite du défi théorique et computationnel consistant à réconcilier deux cadres de modélisation distincts pour la plasticité médiée par les dislocations : la mécanique des dislocations de champ (FDM - Field Dislocation Mechanics) et le modèle de cristal à champ de phase (PFC - Phase Field Crystal).

• La FDM résout avec succès le problème de Cauchy (IBVP) pour les corps continus avec des champs de dislocations, capturant l'équilibre des contraintes élastiques et la dynamique de transport. Cependant, elle manque de cristallographie intrinsèque ; la localisation du cœur de dislocation doit être prescrite via des termes d'énergie non convexes dépendant de variables de distorsion plastique dépendantes du chemin, ce qui introduit des incohérences thermodynamiques.
• Le PFC capture la structure cristallographique sous-jacente à l'échelle atomique et maintient naturellement des cœurs de dislocations compacts grâce à une fonctionnelle d'énergie libre de Swift-Hohenberg non convexe. Cependant, il échoue à capturer un comportement élastique correct car ses équations directrices sont dérivées uniquement de considérations dissipatives, menant à une relaxation mécanique diffusive non physique plutôt qu'à une propagation d'ondes élastiques.

Un cadre couplé (PFCFDM), proposé par Acharya et Viñals (2020), a tenté d'unifier ces modèles en introduisant un terme de pénalité énergétique. Ce terme pénalise la distance $L^2$ entre la distorsion élastique ($U^e$) de la FDM et la distorsion configurationnelle ($Q$) dérivée du paramètre d'ordre du PFC. L'objectif est de doter le PFC d'une élasticité correcte et la FDM d'une cristallographie intrinsèque sans introduire de variables dépendantes du chemin dans l'énergie libre. Cet article examine si ce couplage atteint les résultats escomptés.

Méthodologie
Les auteurs utilisent une combinaison d'analyse variationnelle et de simulations numériques pour déconstruire le couplage PFCFDM :

1. Analyse variationnelle : Les auteurs dérivent les dérivées fonctionnelles de l'énergie libre couplée par rapport au paramètre d'ordre $\psi$ et à la distorsion élastique $U^e$. Ils analysent la structure des termes de forçage résultants dans l'équation d'évolution du champ de phase et le tenseur de contrainte modifié.
2. Décomposition mathématique : En utilisant la décomposition de Stokes-Helmholtz, les auteurs séparent les champs de distorsion en parties compatibles (sans rotation) et incompatibles (sans divergence). Ils examinent comment le terme de pénalité interagit avec ces composantes, en se concentrant spécifiquement sur l'incompatibilité qui porte la topologie de la dislocation (vecteur de Burgers).
3. Simulations numériques : Le modèle couplé est implémenté à l'aide de méthodes spectrales pour un réseau hexagonal 2D. Les simulations couvrent trois scénarios :
   • Un dipôle de dislocation de bord immobile pour examiner les états d'équilibre entre $U^e$ et $Q$.
   • L'évolution du cœur en l'absence de charge externe pour tester la capacité à empêcher l'étalement du cœur.
   • Le glissement et l'annihilation de dislocations sous une contrainte de cisaillement macroscopique, en comparant deux paradigmes de pilotage : piloté par la FDM (résolution de l'équation de transport de la densité de dislocations) et piloté par le PFC (définition du taux de distorsion plastique égal au flux de dislocation du PFC).

Contributions clés et résultats théoriques
L'article démontre que le couplage énergétique $L^2$ souffre de limitations structurelles fondamentales :

• Insensibilité aux distorsions incompatibles : La dérivée variationnelle du terme de pénalité par rapport au paramètre d'ordre $\psi$ dépend de la divergence du désaccord ($\nabla \cdot (Q - U^e)$). Par conséquent, la force de couplage n'agit que sur les parties compatibles des champs de distorsion. Elle est aveugle aux parties incompatibles (sans divergence), qui contiennent toute l'information concernant la topologie de la dislocation (rotation de la distorsion). Par conséquent, le couplage ne peut pas imposer l'égalité des densités de dislocations polaires ($\alpha = \nabla \times U^e$ et $\tilde{\alpha} = \nabla \times Q$).
• Transmission diffusive de la mécanique : Le couplage transmet les conditions limites mécaniques du domaine FDM vers le domaine PFC via la diffusion de phase plutôt que par relaxation élastique. Le terme de pénalité agit comme un opérateur de diffusion sur les phases des amplitudes complexes. Près des cœurs de dislocations, où les amplitudes s'annulent, ce coefficient de diffusion diverge, entraînant un comportement non physique et empêchant la transmission de l'information élastique sur les échelles de temps pertinentes.
• Étalement et ancrage du cœur : Le couplage ne parvient pas à empêcher l'étalement non naturel des cœurs de dislocations dans la FDM. De plus, le terme de pénalité agit comme un « ressort de rappel » qui verrouille la dislocation du PFC sur place. Les auteurs identifient une contradiction dans le paramètre de couplage $c_{pen}$ : il doit être suffisamment grand pour surmonter la relaxation de Swift-Hohenberg (nécessitant $c_{pen} \gg 1$) mais suffisamment petit pour éviter de verrouiller la dislocation (nécessitant $c_{pen} \ll 1$). Cela crée une fenêtre opérationnelle vide où le modèle ne fonctionne ni comme un relaxateur rapide, ni comme un moteur de dislocation mobile.
• Incohérence thermodynamique dans l'approche pilotée par le PFC : Alors que l'approche pilotée par la FDM souffre de deux densités de dislocations distinctes et irréconciliables, l'approche pilotée par le PFC (imposant $\dot{U}^p = J^\psi$) ne garantit pas une dissipation plastique non négative. Les résultats numériques montrent que la dissipation plastique peut devenir négative, indiquant que le mouvement induit par le PFC n'est pas purement piloté par la contrainte et viole la cohérence thermodynamique dans certains régimes.

Résultats numériques

• Équilibre : Dans les simulations statiques, le système atteint un minimum local où les divergences de $U^e$ et $Q$ correspondent ($\nabla \cdot U^e = \nabla \cdot Q$), mais les parties incompatibles (et donc les cœurs de dislocations) restent désalignées.
• Évolution du cœur : En l'absence de charge externe, la densité FDM $\alpha$ s'étale de manière diffusive, quel que soit la force du couplage. La densité PFC $\tilde{\alpha}$ reste compacte, mais les deux densités évoluent indépendamment ($\alpha \neq \tilde{\alpha}$).
• Annihilation : Sous contrainte appliquée, le mouvement des dislocations est sévèrement entravé. Dans le cas piloté par la FDM, les cœurs s'étalent et ne parviennent pas à s'annihiler efficacement. Dans le cas piloté par le PFC, le mouvement est soit verrouillé par le « ressort » de pénalité, soit se produit à des vitesses plus lentes que l'évolution diffusive pure du PFC, sans récupération de la relaxation élastique.

Signification et affirmations
L'article conclut que le couplage énergétique $L^2$ proposé ne parvient pas à intégrer la mécanique des dislocations et la cristallographie.

• Le terme de couplage agit comme un forçage piloté par la divergence qui ne fait correspondre que les parties compatibles des champs de distorsion, laissant l'information topologique (parties incompatibles) non couplée.
• L'information mécanique est transmise diffusivement plutôt qu'élastiquement, conduisant à un comportement de cœur non physique et à une incapacité à satisfaire l'équilibre élastique.
• Le modèle souffre d'un compromis fondamental dans le paramètre de couplage qui empêche la localisation simultanée du cœur et la mobilité de la dislocation.
• Les deux implémentations, pilotée par la FDM ou par le PFC, présentent des défauts critiques : la première échoue à réconcilier les densités de dislocations, tandis que la seconde viole le second principe de la thermodynamique (dissipation plastique négative).

Les auteurs affirment que bien que le couplage soit thermodynamiquement cohérent dans sa formulation, sa structure variationnelle limite intrinsèquement sa capacité à réconcilier la mécanique des milieux continus de la FDM avec la résolution cristallographique du PFC. L'étude définit les domaines d'application de ce couplage spécifique, montrant efficacement qu'il ne peut servir de pont robuste entre les deux théories sous sa forme actuelle.