Nonconserved critical dynamics of the two-dimensional Ising model as a surface kinetic roughening process

We have revisited the nonconserved (or model A) critical dynamics of the two-dimensional Ising model through numerical simulations of its lattice and continuum formulations—Glauber dynamics and the time-dependent Ginzburg-Landau (TDGL) equation, respectively—to analyze them with current tools from surface kinetic roughening. Our study examines two critical quenches, one from an ordered and other one from a disordered initial state, for both of which we assess the dynamic scaling ansatz, the critical exponent values, and the fluctuation field statistics that occur. Notably, the dynamic scaling ansatz followed by the system strongly depends on the initial condition: A critical quench from the ordered phase follows Family-Vicsek (FV) scaling, while a critical quench from the disordered phase shows an initial overgrowth regime with intrinsic anomalous roughening, followed by relaxation to equilibrium. This behavior is explained in terms of the dynamical instability of the stochastic Ginzburg-Landau equation at the critical temperature, whereby the linearly unstable term is eventually stabilized by nonlinear interactions. For both quenches, we have determined the occurrence of probability distribution functions for the field fluctuations, which are time independent along the nonequilibrium dynamics when suitably normalized by the time-dependent fluctuation amplitude. Additionally, we have developed a related interface model for a field that scales as the space integral of the TDGL field (integral GL model). Numerical simulations of this model reveal either FV or faceted anomalous roughening, depending on the critical quenched performed, as well as an emergent symmetry in the fluctuation statistics for a critical quench from the ordered phase.