Abstract: This talk is concerned with the numerical simulations for the dynamics of the molecular beam epitaxy (MBE) models. The numerical simulations of the MBE models require long time computations, and therefore large time-stepping methods become necessary. In this talk, we will discuss on some recent developments on numerical simulations of two MBE models.