Surface diffusion was studied in a two barrier system, i.e. a system with two different activation energies. The dynamic and thermodynamic properties of the system of interacting ad-atoms on an fcc(111) surface were studied using Monte Carlo simulations, the Bortz-Kalos-Lebowitz algorithm and the Metropolis algorithm. The results were compared with previous analytical calculations in the same system. Numerical results fitted well the analytical results for high temperatures. At low temperatures, contrary to previous analytical calculation (based on quasi-chemical approximation), a sharp local minimum was observed at a coverage Θ = 0.5 in the case with repulsive interaction and a large difference between activation energies. The Monte Carlo simulations result in a higher occupation of deep sites, formation of an ordered phase and in lower diffusion coefficient in comparison with the analytical results.

Monte-Carlo Simulation of Diffusion in a Two Barriers System. M.Mašín, Z.Chvoj, P.Jelínek, J.Heinonen: Surface Science, 2002, 516[1-2], 169-78