A statistical mechanics derivation and implementation of non-conservative phase field models for front propagation in elastic media

Over the past several decades, phase field modeling has been established as a standard simulation technique for mesoscopic science, allowing for seamless boundary tracking of moving interfaces and relatively easy coupling to other physical phenomena. However, despite its widespread success, phase field modeling remains largely driven by phenomenological justifications except in a handful of instances. In this work, we leverage a recently developed statistical mechanics framework for non-equilibrium phenomena, called Stochastic Thermodynamics with Internal Variables (STIV), to provide the first derivation of a phase field model for front propagation in a one dimensional elastic medium without appeal to phenomenology or fitting to experiments or simulation data. In the resulting model, the variables obey a gradient flow with respect to a non-equilibrium free energy, although notably, the dynamics of the strain and phase variables are coupled, and while the free energy functional is non-local in the phase field variable, it deviates from the traditional Landau-Ginzburg form. Moreover, in the systems analyzed here, the model accurately captures stress induced nucleation of transition fronts without the need to incorporate additional physics. We find that the STIV phase field model compares favorably to Langevin simulations of the microscopic system and we provide two numerical implementations enabling one to simulate arbitrary interatomic potentials.