Collective surface diffusion: n -fold way kinetic Monte Carlo simulation

Collective surface diffusion of strongly interacting particles is simulated on the basis of an n-fold way kinetic Monte Carlo scheme. The coverage dependence of the jump and tracer diffusion coefficients is calculated for one- and two-dimensional lattice gases at very low (subcritical) temperatures. Results are compared with exact analytical ones and Monte Carlo simulations using the standard Metropolis algorithm. The method proves to be highly reliable to investigate surface diffusion at subcritical temperatures where phase coexistence in the adlayer occurs.