---
abstract: 'The role played by a kinetic barrier originated by out-of-plane step edge diffusion, introduced in \[Leal *et al.*, [J. Phys. Condens. Matter **23**, 292201 (2011)](https://doi.org/10.1088/0953-8984/23/29/292201)\], is investigated in the Wolf-Villain and Das Sarma-Tamborenea models with short range diffusion. Using large-scale simulations, we observe that this barrier is sufficient to produce growth instability, forming quasiregular mounds in one and two dimensions. The characteristic surface length saturates quickly indicating a uncorrelated growth of the three-dimensional structures, which is also confirmed by a growth exponent $\beta=1/2$. The out-of-plane particle current shows a large reduction of the downward flux in the presence of the kinetic barrier enhancing, consequently, the net upward diffusion and the formation of three-dimensional self-assembled structures.'
address: 'Departamento de Estatística, Física e Matemática, Campus Alto Paraopeba, Universidade Federal de São João Del-Rei, 36420-000, Ouro Branco, MG, Brazil'
author:
- 'Anderson J. Pereira'
- 'Sidiney G. Alves'
- 'Silvio C. Ferreira'
title: 'Effects of a kinetic barrier on limited-mobility interface growth models'
---

=1

Introduction
============

A rich variety of morphologies can be observed during far-from-equilibrium growth processes and many of them with potential for technological applications [@michely2004islands; @Evans2006; @barabasi; @meakin]. Growth instability can induce three-dimensional mound-like patterns in different types of films such as metals [@Jorritsma; @Caspersen; @Han], inorganic [@Johnson; @Tadayyon] and organic [@Zorba; @Hlawacek] semiconductors materials to cite only a few examples. Such a growth instability has been mainly attributed to the presence of Ehrlich-Schwoebel (ES) step barriers [@Ehrlich; @Schwoebel] that reduce the rate with which atoms move downwardly on the edges of terraces leading to net uphill flows. Growth instabilities can also emerge from topologically induced uphill currents which depend on the crystalline structure [@Kanjanaput2010] or from fast diffusion on terrace edges [@Murty2003; @Pierre-Louis1999] among other mechanisms [@Evans2006; @michely2004islands]. The existence of ES barriers is supported by molecular dynamic simulations [@Yang].

![image](film_wv_1d.pdf){width="0.35\linewidth"}     ![image](film_dt_1d.pdf){width="0.35\linewidth"}

Discrete solid-on-solid (SOS) growth models constitute an important approach to investigate the dynamic of kinetic roughening and morphological properties of interfaces. The rules are easily implemented in a discrete space (lattices) rid of overhangs and bulk voids. The role played by ES barriers has been investigated in models with thermally activated diffusion [@Evans2006; @michely2004islands] being the Clark-Vvedenski (CV) model [@CV_PRL; @Clarke1988] one of the simplest examples, in which any surface adatom can move according to an Arrhenius diffusion coefficient $D\sim
\exp(-E/k_B T)$ [@barabasi] where $E$ is an energy activation barrier to be overcome in a diffusion hopping. An ES barrier can be included as an additional activation energy for diffusion at the edges of terraces [@Evans2006]. The effects of a step barrier of purely kinetic origin, namely simple diffusion, were investigated in an epitaxial growth model with thermally activated diffusion [@Leal_JPCM]. In this model, a particle performing an interlayer movement through steps with more than one monolayer has to diffuse along the columns, perpendicularly to the substrate, instead of attaching directly at the bottom or top of a terrace. This kinetic barrier reduces downhill currents and three-dimensional structures in the form of mounds are obtained at short-time scales even in the case of weak ES barriers where the conventional rule would not lead to mound formation.

