Particle diffusion in large eddy structures with various Stokes number was analyzed numerically for the two-dimensional gas-particle turbulent flow over a backward-facing step. The continuous phase simulation was analyzed using the large eddy simulation (LES) method while the particle phase was analyzed by the Lagrangian method. The spatial and temporal results were used to analyze the evolution of the large eddy coherent structures and the vortex-particle interactions. The effect of the particle Stokes number on the instantaneous particle concentration distributions was also discussed. The simulations of the single-phase flow reproduced the character of the separation and reattachment flow and the essential features of the coherent structures. Numerous and complex vortex structures appeared at the high Reynolds number. Furthermore, the simulation shows that the Stokes number plays an important role in the particle dispersion and the instantaneous particle distribution. The continuous phase time-averaged results and the particle phase time-averaged results obtained from the LES analysis agree well with previous experimental data.