Crystalline motion of discrete interfaces in the Blume-Emery-Griffiths model
Marco Cicalese, Giuliana Fusco, Giovanni Savaré
TL;DR
The work rigorously derives the discrete-to-continuum limit for the Blume–Emery–Griffiths surfactant model in a two-dimensional lattice under a minimizing-movements scheme with dissipation split into interface-change and surfactant-mass-variation parts. The limiting dynamics depend critically on the surfactant evaporation parameter $\\gamma$: for $\\gamma>2$ the system undergoes a crystalline mean curvature-type flow where an initial octagonal crystal shrinks while diagonals remain pinned; for $\\gamma<2$ mass conservation induces pinning and a two-stage evolution, first with pinned diagonals and then with a transition to an $\\varepsilon$-quasi-rectangle whose continuum limit is described by a coupled octagon-rectangle geometry with explicit velocity rules. The analysis hinges on preserving discrete geometric structures (octagons, staircase sets) through the minimization steps and establishing connectedness and containment properties, ultimately yielding precise continuum descriptions of the evolving interfaces. This provides a rigorous bridge between a three-phase lattice model with surfactant and its geometric continuum counterpart, clarifying how evaporation vs conservation of surfactant mass shapes interfacial motion.
Abstract
We study the discrete-to-continuum evolution of a lattice system consisting of two immiscible phases labelled by -1 and +1 in presence of a surfactant phase labelled by 0. The system's energy is described by the classical Blume-Emery-Griffith model on the lattice epsilon Z^2, and its continuum evolution is obtained as epsilon tends to zero through a minimizing-movements scheme with a time step proportional to epsilon. The dissipation functional we choose contains two contributions: a standard Almgren-Taylor-Wang type term penalizing the distance between successive configurations of the +1 phase, and a term penalizing the variation of the surfactant mass and modeling surfactant evaporation. The latter term depends on a scaling parameter gamma > 0, which determines whether the surfactant mass is conserved at each time step. We focus on the case in which the initial configuration consists of a single crystal of phase 1 completely wetted by the surfactant. For gamma > 2 the surfactant can lose mass and the evolution reduces to the crystalline mean curvature flow of an Ising-type model, while for gamma < 2 the conservation of the surfactant mass leads to a more complex evolution characterized by stronger non-uniqueness and partial pinning.
