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

Resumo Técnico: Limites do Acoplamento Energético entre a Mecânica de Dislocações de Campo e o Cristal de Campo de Fase

Problema
O artigo aborda o desafio teórico e computacional de reconciliar dois frameworks de modelagem distintos para a plasticidade mediada por dislocações: a Mecânica de Dislocações de Campo (FDM) e o modelo de Cristal de Campo de Fase (PFC).

• FDM resolve com sucesso o problema de valor inicial (IBVP) para corpos contínuos com campos de dislocação, capturando corretamente o equilíbrio de tensão elástica e a dinâmica de transporte. No entanto, carece de uma cristalografia intrínseca; a localização do núcleo da dislocação deve ser prescrita via termos de energia não convexos dependentes de variáveis de distorção plástica dependentes do caminho, o que introduz inconsistências termodinâmicas.
• PFC captura a estrutura cristalográfica subjacente em escala atômica e mantém naturalmente núcleos de dislocação compactos através de um funcional de energia livre de Swift-Hohenberg não convexo. Contudo, falha em capturar o comportamento elástico correto porque suas equações governantes são derivadas puramente de considerações dissipativas, levando a uma relaxação mecânica difusiva não física em vez de propagação de ondas elásticas.

Um framework acoplado (PFCFDM), proposto por Acharya e Viñals (2020), tentou unificar esses modelos introduzindo um termo de penalidade energética. Este termo penaliza a distância $L^2$ entre a distorção elástica ($U^e$) da FDM e a distorção configuracional ($Q$) derivada do parâmetro de ordem PFC. O objetivo é dotar o PFC de elasticidade correta e a FDM de cristalografia intrínseca sem introduzir variáveis dependentes do caminho na energia livre. Este artigo investiga se este acoplamento alcança os resultados pretendidos.

Metodologia
Os autores utilizam uma combinação de análise variacional e simulações numéricas para desconstruir o acoplamento PFCFDM:

1. Análise Variacional: Os autores derivam as derivadas funcionais da energia livre acoplada em relação ao parâmetro de ordem $\psi$ e à distorção elástica $U^e$. Eles analisam a estrutura dos termos de força resultantes na equação de evolução do campo de fase e no tensor de tensão modificado.
2. Decomposição Matemática: Utilizando a decomposição de Stokes-Helmholtz, os autores separam os campos de distorção em partes compatíveis (livres de rotacional/curl-free) e incompatíveis (livres de divergência/divergence-free). Eles examinam como o termo de penalidade interage com essas componentes, focando especificamente na incompatibilidade que carrega a topologia da dislocação (vetor de Burgers).
3. Simulações Numéricas: O modelo acoplado é implementado usando métodos espectrais para uma rede hexagonal 2D. As simulações cobrem três cenários:
   • Um dipolo de dislocação de aresta imóvel para examinar estados de equilíbrio entre $U^e$ e $Q$.
   • Evolução do núcleo na ausência de carga externa para testar a capacidade de prevenir o espalhamento do núcleo.
   • Deslizamento e aniquilação de dislocação sob tensão de cisalhamento macroscópica, comparando dois paradigmas de condução: conduzido por FDM (resolvendo a equação de transporte para a densidade de dislocação) e conduzido por PFC (definindo a taxa de distorção plástica igual ao fluxo de dislocação PFC).

Contribuições Principais e Achados Teóricos
O artigo demonstra que o acoplamento energético $L^2$ sofre de limitações estruturais fundamentais:

