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

Technical Summary: Limits of Energetic Coupling between Field Dislocation Mechanics and Phase Field Crystal

Problem Statement
The paper addresses the theoretical and computational challenge of reconciling two distinct modeling frameworks for dislocation-mediated plasticity: Field Dislocation Mechanics (FDM) and the Phase Field Crystal (PFC) model.

• FDM successfully solves the initial boundary value problem (IBVP) for continuum bodies with dislocation fields, capturing correct elastic stress equilibrium and transport dynamics. However, it lacks intrinsic crystallography; dislocation core localization must be prescribed via non-convex energy terms dependent on path-dependent plastic distortion variables, which introduces thermodynamic inconsistencies.
• PFC captures the underlying atomic-scale crystallographic structure and naturally maintains compact dislocation cores through a non-convex Swift-Hohenberg free energy functional. However, it fails to capture correct elastic behavior because its governing equations are derived solely from dissipative considerations, leading to unphysical diffusive mechanical relaxation rather than elastic wave propagation.

A coupled framework (PFCFDM), proposed by Acharya and Viñals (2020), attempted to unify these models by introducing an energetic penalty term. This term penalizes the $L^2$ distance between the elastic distortion ($U^e$) from FDM and the configurational distortion ($Q$) derived from the PFC order parameter. The goal was to endow PFC with correct elasticity and FDM with intrinsic crystallography without introducing path-dependent variables into the free energy. This paper investigates whether this coupling achieves its intended outcomes.

Methodology
The authors employ a combination of variational analysis and numerical simulations to deconstruct the PFCFDM coupling:

1. Variational Analysis: The authors derive the functional derivatives of the coupled free energy with respect to the order parameter $\psi$ and the elastic distortion $U^e$. They analyze the structure of the resulting forcing terms in the phase-field evolution equation and the modified stress tensor.
2. Mathematical Decomposition: Using the Stokes-Helmholtz decomposition, the authors separate distortion fields into compatible (curl-free) and incompatible (divergence-free) parts. They examine how the penalty term interacts with these components, specifically focusing on the incompatibility that carries dislocation topology (Burgers vector).
3. Numerical Simulations: The coupled model is implemented using spectral methods for a 2D hexagonal lattice. Simulations cover three scenarios:
   • An immobile edge dislocation dipole to examine equilibrium states between $U^e$ and $Q$.
   • Core evolution in the absence of external loading to test the ability to prevent core spreading.
   • Dislocation glide and annihilation under macroscopic shear stress, comparing two driving paradigms: FDM-driven (solving the transport equation for dislocation density) and PFC-driven (setting the plastic distortion rate equal to the PFC dislocation flux).

Key Contributions and Theoretical Findings
The paper demonstrates that the $L^2$ energetic coupling suffers from fundamental structural limitations:

• Insensitivity to Incompatible Distortions: The variational derivative of the penalty term with respect to the order parameter $\psi$ depends on the divergence of the mismatch ($\nabla \cdot (Q - U^e)$). Consequently, the coupling force acts only on the compatible parts of the distortion fields. It is blind to the incompatible (divergence-free) parts, which contain all information regarding dislocation topology (curl of the distortion). Therefore, the coupling cannot enforce the equality of the polar dislocation densities ($\alpha = \nabla \times U^e$ and $\tilde{\alpha} = \nabla \times Q$).
• Diffusive Transmission of Mechanics: The coupling transmits mechanical boundary conditions from the FDM domain to the PFC domain via phase diffusion rather than elastic relaxation. The penalty term acts as a diffusion operator on the phases of the complex amplitudes. Near dislocation cores, where amplitudes vanish, this diffusion coefficient diverges, leading to unphysical behavior and preventing the transmission of elastic information on relevant time scales.
• Core Spreading and Pinning: The coupling fails to prevent the unnatural spreading of dislocation cores in FDM. Furthermore, the penalty term acts as a "restoring spring" that locks the PFC dislocation in place. The authors identify a contradiction in the coupling parameter $c_{pen}$: it must be large enough to overcome the Swift-Hohenberg relaxation (requiring $c_{pen} \gg 1$) but small enough to avoid locking the dislocation (requiring $c_{pen} \ll 1$). This creates an empty operating window where the model functions neither as a fast relaxer nor as a mobile dislocation driver.
• Thermodynamic Inconsistency in PFC-Driven Approach: While the FDM-driven approach suffers from two distinct, irreconcilable dislocation densities, the PFC-driven approach (forcing $\dot{U}^p = J^\psi$) fails to guarantee non-negative plastic dissipation. Numerical results show that the plastic dissipation can become negative, indicating that the PFC-induced motion is not purely stress-driven and violates thermodynamic consistency in certain regimes.

Numerical Results

• Equilibrium: In static simulations, the system reaches a local minimum where the divergences of $U^e$ and $Q$ match ($\nabla \cdot U^e = \nabla \cdot Q$), but the incompatible parts (and thus the dislocation cores) remain mismatched.
• Core Evolution: In the absence of external load, the FDM density $\alpha$ spreads diffusively regardless of the coupling strength. The PFC density $\tilde{\alpha}$ remains compact, but the two densities evolve independently ($\alpha \neq \tilde{\alpha}$).
• Annihilation: Under applied stress, dislocation motion is severely hindered. In the FDM-driven case, cores spread and fail to annihilate efficiently. In the PFC-driven case, motion is either locked by the penalty "spring" or occurs at speeds slower than pure diffusive PFC evolution, with no recovery of elastic relaxation.

Significance and Claims
The paper concludes that the proposed $L^2$ energetic coupling does not successfully integrate dislocation mechanics with crystallography.

• The coupling term acts as a divergence-driven forcing that matches only the compatible parts of the distortion fields, leaving the topological information (incompatible parts) uncoupled.
• Mechanical information is transmitted diffusively rather than elastically, leading to unphysical core behavior and an inability to satisfy elastic equilibrium.
• The model suffers from a fundamental trade-off in the coupling parameter that prevents simultaneous core localization and dislocation mobility.
• Both FDM-driven and PFC-driven implementations exhibit critical flaws: the former fails to reconcile dislocation densities, while the latter violates the second law of thermodynamics (negative plastic dissipation).

The authors assert that while the coupling is thermodynamically consistent in its formulation, its variational structure inherently limits its ability to reconcile the continuum mechanics of FDM with the crystallographic resolution of PFC. The study defines the domains of applicability for this specific coupling, effectively showing that it cannot serve as a robust bridge between the two theories in its current form.