“…Phase field models of spinodal decomposition and stress-induced MT are formulated based on a combination of Landau theory [20], gradient thermodynamics [21], and the modified Khachaturyan-Shatalov's microelasticity theory [22,23]. The concentration evolution during spinodal decomposition is governed by the Cahn-Hilliard equation [24], and the temporal and spatial evolution of the structural order parameters (which represent the parent and martensitic phases) during MT are described by the time-dependent Ginzburg-Landau equation [25][26][27][28]. The symmetry breaking accompanying the β (BCC, point group 𝑚3𝑚) → α (orthorhombic, point group mmm) MT [29,30] (see Fig.…”