Both　the　formation　energies　and　the　intra-　and　inter-layer　diffuse　activation　energies　of　a　vacancy　in　the　first　six　lattice　planes　of　Mg　(001)　surface　have　been　calculated　by　combining　the　modified　analytical　embedded-atom　method　(MAEAM)　with　molecular　dynamics　(MID).　The　results　show　that　the　effect　of　the　surface　on　the　formation　and　migration　of　the　vacancy　is　only　down　to　the　third-layer.　It　is　easer　for　a　single　vacancy　to　form　and　to　migrate　in　the　first　layer.　Furthermore,　the　vacancy　in　the　second　layer　is　favorable　to　migrate　to　the　first　layer.　This　is　in　agreement　with　the　experimental　results　that　the　first　layer　has　the　highest　concentration　of　the　vacancy.　(C)　2008　WILEY-VCH　GmbH　&　Co.　KGaA,　Weinbeim.