• Insensibilidade a Distorções Incompatíveis: A derivada variacional do termo de penalidade em relação ao parâmetro de ordem $\psi$ depende da divergência do desajuste ($\nabla \cdot (Q - U^e)$). Consequentemente, a força de acoplamento atua apenas nas partes compatíveis dos campos de distorção. Ela é cega para as partes incompatíveis (livres de divergência), que contêm toda a informação relativa à topologia da dislocação (rotacional da distorção). Portanto, o acoplamento não consegue impor a igualdade das densidades polares de dislocação ($\alpha = \nabla \times U^e$ e $\tilde{\alpha} = \nabla \times Q$).
• Transmissão Difusiva de Mecânica: O acoplamento transmite condições de contorno mecânicas do domínio FDM para o domínio PFC via difusão de fase, em vez de relaxação elástica. O termo de penalidade atua como um operador de difusão nas fases das amplitudes complexas. Próximo aos núcleos de dislocação, onde as amplitudes desaparecem, esse coeficiente de difusão diverge, levando a comportamentos não físicos e impedindo a transmissão de informação elástica em escalas de tempo relevantes.
• Espalhamento e Travamento do Núcleo: O acoplamento falha em prevenir o espalhamento não natural dos núcleos de dislocação na FDM. Além disso, o termo de penalidade atua como uma "mola de restauração" que trava a dislocação do PFC no lugar. Os autores identificam uma contradição no parâmetro de acoplamento $c_{pen}$: ele deve ser grande o suficiente para superar a relaxação de Swift-Hohenberg (exigindo $c_{pen} \gg 1$), mas pequeno o suficiente para evitar o travamento da dislocação (exigindo $c_{pen} \ll 1$). Isso cria uma janela operacional vazia onde o modelo não funciona nem como um relaxador rápido, nem como um motor de mobilidade de dislocação.
• Inconsistência Termodinâmica na Abordagem Conduzida por PFC: Enquanto a abordagem conduzida por FDM sofre com duas densidades de dislocação distintas e irreconciliáveis, a abordagem conduzida por PFC (forçando $\dot{U}^p = J^\psi$) falha em garantir a não negatividade da dissipação plástica. Resultados numéricos mostram que a dissipação plástica pode tornar-se negativa, indicando que o movimento induzido pelo PFC não é puramente movido por tensão e viola a consistência termodinâmica em certos regimes.

Resultados Numéricos

• Equilíbrio: Em simulações estáticas, o sistema atinge um mínimo local onde as divergências de $U^e$ e $Q$ coincidem ($\nabla \cdot U^e = \nabla \cdot Q$), mas as partes incompatíveis (e, portanto, os núcleos das dislocações) permanecem desalinhadas.
• Evolução do Núcleo: Na ausência de carga externa, a densidade FDM $\alpha$ se espalha diffusivamente, independentemente da força de acoplamento. A densidade PFC $\tilde{\alpha}$ permanece compacta, mas as duas densidades evoluem independentemente ($\alpha \neq \tilde{\alpha}$).
• Aniquilação: Sob tensão aplicada, o movimento da dislocação é severamente dificultado. No caso conduzido por FDM, os núcleos se espalham e falham em aniquilar eficientemente. No caso conduzido por PFC, o movimento é ou travado pela "mola" de penalidade ou ocorre em velocidades menores do que a evolução puramente difusiva do PFC, sem recuperação da relaxação elástica.

Significância e Alegações
O artigo conclui que o proposto acoplamento energético $L^2$ não integra com sucesso a mecânica de dislocação com a cristalografia.

• O termo de acoplamento atua como uma força dirigida pela divergência que combina apenas as partes compatíveis dos campos de distorção, deixando a informação topológica (partes incompatíveis) sem acoplamento.
• A informação mecânica é transmitida diffusivamente em vez de elasticamente, levando a comportamentos de núcleo não físicos e à incapacidade de satisfazer o equilíbrio elástico.
• O modelo sofre de um conflito fundamental no parâmetro de acoplamento que impede a simultânea localização do núcleo e a mobilidade da dislocação.
• Ambas as implementações, conduzidas por FDM e por PFC, exibem falhas críticas: a primeira falha em reconciliar as densidades de dislocação, enquanto a segunda viola a segunda lei da termodinâmica (dissipação plástica negativa).

Os autores afirmam que, embora o acoplamento seja termodinamicamente consistente em sua formulação, sua estrutura variacional limita inerentemente sua capacidade de reconciliar a mecânica contínua da FDM com a resolução cristalográfica do PFC. O estudo define os domínios de aplicabilidade para este acoplamento específico, mostrando efetivamente que ele não pode servir como uma ponte robusta entre as duas teorias em sua forma atual.