The diffusion of Ga and N adatoms has been studied for the technologically relevant wurtzite (0001) and (0001) surfaces employing density-functional theory. Our calculations reveal a very different diffusivity for Ga and N adatoms on the equilibrium surfaces: While Ga is very mobile at typical growth temperatures, the diffusion of N is by orders of magnitudes slower. These results give a very detailed insight of how and under which growth conditions N adatoms can be stabilized and efficiently incorporated at the surface. We further find that the presence of excess N strongly increases the Ga diffusion barrier and discuss the consequences for the growth of GaN.Recently, great progress in fabricating highly efficient GaN-based devices has been achieved [1,2]. Nevertheless, there are still substantial problems concerning growth optimization and insight into the fundamental mechanisms is rather shallow. Investigations of molecular beam epitaxy (MBE) growth have shown that the film structure and morphology are very sensitive to the III/V ratio [3,4]. GaN films grown under more N-rich conditions are rough and faceted while going towards more Ga-rich conditions a smoother surface morphology and better film quality are obtained [3].In order to improve growth in a systematic way it is essential to understand the underlying kinetic processes such as adsorption, desorption, and surface diffusion. In particular, adatom diffusion on surfaces is considered to be a key parameter controlling the growth rate, the material quality, and the surface morphology. Experimentally, an analysis of surface diffusion is difficult: So far only effective diffusion barriers have been obtained for GaN [4,5]. It is also not clear whether the cation or anion surface diffusion is the rate limiting process [4]. A problem in growing GaN is also the high N vapor pressure and considerable efforts have been made to enhance the N incorporation. However, the mechanism by which N is incorporated at the surface could not be identified. Theoretically, only few studies to describe GaN growth have been reported, which were either based on thermodynamic models [6] or Monte Carlo simulations [7]. These approaches are useful to model growth on a mesoscopic scale. However, the microscopic meaning of the effective parameters and the underlying microscopic processes remain unclear. In order to identify the fundamental growth mechanisms, but also to improve the above mentioned methods, the correct microscopic parameters are needed [8].We have therefore performed a comprehensive and detailed study of the migration and energetics of Ga and N adatoms on GaN surfaces, employing totalenergy density-functional theory calculations with ab initio pseudopotentials. We will focus on the two technologically relevant orientations for wurtzite GaN, the (0001) and the (0001) surfaces. The adatom substrate system is modeled by supercell geometries with at least 9 layers of GaN, a 14 bohr vacuum region and (2 × 2) periodicity. One side of the slab is passivated by fractional pseud...