Long time evolution of meandered steps during the crystal growth   processes
  Step meandering during growth of gallium nitride (0001) surface is studied using kinetic Monte Carlo method. Simulated growth process, conducted in N-rich conditions are therefore controlled by Ga atoms surface diffusion. The model employs dominating four-body interactions of Ga atoms that cause step flow anisotropy during growth. Overall kinetics and shape selection features of step meandering are analyzed assuming their dependence on the external particle flux and on the temperature. It appears that at relatively high temperatures and low fluxes steps move regularly preserving their initial shapes of straight, parallel lines. For higher fluxes and at wide range of temperatures step meandering happens. It is shown that, depending on the initial surface parameters, two different scenarios of step meandering are realized. In both these regimes meandering has different character as a function of time. For relatively high fluxes meanders have wavelengths shorter than the terrace width and they grow independently. Eventually surface ends up as a rough structure. When flux is lower regular pattern of meanders emerges. Step meander even if Schwoebel barrier is absent. Too high barrier destroys step stability. The amplitude of wavelike step meanders increases in time up to a saturation value. The mechanism of such meander development is elucidated. 
Magdalena A. Za{\l}uska-Kotur, Filip Krzy\.zewski and Stanis{\l}aw   Krukowski
cond-mat.stat-mech cond-mat.mtrl-sci