Simple models with limited mobility can be used to investigate kinetic roughening [@barabasi; @meakin]. Wolf-Villain (WV) [@WV] and Das Sarma-Tamborenea (DT) [@DT] models, introduced to investigate molecular-beam-epitaxy (MBE) growth, are benchmarks of this class and have been intensively investigated [@Smilauer; @Milan; @Huang; @HaselwandterPRL; @HaselwandterPRE; @Sarma; @Punyindu; @wvbogo; @Xun; @Luis2019]. A variation of the CV model with limited mobility has been considered [@Aarao2010; @Aarao2013] and many features of the original model have been reproduced with this simplified version [@To2018]. Effects of a step barrier were investigated in both WV [@Rangdee] and DT [@DasSarma_SC] models introducing two additional probabilities for downward and upward interlayer diffusion with the former larger than the latter, and mound formation was observed in both models. WV and DT models without step barrier were investigated in several lattices [@Chatraphorn2001; @Kanjanaput2010] and it was found that the WV model can present topologically induced mound morphologies on some lattices but not in others while no clear evidence for three-dimensional structures was observed for DT. In one-dimension, it is widely accepted that both DT and WV models asymptotically produce self-affine surfaces belonging to nonlinear MBE [@Luis2019] and Edwards-Wilkinson [@Vvedensky] universality classes, respectively.

It was reported that a kinetic barrier alone does not induce mound morphologies in thermally activated CV-like models [@Leal_JPCM] but, instead, they exhibit kinetic roughening with exponents consistent with the nonlinear MBE universality class [@DT; @Villain; @LSarma]. Therefore, given the simplicity of limited-mobility growth models and the non-trivial effects of topologically induced uphill currents in DT and WV models, one would wonder how they respond to a barrier of purely kinetic origin. In order to fill this gap, we investigate WV and DT models with the introduction of the kinetic barrier proposed in Ref. [@Leal_JPCM]. We observed mounds in both models in 1+1 and 2+1 dimensions, being much more evident for WV model. The surface coarsening ceases quickly with the saturation of the characteristic surface length and regimes of uncorrelated mound growth are asymptotically observed. Analysis of the out-of-plane currents shows a large reduction of the downhill flux of particles, enhancing surface instabilities and mound formation.

The remaining of the paper is organized as follows. The model implementation details are presented in section \[sec:model\]. In section \[results\], we discuss the results obtained in the simulations. Our conclusions and some perspectives are drawn in the section \[conclusion\].

