We present a dynamical simulation approach by Monte Carlo calculation to describe atomic mixing phenomena in depth profiling of multilayered systems. This approach is based on the binary encounter model, taking into account (1) generation of both the interstitial atoms and vacancies and (2) annihilation of the vacancies. The results indicate that the simulation describes very well the depth profiles of AlAs/GaAs multilayered systems obtained by Auger electron spectroscopy. It predicts the existence of preferential sputtering of Al in the AlAs layer, where the Auger intensity of the As (MVV-32 eV) in the AlAs layer is about 1.2 times larger than that of the GaAs layer for 0.5 keV Ar+ sputtering.