“…The two coefficients A and B are determined by the boundary conditions at the steps: if, for the sake of simplicity, we assume no Ehrlich-Schwöbel effect [31], the symmetry imposes B = 0 if the position is measured from the center of the terrace. The parameter A is determined by the fact that the current related to the gradient of adatoms balances the incorporation of adatoms on both sides at the step: ∓D∇n(±l s ) = β[n(±l s ) − n eq ], where the parameter β sets the efficiency of trapping at the step [22,35], and n eq is determined by the sublimation rate V sub of the material, n eq = V sub τ . The sublimation process involves the release of an atom already incorporated, generally at a step edge, adding thus an adatom to the neighboring terrace [33], which then diffuses and may be desorbed or re-incorporated.…”