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

Resumen Técnico: Límites del Acoplamiento Energético entre la Mecánica de Dislocaciones de Campo y el Cristal de Campo de Fase

Planteamiento del Problema
El artículo aborda el desafío teórico y computacional de reconciliar dos marcos de modelado distintos para la plasticidad mediada por dislocaciones: la Mecánica de Dislocaciones de Campo (FDM, por sus siglas en inglés) y el modelo de Cristal de Campo de Fase (PFC).

• FDM resuelve con éxito el problema de valor inicial y de contorno (IBVP) para cuerpos continuos con campos de dislocación, capturando correctamente el equilibrio de tensiones elásticas y la dinámica de transporte. Sin embargo, carece de una cristaloografía intrínseca; la localización del núcleo de la dislocación debe prescribirse mediante términos de energía no convexos que dependen de variables de distorsión plástica dependientes de la trayectoria, lo que introduce inconsistencias termodinámicas.
• PFC captura la estructura cristalográfica subyacente a escala atómica y mantiene naturalmente núcleos de dislocación compactos mediante un funcional de energía libre de Swift-Hohenberg no convexo. No obstante, falla al capturar un comportamiento elástico correcto porque sus ecuaciones gobernantes se derivan únicamente de consideraciones disipativas, lo que conduce a una relajación mecánica difusiva no física en lugar de una propagación de ondas elásticas.

Un marco acoplado (PFCFDM), propuesto por Acharya y Viñals (2020), intentó unificar estos modelos introduciendo un término de penalización energética. Este término penaliza la distancia $L^2$ entre la distorsión elástica ($U^e$) de FDM y la distorsión configuracional ($Q$) derivada del parámetro de orden de PFC. El objetivo era dotar a PFC de una elasticidad correcta y a FDM de una cristaloografía intrínseca sin introducir variables dependientes de la trayectoria en la energía libre. Este artículo investiga si este acoplamiento logra los resultados previstos.

Metodología
Los autores emplean una combinación de análisis variacional y simulaciones numéricas para deconstruir el acoplamiento PFCFDM:

1. Análisis Variacional: Los autores derivan las derivadas funcionales de la energía libre acoplada con respecto al parámetro de orden $\psi$ y la distorsión elástica $U^e$. Analizan la estructura de los términos de forzamiento resultantes en la ecuación de evolución del campo de fase y el tensor de tensiones modificado.
2. Descomposición Matemática: Utilizando la descomposición de Stokes-Helmholtz, los autores separan los campos de distorsión en partes compatibles (libres de rotacional/curl-free) e incompatibles (libres de divergencia/divergence-free). Examinan cómo el término de penalización interactúa con estos componentes, centrándose específicamente en la incompatibilidad que transporta la topología de la dislocación (vector de Burgers).
3. Simulaciones Numéricas: El modelo acoplado se implementa utilizando métodos espectrales para una red hexagonal 2D. Las simulaciones cubren tres escenarios:
   • Un dipolo de dislocación de borde inmóvil para examinar los estados de equilibrio entre $U^e$ y $Q$.
   • La evolución del núcleo en ausencia de carga externa para probar la capacidad de prevenir la dispersión del núcleo.
   • Deslizamiento y aniquilación de dislocaciones bajo tensión de cizalla macroscópica, comparando dos paradigmas de conducción: impulsado por FDM (resolviendo la ecuación de transporte para la densidad de dislocaciones) e impulsado por PFC (estableciendo la tasa de distorsión plástica igual al flujo de dislocación de PFC).

Contribuciones Clave y Hallazgos Teóricos
El artículo demuestra que el acoplamiento energético $L^2$ sufre limitaciones estructurales fundamentales:

