Mechanics of incompatible asymmetric grain boundary migration

Grain boundary (GB) migration governs microstructure evolution and can mediate plastic deformation through sliding or shear coupling. Numerous experimental and numerical studies have reported a wide range of behaviors associated with boundary migration, such as defect emission or mode switching. Notably, recent studies have reported directionally asymmetric migration rates under symmetric loading, attributing this behavior to intrinsically asymmetric mobility; however, a mechanistic mesoscale explanation for this behavior remains lacking. In this work, we introduce a constitutive flow rule for grain-boundary eigendeformation within a multiphase-field framework, in which interfacial shear evolves in response to its mechanically conjugate driving force through the phase field Allen-Cahn equations. The formulation systematically employs regularized grain boundary shear kinematics informed by crystallography, and enables elastic compatibility to modulate boundary motion. Migration thresholds, residual back-stress, and apparent directional asymmetry appear naturally as emergent mechanical behavior. Simulations of symmetric and asymmetric tilt grain boundaries under mechanical, synthetic, and curvature-driven loading reveal persistent defect-like residuals following incompatible migration, transitions from planar motion to lamination at large inclinations, and even “ratcheting” behavior. These results provide a mechanically transparent explanation for behaviors such as effective mobility asymmetry and establish elastic compatibility as a constitutive mechanism in mesoscale models of boundary-mediated plasticity.

💡 Research Summary

The paper addresses a long‑standing gap in mesoscale modeling of grain‑boundary (GB) mediated plasticity: the inability of existing phase‑field (PF) frameworks to capture directionally asymmetric migration rates that have been observed experimentally under symmetric loading conditions. Traditional PF approaches prescribe a fixed shear‑coupling factor (β) for each boundary, which can reproduce shear coupling but cannot generate emergent phenomena such as migration thresholds, back‑stress, or effective mobility asymmetry without ad‑hoc parameter tuning.

To overcome this limitation, the authors embed a grain‑boundary eigendeformation field (F_gb) into a multiphase‑field description. The deformation gradient is multiplicatively decomposed as F = F_e F_gb, where F_e is the elastic part and F_gb represents the shear transformation associated with boundary motion. The eigendeformation is not a static grain‑wise constant; instead, it evolves according to a constitutive flow rule derived from the principle of minimum dissipation. The key kinematic relation is

 d F_gb_i = ½ ΔF_gb_ij ∂g_j/∂η_k dη_k,

where ΔF_gb_ij = b ⊗ h is a rank‑one tensor built from the Burgers vector b and step height vector h of the underlying disconnection, and g_j(η) are smooth indicator functions of the phase‑field order parameters η. The scalar coefficients α_n weight the contribution of each possible shear‑coupling mode, allowing the model to select among many crystallographically allowed modes.

The free energy functional combines a diffuse‑interface term (½ l_gb σ|∇η|²) with an elastic energy that is interpolated over grains:

 U(η,F) = Σ_j g_j(η) U_j(F F_gb_j⁻¹).

Taking the variational derivative of U with respect to η yields two distinct driving forces: (1) a term that pushes the boundary to reduce elastic‑modulus mismatch between neighboring grains, and (2) a term that drives the accumulation of shear to relieve incompatibility between the elastic deformation and the prescribed ΔF_gb_ij. The second term naturally produces a residual back‑stress whenever the shear‑coupling deformation is not rank‑one compatible with the interface normal, thereby generating migration thresholds and hysteresis without explicit prescription.

Numerical experiments are presented in three families. First, symmetric and asymmetric tilt boundaries are loaded mechanically. Because the flow rule allows α_n to differ for opposite senses of motion, the same symmetric load produces different migration velocities in opposite directions, reproducing the “effective mobility asymmetry” reported in recent experiments. Second, curvature‑driven migration of highly inclined boundaries shows a transition from planar motion to lamination (layered structures) when the shear coupling cannot be accommodated compatibly; the lamination is a direct manifestation of the residual shear stored in the bulk. Third, a cyclic loading protocol reveals a “ratcheting” effect: residual defect‑like shear fields left behind after an incompatible migration episode act as nucleation sites for subsequent motion, leading to accumulated plastic strain even though each individual loading step is symmetric.

The model’s salient features are: (i) a physically based, crystallography‑informed shear‑coupling tensor rather than an ad‑hoc scalar; (ii) a history‑dependent eigendeformation field that captures the accumulation of shear and back‑stress; (iii) emergent migration thresholds and directional asymmetry arising naturally from elastic compatibility constraints; and (iv) a straightforward implementation within existing PF codes because the driving force is obtained analytically from the energy functional.

Limitations are acknowledged: anisotropic grain‑boundary energy, higher‑order gradient regularization, and explicit coupling to bulk plasticity (F_p) are not included, and the current formulation assumes a single dominant disconnection mode with a common normal. Nonetheless, the framework provides a transparent mechanical explanation for a suite of complex GB behaviors—mobility asymmetry, lamination, defect‑like residuals, and ratcheting—that have previously required phenomenological modeling.

In conclusion, by integrating a constitutive flow rule for GB eigendeformation into a multiphase‑field setting, the authors deliver a mesoscale model that captures both the thermodynamic and mechanical drivers of GB migration. The approach bridges the gap between atomistic insight (disconnections, shear coupling) and continuum plasticity, offering a powerful tool for predicting microstructure evolution and boundary‑mediated plasticity in engineering alloys.