Monte Carlo Simulations of a Generalized n--spin facilitated kinetic   Ising Model
  A kinetic Ising model is analyzed where spin variables correspond to lattice cells with mobile or immobile particles. Introducing additional restrictions for the flip processes according to the n-spin facilitated kinetic Ising model and using Monte Carlo methods we study the freezing process under the influence of an additional nearest-neighbor interaction. The stretched exponential decay of the auto-correlation function is observed and the exponent $\gamma$ as well as the relaxation time are determined depending on the activation energy $h$ and the short range coupling $J$. The magnetization corresponding to the density of immobile particles is found to be the controlling parameter for the dynamic evolution. 
B. Zheng, M. Schulz, and S. Trimper