• Insensibilidad a las Distorsiones Incompatibles: La derivada variacional del término de penalización con respecto al parámetro de orden $\psi$ depende de la divergencia del desajuste ($\nabla \cdot (Q - U^e)$). En consecuencia, la fuerza de acoplamiento actúa solo sobre las partes compatibles de los campos de distorsión. Es ciega a las partes incompatibles (libres de divergencia), que contienen toda la información relativa a la topología de la dislocación (rotacional de la distorsión). Por lo tanto, el acoplamiento no puede imponer la igualdad de las densidades polares de dislocación ($\alpha = \nabla \times U^e$ y $\tilde{\alpha} = \nabla \times Q$).
• Transmisión Mecánica Difusiva: El acoplamiento transmite las condiciones de contorno mecánicas del dominio FDM al dominio PFC mediante difusión de fase en lugar de relajación elástica. El término de penalización actúa como un operador de difusión sobre las fases de las amplitudes complejas. Cerca de los núcleos de dislocación, donde las amplitudes se anulan, este coeficiente de difusión diverge, lo que provoca un comportamiento no físico y evita la transmisión de información elástica en escalas de tiempo relevantes.
• Dispersión y Bloqueo del Núcleo: El acoplamiento no logra prevenir la dispersión antinatural de los núcleos de dislocación en FDM. Además, el término de penalización actúa como un "resorte de restauración" que bloquea la dislocación de PFC en su lugar. Los autores identifican una contradicción en el parámetro de acoplamiento $c_{pen}$: debe ser lo suficientemente grande como para superar la relajación de Swift-Hohenberg (requiriendo $c_{pen} \gg 1$) pero lo suficientemente pequeño como para evitar el bloqueo de la dislocación (requiriendo $c_{pen} \ll 1$). Esto crea una ventana operativa vacía donde el modelo no funciona ni como un relajador rápido ni como un motor de dislocación móvil.
• Inconsistencia Termodinámica en el Enfoque Impulsado por PFC: Mientras que el enfoque impulsado por FDM sufre de dos densidades de dislocación distintas e irreconciliables, el enfoque impulsado por PFC (forzando $\dot{U}^p = J^\psi$) no garantiza una disipación plástica no negativa. Los resultados numéricos muestran que la disipación plástica puede volverse negativa, lo que indica que el movimiento inducido por PFC no es puramente impulsado por tensiones y viola la consistencia termodinámica en ciertos regímenes.

Resultados Numéricos

• Equilibrio: En simulaciones estáticas, el sistema alcanza un mínimo local donde las divergencias de $U^e$ y $Q$ coinciden ($\nabla \cdot U^e = \nabla \cdot Q$), pero las partes incompatibles (y, por ende, los núcleos de dislocación) permanecen desajustadas.
• Evolución del Núcleo: En ausencia de carga externa, la densidad de FDM $\alpha$ se dispersa difusivamente independientemente de la fuerza de acoplamiento. La densidad de PFC $\tilde{\alpha}$ permanece compacta, pero las dos densidades evolucionan de forma independiente ($\alpha \neq \tilde{\alpha}$).
• Aniquilación: Bajo tensión aplicada, el movimiento de las dislocaciones se ve severamente obstaculizado. En el caso impulsado por FDM, los núcleos se dispersan y no logran aniquilarse eficientemente. En el caso impulsado por PFC, el movimiento o bien es bloqueado por el "resorte" de la penalización, o bien ocurre a velocidades menores que la evolución puramente difusiva de PFC, sin recuperación de la relajación elástica.

Significancia y Reivindicaciones
El artículo concluye que el propuesto acoplamiento energético $L^2$ no integra con éxito la mecánica de dislocaciones con la cristaloografía.

• El término de acoplamiento actúa como un forzamiento impulsado por la divergencia que iguala solo las partes compatibles de los campos de distorsión, dejando la información topológica (partes incompatibles) sin acoplar.
• La información mecánica se transmite difusivamente en lugar de elásticamente, lo que conduce a un comportamiento de núcleo no físico y a la incapacidad de satisfacer el equilibrio elástico.
• El modelo sufre de un compromiso fundamental en el parámetro de acoplamiento que impide la localización simultánea del núcleo y la movilidad de la dislocación.
• Tanto las implementaciones impulsadas por FDM como las impulsadas por PFC exhiben fallos críticos: la primera no logra reconciliar las densidades de dislocación, mientras que la segunda viola la segunda ley de la termodinámica (disipación plástica negativa).

Los autores sostienen que, aunque el acoplamiento es termodinámicamente consistente en su formulación, su estructura variacional limita inherentemente su capacidad para reconciliar la mecánica de medios continuos de FDM con la resolución cristalográfica de PFC. El estudio define los dominios de aplicabilidad para este acoplamiento específico, demostrando efectivamente que no puede servir como un puente robusto entre ambas teorías en su forma actual.