Models {#sec:model}
======

In all investigated models, the particles are randomly deposited on a $d$-dimensional lattice of linear size $L$ with periodic boundary conditions under the SOS condition. Results presented in this work correspond to regular chains in $d=1$ and square lattices in $d=2$. Other lattices were tested and the central conclusions remain unaltered. The height of the interface at site $i$ and time $t$ is represented by $h_i(t)$ and the initial condition is given by $h_i(0)=0$ such that the initial interface is flat.

In the WV model with a kinetic barrier [investigated in the present work]{}, the growth rule is implemented as follows. At each time step, a position $i$ is randomly chosen. A location $i'$ with the largest number of bonds that a new deposited adatom would have is determined within a set containing $i$ and its nearest-neighbors. If the initial position corresponds to the largest number of bonds ($i'\equiv i$), it is chosen as the deposition place and the simulation runs to the next step. In case of multiple options, one is chosen at random. Otherwise, the particle tries to diffuse to the neighbor $i'$ with a probability given by [@Leal_JPCM] $$\label{prob}
P_{\delta h}(i,i')=
\left\{
\begin{array}{cl}
1, & \textrm{if } ~ |\delta h|<2\\ 
\frac {1}{|\delta h|}, & \textrm{if } ~ |\delta h|\geq 2 
\end{array}
\right.$$ where $\delta h = h_i-h_{i'}$. With probability $1 - P_{\delta h}(i,i')$ the particle remains at the site $i$. It is important to mention that Eq.  is obtained [ assuming that the adatom first moves to top kink of the terrace and then start a unbiased one-dimensional random-walk normally to the initial substrate, stopping the movement if it either arrives at the bottom or return to top of the terrace. The result is the]{} solution of a non-directed one-dimensional random walk with absorbing boundaries separated by a distance $|\delta h|$ [@Shehawey]; see Fig. 1 of Ref. [@Leal_JPCM] for further details of this diffusion rule. This diffusion attempt is successively applied $N_s$ times (representing a $N_s$ diffusive steps) departing from the last position of the adatom. A unit time is defined as the deposition of $L^d$ particles.

The implementation of the DT model with kinetic barrier is similar. The difference is that diffusion to the nearest-neighbors are performed only if the adatom does not have lateral bounds and any neighbor with a number of bonds higher than 1 can be chosen with equal chance as the target site.

![image](film_wvO_2d.pdf){width="0.32\linewidth"} ![image](film_wvNs1_2d.pdf){width="0.32\linewidth"} ![image](film_wvNs10_2d.pdf){width="0.32\linewidth"}\
![image](film_dtO_2d.pdf){width="0.32\linewidth"} ![image](film_dtNs1_2d.pdf){width="0.32\linewidth"} ![image](film_dtNs10_2d.pdf){width="0.32\linewidth"}

![\[hh\_corr\_time\]Main panels: Height-height correlation function for the WV model at distinct times indicated in the legends for (a) one- and (b) two-dimensional substrates. The number of steps is $N_s=1$. The averages were computed over 100 independent runs. Insets: Correlation functions averaged over 1 and 10 samples for the original WV model at time $t=10^5$ showing that the oscillations observed in single samples are not due to regular structures.](gamma_wv1d.pdf "fig:"){width="0.8\linewidth"}\
![\[hh\_corr\_time\]Main panels: Height-height correlation function for the WV model at distinct times indicated in the legends for (a) one- and (b) two-dimensional substrates. The number of steps is $N_s=1$. The averages were computed over 100 independent runs. Insets: Correlation functions averaged over 1 and 10 samples for the original WV model at time $t=10^5$ showing that the oscillations observed in single samples are not due to regular structures.](gamma_wv2dNs1.pdf "fig:"){width="0.8\linewidth"}

Results
=======

The one-dimensional simulations were carried out on chains with up to $L=2^{14}$ sites and evolution times of up to $t = 10^7$. In the two-dimensional case, the simulations were done in systems of size up to $L=2^{10}$ and time up to $t = 10^6$. The averages were performed over $100$ independent runs.

Figures \[snapshot\_1d\] and \[snapshot\_2d\] show interfaces obtained in simulations in one- and two-dimensional substrates, respectively. Surfaces for the original WV and DT models without and with ($N_s=1$ or $N_s=10$) kinetic barriers are compared. In both dimensions, the irregular morphologies without a characteristic length observed in the original versions change to structures separated by valleys that present a well-defined characteristic length. We also observe that an increase in the value of $N_s$ reduces valley deepness and increases the characteristic width of the mounds. The effects of the kinetic barrier seem to be stronger in two- than one-dimension. A remarkable change in the profiles happens when just one hop to nearest-neighbors is allowed in the DT model with kinetic barrier, as can been seen in Fig. \[snapshot\_1d\](e). Surfaces become columnar with a high aspect ratio (height/width). Such a behavior is reminiscent of the very strict rule for diffusion in DT when a single lateral bound is enough to irreversibly stick the adatom on a site. In the WV case, where diffusion happens more readily, mound morphologies with quasiregular structures emerge more clearly.

![\[hh\_corr\_Ns\] Main plot: Height-height correlation function dependence with the parameter $N_s$ (indicated in the legends) for the WV model in two-dimensional substrates at a time $t=10^5$. Inset: Same as the main plot for one dimension at a time $t=10^7$. Curves correspond to averages over 100 independent samples.](gamma_wv2d1d.pdf){width="0.8\linewidth"}

A standard tool to characterize the morphology of interfaces in growth process is the height-height correlation function defined as [@Evans2006; @Murty2003; @Chatraphorn2001] $$\Gamma({\mathbf}r) = \left\langle \tilde h({\mathbf}x) \tilde h({\mathbf}x + {\mathbf}r) \right\rangle_x,$$ here $\tilde h({\mathbf}x)$ is the height interface at position ${\mathbf}x$ relative to the mean height and $\left\langle \ldots \right\rangle_x$ denotes an average over the surface. The height-height correlation for ${\mathbf}r = 0$ is related to the interface width by $$\sqrt{\left\langle\Gamma(0)\right\rangle} = w$$ here $\left\langle\ldots\right\rangle$ denotes an average over independent runs. A self-affine interface is characterized by a height-height correlation function that goes monotonically to zero while those characterized by mounds exhibit oscillatory behavior around $0$. In the latter case, the first zero of $\Gamma(r)$, denoted by $\xi$, is a characteristic lateral length of the mounds in the surface.

Figure \[hh\_corr\_time\] shows the height-height correlation function for the WV model with kinetic barrier in one- and two-dimensional substrates. The curves clearly exhibit oscillatory behavior even for averages over 100 independent samples. Conversely, the irregular oscillatory behavior observed for the original WV model shown in insets of Fig. \[hh\_corr\_time\] is lumped after averaging. Therefore, interfaces obtained with kinetic barrier are characterized by the formation of quasiregular mound structures differently from those obtained using the original model that exhibits irregular structures within the intervals of size and time we investigated. These plots also show a coarsening of the mounds represented by the first minimum displacement at the early growth times.

![Main plot: Correlation function for the DT model in two-dimensional substrates for distinct times shown in the legends and fixed $N_s=10$. Inset: Correlation function for DT model in two dimensions at a fixed time $t=10^5$ and different values of $N_s$ shown in legends. Curves correspond to averages over 100 independent samples.[]{data-label="h_corr_DT"}](gamma_dt2d_ts_Ns.pdf){width="0.8\linewidth"}

The effect of the parameter $N_s$ in WV model is shown in Fig. \[hh\_corr\_Ns\]. As indicated by the interface profiles shown in Figs. \[snapshot\_1d\] and \[snapshot\_2d\], the characteristic lateral length increases with $N_s$ in both dimensions. The correlation function for DT model follows a qualitative similar dependence with $N_s$, as can be seen in Fig. \[h\_corr\_DT\] where the effects of time and number of diffusion steps in the correlation function of the DT model are shown. However, the mounds are much less evident than those obtained in the WV model. However, the correlation functions still present the typical oscillatory behavior of mounded structures that is preserved after the averaging over 100 independent samples. Besides, the typical width of the mounds in the DT model are much smaller than those of WV. It is important to note that the correlation function of the original DT model also presents an irregular behavior as does the WV model.

![\[width\]Time evolution of the interface width $w$ for WV (main panels) and DT (insets) models grown on (a) one- and (b) two-dimensional substrates. Both simulations with the kinetic barrier (using $N_s$ values indicated in the legend) and the original version are shown. In (a), dashed and solid lines are power-laws with exponents $3/8$ and $1/2$, respectively, in both main panels and insets. In (b), the slopes of the dashed and solid lines are $1/4$ and $1/2$, respectively.](w_wvdt_1d.pdf "fig:"){width="0.8\linewidth"}\
![\[width\]Time evolution of the interface width $w$ for WV (main panels) and DT (insets) models grown on (a) one- and (b) two-dimensional substrates. Both simulations with the kinetic barrier (using $N_s$ values indicated in the legend) and the original version are shown. In (a), dashed and solid lines are power-laws with exponents $3/8$ and $1/2$, respectively, in both main panels and insets. In (b), the slopes of the dashed and solid lines are $1/4$ and $1/2$, respectively.](w_wvdt_2d.pdf "fig:"){width="0.8\linewidth"}

Figure \[width\] shows the time evolution of the interface width for both models in one and two dimensions. The main panels and insets present the results for the WV and DT models, respectively, including or not the kinetic barrier. The interface width is expected to scale as $w\sim t^\beta$ where $\beta$ is the growth exponent [@barabasi]. The short time dynamics of both WV and DT models is well described by the linear version of the MBE equation [@Villain; @LSarma] $$\frac{\partial h}{\partial t} = -\nu \nabla^{4}h + \lambda \nabla^{2} (\nabla h)^{2} + \eta,
\label{Eq1}$$ with $\lambda=0$ where $\eta$ is a non-conservative Gaussian noise [@Villain; @LSarma; @DasSarma1992]. This result is confirmed in Fig. \[width\] where the short time behavior is consistent with the growth exponents $\beta=3/8$ in $d=1$ and $\beta=1/4$ in $d=2$ expected for the linear MBE universality class [@barabasi]. It is worth to mention that these models may undergo crossovers to different universality classes in the asymptotic, depending on the dimension and model [@wvbogo; @Vvedensky; @Xun; @Chen2017; @AaraoReis2004b; @Punyindu]. The curves in Fig. \[width\] are consistent with crossovers to different universality classes at long times. One expects that DT is asymptotically consistent with the non-linear MBE equation with $\lambda>0$ [@Xun; @Luis2019; @Luis2017], for which $\beta\approx 1/3$ and $1/5$ in $d=1$ and $d=2$, respectively[^1], while crossovers to the Edwards-Wilkinson universality class with $\beta=1/4$ in $d=1$ and $\beta=0$ (logarithmic growth) in $d=2$ are expected for the WV model [@wvbogo; @Vvedensky]. The simulations with the kinetic barrier, however, departs from the original dynamics after a transient which increases with the diffusion of particles. For long times, an evolution consistent with an uncorrelated growth [described by $\frac{\partial h}{\partial t}=\eta$, characterized by a growth exponent $\beta=1/2$ [@barabasi], is observed.]{} [This observation can be rationalized as follows. At long times, mounds interact weakly since the kinetic barrier reduces drastically inter-mound diffusion. Consider the idealized case of plateaus of size $L_0$ with an infinity barrier at their edges. A particle initially adsorbed on the top of a plateau will never slide down to its bottom. So, the probability that this plateau receives $R$ particles after one unity of time (deposition of $L$ particles) is a binomial distribution $$P(R)=\binom{L}{R} p^R (1-p)^{L-R}\simeq 
    \frac{1}{\sqrt{2\pi L_0}} e^{-\frac{(R-L_0)^2}{2L_0}},
    \label{eq:PR}$$ where $p=L_0/L$ is the probability that a particle is deposited on this terrace and $1\ll L_0\ll L$ is assumed in the Gaussian limit in right-hand side of Eq. . We argue that this situation is similar to the weakly interacting mound observed in our simulations.]{}

![Characteristic length of mounds $\xi$ for WV (main plots) and DT (insets) models with and without the kinetic barrier in (a) one- and (b) two-dimensional substrates for different values of the parameter $N_s$ indicated in the legend.\[xi\]](xi_wvdt_1d.pdf "fig:"){width="0.8\linewidth"}\
![Characteristic length of mounds $\xi$ for WV (main plots) and DT (insets) models with and without the kinetic barrier in (a) one- and (b) two-dimensional substrates for different values of the parameter $N_s$ indicated in the legend.\[xi\]](xi_wvdt_2d.pdf "fig:"){width="0.8\linewidth"}

In addition, as can be seen in Fig. \[xi\], the characteristic lateral lengths of simulations with kinetic barrier saturate after an initial transient in values that increase with the parameter $N_s$ while the models without barrier present coarsening with $\xi\sim t^{1/z}$ [@barabasi]. The saturation implies that the aspect ratio (height/width) of the mounds remains increasing with time and the surface does not present slope selection forming columnar growth. This property is also reflected in the asymptotic interface width scaling as $w\sim t^{1/2}$. As explained previously, it can be interpreted as an uncorrelated evolution of the columns, in which the 1/2 exponent comes out. The results shown in the insets of Figs. \[width\] and \[xi\] corroborate that the DT model presents the same behavior of the WV model despite of the mounds are less evident in the former.

![\[current\_1d\] Evolution of the out-of-plane current for (a) WV and (b) DT models grown in one-dimensional substrates. Models with the kinetic barrier using $N_s =1$, $2$ and $10$ steps (indicated in the legend) and the original version are shown. ](jz_wv1d.pdf "fig:"){width="0.8\linewidth"}\
![\[current\_1d\] Evolution of the out-of-plane current for (a) WV and (b) DT models grown in one-dimensional substrates. Models with the kinetic barrier using $N_s =1$, $2$ and $10$ steps (indicated in the legend) and the original version are shown. ](jz_dt1d.pdf "fig:"){width="0.8\linewidth"}

Instability and mound formation can be investigated considering the surface currents [@Siegert1994; @Krug1993]; see [@Krug1997] for details. In this work, we investigated the out-of-plane component of the current defined as [@Leal_Jstat] $$J_z = \frac{1}{N}\sum_{(i,j)}  \mathrm{sgn}(\delta h) D(i,j) P_{\delta h}(i,j)$$ where $\mathrm{sgn}(x)=1$ for $x>0$, $\mathrm{sgn}(x)=-1$ for $x<0$, and $\mathrm{sgn}(0)=0$ is the definition of sign function, $P_{\delta h}(i,j)$ is given by Eq. , and $D(i,j)$ is the rate of hopping attempts from site $i$ to $j$ and depend on the investigated model. The sum runs over all $N$ pairs of nearest-neighbors of the lattice. Let $n_i$ be the number of lateral bonds of site $i$ and $n^\text{max}_i$ the largest number of bonds among the nearest-neighbors of $i$. For the WV model, $D(i,j)$ is given by $$D(i,j)=
\left\{
\begin{array}{cl}
1/q_i^\text{WV}, & \textrm{if } ~ n_j = n^\text{max}_i \textrm{  and  } n_i < n^\text{max}_i\\ 
0,  & \textrm{otherwise.}
\end{array}
\right.,$$ where $q_i^\text{WV}$ is the number of nearest-neighbors with $n^\text{max}_i$ lateral bonds. We can express $D(i,j)$ for the DT as $$D(i,j)=
\left\{
\begin{array}{cl}
1/q_i^\text{DT}, & \textrm{if } ~ n_j > 0  \textrm{ and } n_i = 0\\ 
0,  & \textrm{otherwise.}
\end{array}
\right.,$$ where $q_i^{DT}$ is number of nearest-neighbors with at least one lateral bond. The quantity $J_z$ is the average interlayer diffusion rate per site.

=\_20pt\^3pt

  -- --------- --------- -------- --------------------
                                  
        WV        DT        WV             DT
      -0.0034   -0.0042   -0.015   -5$\times 10^{-5}$
      -0.0011   -0.0052   -0.019   -3$\times 10^{-4}$
      -0.014    -0.0053   -0.020   -5$\times 10^{-4}$
      -0.090    -0.050    -0.047         -0.030
  -- --------- --------- -------- --------------------

  : \[fit\_parameters\] Parameters $J_\infty$ obtained in the regression using the Eq. (\[current\_fit\]) in the in the last decade of data of the out-plane current curves ($t>10^6$ for $d=1$ and $t>10^5$ for $d=2$).

The currents for simulations in $d=1$ are presented in Fig. \[current\_1d\]. All versions in both $1+1$ and $2+1$ dimensions are characterized by a current with a downward (negative) flux with the intensity decreasing monotonically. Considering the last decade of time, we estimated the current $J_\infty$ for $t\rightarrow\infty$ using a regression with a simple allometric function in the form. $$J_z = J_\infty + a t^{-\gamma}\label{current_fit},$$ where $a$ and $\gamma$ are parameters. In all cases with step barrier, we obtained asymptotic small negative currents with a non-universal value of $\gamma$. The results can be seen in table \[fit\_parameters\]. The currents for the standard models are considerably larger than in the cases with barrier. The values for the DT model with barrier are very small indicating that this current could be actually null in the asymptotic limit as observed in thermally activated diffusion models with ES step barriers [@Leal_Jstat]. In the case of the WV model, the current values may indicate the same asymptotic behavior, but our present accuracy does not allow a conclusion on this issue.

Conclusions {#conclusion}
===========

In this work, we investigate the effects of a purely kinetic barrier caused by the out-of-plane step edge diffusion [@Leal_JPCM] on limited-mobility growth models. The cases of studies were the benchmark models of Wolf-Villain [@WV] and Das Sarma-Tamborenea [@DT]. Large-scale simulations were performed considering one- and two-dimensional substrates. It was observed that the introduction of the kinetic barrier induces the formation of quasiregular mound structures differently from those obtained with the original models that forms irregular (self-affine) structures in the interface. The kinetic barrier stabilizes the mound width, leading to the formation of quasiregular structures. The interface width in models with kinetic barriers has an initial regime similar to the original models. However, a growth exponent very close to $\beta = 1/2$ is observed for asymptotically long times. Also, the characteristic lateral length saturates after a transient that depends on the number of steps that an adatom can perform before irreversibly stick in a position. These results are consistent with mounds evolving independently. The dynamics in both one- and two-dimensional substrates are characterized by a strong reduction of downward current with respect to the original models. The downward flux have an intensity decreasing monotonically to a asymptotic value that seems to be null for DT model and small for WV, being the latter possibly still subject to strong crossover effects in the present analysis.

A central contribution of this work is to show that a very simple mechanism neglected in previous analysis, in which particles also diffuse in the direction perpendicular to the substrate, is able to change markedly the surface morphology of basic growth models with limited mobility. Our results are qualitatively very similar to those obtained when an explicit step barrier, with a smaller probability to move downward, is considered [@Rangdee]. Particularly, asymptotic mound morphology has been reported for limited mobility models in $d=2$ without barriers with the application of the noise reduction method [@Chatraphorn2001]. Our results corroborate this scenario since a small perturbation induces mound instability in this kind of processes while it alone does not produce mounds in models with thermally activated diffusion [@Leal_Jstat].

We expect that the concepts investigated in this work will be applied to more sophisticated models and aid the understanding of pattern formation in film growth and the production of self-assembled structures for technological applications.

SGA and SCF thank the financial support of Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq) and Fundação de Pesquisa do Estado de Minas Gerais (FAPEMIG). This study was financed in part by the Coordenação de Aperfeiçoamento de Pessoal de Nível Superior - Brasil (CAPES) - Finance Code 001.

[51]{} ifxundefined \[1\][ ifx[\#1]{} ]{} ifnum \[1\][ \#1firstoftwo secondoftwo ]{} ifx \[1\][ \#1firstoftwo secondoftwo ]{} ““\#1”” @noop \[0\][secondoftwo]{} sanitize@url \[0\][‘\
12‘\$12 ‘&12‘\#12‘12‘\_12‘%12]{} @startlink\[1\] @endlink\[0\] @bib@innerbibempty [**](https://www.springer.com/us/book/9783540407287) (, , ) [****,  ()](http://linkinghub.elsevier.com/retrieve/pii/S0167572906000021) [**](https://www.cambridge.org/core/books/fractal-concepts-in-surface-growth/0D9076FC287B60B2B1126BB165112F13) (, , ) [**](https://onlinelibrary.wiley.com/doi/abs/10.1002/%28SICI%291099-1204%28199911/12%2912%3A6%3C493%3A%3AAID-JNM346%3E3.0.CO%3B2-7) (, , ) [****,  ()](\doibase 10.1103/PhysRevLett.78.911) [****,  ()](\doibase 10.1103/PhysRevB.65.193407) [****,  ()](\doibase
        10.1103/PhysRevB.81.115462) [****,  ()](\doibase 10.1103/PhysRevLett.72.116) [****,  ()](\doibase
        10.1103/PhysRevLett.97.126101) [****, ()](\doibase 10.1103/PhysRevB.74.245410) [****,  ()](\doibase 10.1126/science.1159455) [****,  ()](\doibase
        10.1063/1.1726787) [****,  ()](\doibase 10.1063/1.1707904) [****, ()](http://link.aps.org/doi/10.1103/PhysRevE.82.041607) @noop [****,  ()]{} [****, ()](https://link.aps.org/doi/10.1103/PhysRevLett.82.3661) [****, ()](https://link.aps.org/doi/10.1103/PhysRevB.76.115417) [****,  ()](\doibase 10.1103/PhysRevLett.58.2235) [****,  ()](http://scitation.aip.org/content/aip/journal/jap/63/7/10.1063/1.341041) [****,  ()](http://stacks.iop.org/0953-8984/23/i=29/a=292201?key=crossref.9247eb5fac27043f6f6da8363d2658c5) [****,  ()](https://iopscience.iop.org/article/10.1209/0295-5075/13/5/002) [****,  ()](\doibase 10.1103/PhysRevLett.66.325) [****,  ()](\doibase
        10.1103/PhysRevB.49.5769) [****,  ()](\doibase 10.1103/PhysRevE.54.3933) [****,  ()](\doibase 10.1103/PhysRevE.54.5935) [****, ()](\doibase 10.1103/PhysRevLett.98.046102) [****,  ()](\doibase 10.1103/PhysRevE.77.061129) [****,  ()](\doibase
        10.1103/PhysRevE.65.036144) [****, ()](\doibase 10.1103/PhysRevE.57.R4863) [****,  ()](https://iopscience.iop.org/article/10.1088/1742-5468/2011/04/P04022/meta) [****, ()](\doibase 10.1103/PhysRevE.85.041126) [****, ()](https://link.aps.org/doi/10.1103/PhysRevE.99.022801) [****,  ()](\doibase 10.1103/PhysRevE.81.041605) [****,  ()](\doibase
        10.1103/PhysRevE.88.022128) [****,  ()](https://linkinghub.elsevier.com/retrieve/pii/S0378437118308884) [****,  ()](\doibase 10.1016/j.susc.2005.12.021) [****,  ()](\doibase 10.1016/S0039-6028(99)00209-5) [****, ()](http://link.aps.org/doi/10.1103/PhysRevB.64.205407) [****,  ()](\doibase
        10.1103/PhysRevE.68.010601) [****,  ()](\doibase
        10.1051/jp1:1991114) [****, ()](http://link.aps.org/doi/10.1103/PhysRevLett.66.2348) [****, ()](https://iopscience.iop.org/article/10.1088/0305-4470/33/49/301) [****,  ()](\doibase 10.1103/PhysRevLett.69.3762) [****,  ()](http://dx.doi.org/10.1016/j.physa.2016.08.057
        http://linkinghub.elsevier.com/retrieve/pii/S037843711630574X) [****, ()](https://link.aps.org/doi/10.1103/PhysRevE.70.031607) [****, ()](http://link.aps.org/doi/10.1103/PhysRevE.95.042801) [****, ()](http://link.aps.org/doi/10.1103/PhysRevLett.78.1082) [****, ()](http://link.aps.org/doi/10.1103/PhysRevLett.73.1517) [****, ()](https://link.aps.org/doi/10.1103/PhysRevLett.70.3271) [****,  ()](http://www.tandfonline.com/doi/abs/10.1080/00018739700101498) [****,  ()](http://stacks.iop.org/1742-5468/2011/i=09/a=P09018?key=crossref.72e1a82e5e15653a11bd7c44b2902efa)

[^1]: The exponents $\beta=1/3$ and $1/5$ are predictions of the one-loop renormalization group [@Villain; @LSarma]. Two-loop calculations [@Janssen], however, predict corrections where the growth exponents are slightly smaller than these values.
