Derivation of continuum stochastic equations for discrete growth models

Modeling of nonequilibrium surface growth by a limited-mobility model with distributed diffusion length

Kinetic Monte Carlo (KMC) simulations are a well-established numerical tool to investigate the time-dependent surface morphology in molecular beam epitaxy experiments. In parallel, simplified approaches such as limited mobility (LM) models characterized by a fixed diffusion length have been studied. Here we investigate an extended LM model to gain deeper insight into the role of diffusional processes concerning the growth morphology. Our model is based on the stochastic transition rules of the Das Sarma-Tamborena model but differs from the latter via a variable diffusion length. A first guess for this length can be extracted from the saturation value of the mean-squared displacement calculated from short KMC simulations. Comparing the resulting surface morphologies in the sub- and multilayer growth regime to those obtained from KMC simulations, we find deviations which can be cured by adding fluctuations to the diffusion length. This mimics the stochastic nature of particle diffusion on a substrate, an aspect which is usually neglected in LM models. We propose to add fluctuations to the diffusion length by choosing this quantity for each adsorbed particle from a Gaussian distribution, where the variance of the distribution serves as a fitting parameter. We show that the diffusional fluctuations have a huge impact on cluster properties during submonolayer growth as well as on the surface profile in the high coverage regime. The analysis of the surface morphologies on one- and two-dimensional substrates during sub- and multilayer growth shows that the LM model can produce structures that are indistinguishable to the ones from KMC simulations at arbitrary growth conditions.

Local roughness exponent in the nonlinear molecular-beam-epitaxy universality class in one dimension

We report local roughness exponents, α_{loc}, for three interface growth models in one dimension which are believed to belong to the nonlinear molecular-beam-epitaxy (nMBE) universality class represented by the Villain-Lais-Das Sarma (VLDS) stochastic equation. We applied an optimum detrended fluctuation analysis (ODFA) [Luis et al., Phys. Rev. E 95, 042801 (2017)2470-004510.1103/PhysRevE.95.042801] and compared the outcomes with standard detrending methods. We observe in all investigated models that ODFA outperforms the standard methods providing exponents in the narrow interval α_{loc}^{}∈[0.96,0.98] quantitatively consistent with two-loop renormalization group predictions for the VLDS equation. In particular, these exponent values are calculated for the Clarke-Vvdensky and Das Sarma-Tamborenea models characterized by very strong corrections to the scaling, for which large deviations of these values had been reported. Our results strongly support the absence of anomalous scaling in the nMBE universality class and the existence of corrections in the form α_{loc}^{}=1-ε of the one-loop renormalization group analysis of the VLDS equation.

Universality and dependence on initial conditions in the class of the nonlinear molecular beam epitaxy equation

We report extensive numerical simulations of growth models belonging to the nonlinear molecular beam epitaxy (nMBE) class, on flat (fixed-size) and expanding substrates (ES). In both d=1+1 and 2+1, we find that growth regime height distributions (HDs), and spatial and temporal covariances are universal, but are dependent on the initial conditions, while the critical exponents are the same for flat and ES systems. Thus, the nMBE class does split into subclasses, as does the Kardar-Parisi-Zhang (KPZ) class. Applying the "KPZ ansatz" to nMBE models, we estimate the cumulants of the 1+1 HDs. Spatial covariance for the flat subclass is hallmarked by a minimum, which is not present in the ES one. Temporal correlations are shown to decay following well-known conjectures.

Normal dynamic scaling in the class of the nonlinear molecular-beam-epitaxy equation

The scaling of local height fluctuations is studied numerically in lattice growth models of the class of the nonlinear stochastic equation of Villain-Lai-Das Sarma (VLDS) in substrate dimensions d=1 and 2. In d=1, the average local slopes of the conserved restricted solid-on-solid (CRSOS) models converge to a finite value in the long-time limit, with power-law corrections in time whose exponents are close to 0.1. Other VLDS models in d=1, such as that of Das Sarma and Tamborenea, show a divergence of local slopes up to 10(6) monolayers, typical of anomalous roughening, but a comparison of roughness distributions shows that they scale as the linear fourth-order growth equation in those time scales. Normal scaling is also obtained in a modified VLDS equation with instability suppression, in contrast to recent numerical works. In d=2, a CRSOS model and a model with lateral aggregation of diffusing particles show normal scaling of the local slopes, also with small correction exponents. These results consistently show that the VLDS class has normal dynamic scaling in d=1 and 2, in agreement with the theoretical predictions of Phys. Rev. Lett. 94, 166103 (2005), and they show that the apparently anomalous features observed in previous works are effects of large scaling correction terms or crossover effects.

Dynamic scaling in thin-film growth with irreversible step-edge attachment

We study dynamic scaling in a model with collective diffusion (CD) of isolated atoms in terraces and irreversible aggregation at step edges. Simulations are performed in two-dimensional substrates with several diffusion to deposition ratios R identical with D/F. Data collapse of scaled roughness distributions confirms that this model is in the class of the fourth-order nonlinear growth equation by Villain, Lai, and Das Sarma (VLDS) with negligible finite-size effects, while estimates of scaling exponents show some discrepancies. This result is consistent with the prediction of a recent renormalization group approach and improves previous numerical works on related models. The roughness follows dynamic scaling as W=Lalpha/R1/2f(xi/L), with correlation length xi=(Rt)1/z, where z is the dynamic exponent. We also propose a limited mobility (LM) model where the incident atom executes up to G steps before a new atom is adsorbed, and irreversibly aggregates at step edges. This model is also shown to belong to the VLDS class. The size of the plateaus in the film surface increases as G1/2 and the lateral correlation scales as G1/2t1/z. The time evolution of the roughness reproduces that of the CD model if an equivalent parameter G approximately R2/z is chosen. This suggests the possibility of using LM models with tunable diffusion length to simulate processes with simultaneous diffusion of many atoms. A scaling approach is used to justify exponent values and dynamic relations for both models, including the significant decrease of surface roughness as R or G increases.

Maximal- and minimal-height distributions of fluctuating interfaces

Maximal- and minimal-height distributions (MAHD, MIHD) of two-dimensional interfaces grown with the nonlinear equations of Kardar-Parisi-Zhang (KPZ, second order) and of Villain-Lai-Das Sarma (VLDS, fourth order) are shown to be different. Two universal curves may be MAHD or MIHD of each class depending on the sign of the relevant nonlinear term, which is confirmed by results of several lattice models in the KPZ and VLDS classes. The difference between MAHD and MIDH is connected with the asymmetry of the local height distribution. A simple, exactly solvable deposition-erosion model is introduced to illustrate this feature. The average extremal heights scale with the same exponent of the average roughness. In contrast to other correlated systems, generalized Gumbel distributions do not fit those MAHD and MIHD, nor those of Edwards-Wilkinson growth.

Finite-size effects in roughness distribution scaling

We study numerically finite-size corrections in scaling relations for roughness distributions of various interface growth models. The most common relation, which considers the average roughness <w(2)> as a scaling factor, is not obeyed in the steady states of a group of ballisticlike models in 2+1 dimensions, even when very large system sizes are considered. On the other hand, good collapse of the same data is obtained with a scaling relation that involves the root mean square fluctuation of the roughness, which can be explained by finite-size effects on second moments of the scaling functions. We also obtain data collapse with an alternative scaling relation that accounts for the effect of the intrinsic width, which is a constant correction term previously proposed for the scaling of <w(2)> . This illustrates how finite-size corrections can be obtained from roughness distributions scaling. However, we discard the usual interpretation that the intrinsic width is a consequence of high surface steps by analyzing data of restricted solid-on-solid models with various maximal height differences between neighboring columns. We also observe that large finite-size corrections in the roughness distributions are usually accompanied by huge corrections in height distributions and average local slopes, as well as in estimates of scaling exponents. The molecular-beam epitaxy model of Das Sarma and Tamborenea in 1+1 dimensions is a case example in which none of the proposed scaling relations work properly, while the other measured quantities do not converge to the expected asymptotic values. Thus although roughness distributions are clearly better than other quantities to determine the universality class of a growing system, it is not the final solution for this task.

Scaling in the crossover from random to correlated growth

In systems where deposition rates are high compared to diffusion, desorption, and other mechanisms that generate correlations, a crossover from random to correlated growth of surface roughness is expected at a characteristic time t0. This crossover is analyzed in lattice models via scaling arguments, with support from simulation results presented here and in other works. We argue that the amplitudes of the saturation roughness and of the saturation time t(x) scale as t0(1/2) and t0, respectively. For models with lateral aggregation, which typically are in the Kardar-Parisi-Zhang (KPZ) class, we show that t0 approximately p(-1), where p is the probability of the correlated aggregation mechanism to take place. However, t0 approximately p(-2) is obtained in solid-on-solid models with single-particle deposition attempts. This group includes models in various universality classes, with numerical examples being provided in the Edwards-Wilkinson (EW), KPZ, and Villain-Lai-Das Sarma (nonlinear molecular-beam epitaxy) classes. Most applications are for two-component models in which random deposition, with probability 1-p, competes with a correlated aggregation process with probability p. However, our approach can be extended to other systems with the same crossover, such as the generalized restricted solid-on-solid model with maximum height difference S, for large S. Moreover, the scaling approach applies to all dimensions. In the particular case of one-dimensional KPZ processes with this crossover, we show that t0 approximately nu(-1) and nu approximately lambda(2/3), where nu and lambda are the coefficients of the linear and nonlinear terms of the associated KPZ equations. The applicability of previous results to models in the EW and KPZ classes is discussed.

Langevin equations for fluctuating surfaces

Exact Langevin equations are derived for the height fluctuations of surfaces driven by the deposition of material from a molecular beam. We consider two types of model: deposition models, where growth proceeds by the deposition and instantaneous local relaxation of particles, with no subsequent movement, and models with concurrent random deposition and surface diffusion. Starting from a Chapman-Kolmogorov equation the deposition, relaxation, and hopping rules of these models are first expressed as transition rates within a master equation for the joint height probability density function. The Kramers-Moyal-van Kampen expansion of the master equation in terms of an appropriate "largeness" parameter yields, according to a limit theorem due to Kurtz [Stoch. Proc. Appl. 6, 223 (1978)], a Fokker-Planck equation that embodies the statistical properties of the original lattice model. The statistical equivalence of this Fokker-Planck equation, solved in terms of the associated Langevin equation, and solutions of the Chapman-Kolmogorov equation, as determined by kinetic Monte Carlo (KMC) simulations of the lattice transition rules, is demonstrated by comparing the surface roughness and the lateral height correlations obtained from the two formulations for the Edwards-Wilkinson [Proc. R. Soc. London Ser. A 381, 17 (1982)] and Wolf-Villain [Europhys. Lett. 13, 389 (1990)] deposition models, and for a model with random deposition and surface diffusion. In each case, as the largeness parameter is increased, the Langevin equation converges to the surface roughness and lateral height correlations produced by KMC simulations for all times, including the crossover between different scaling regimes. We conclude by examining some of the wider implications of these results, including applications to heteroepitaxial systems and the passage to the continuum limit.

Numerical study of roughness distributions in nonlinear models of interface growth

We analyze the shapes of roughness distributions of discrete models in the Kardar, Parisi, and Zhang (KPZ) and in the Villain, Lai, and Das Sarma (VLDS) classes of interface growth, in one and two dimensions. Three KPZ models in d=2 confirm the expected scaling of the distribution and show a stretched exponential tail approximately as exp(-x0.8), with a significant asymmetry near the maximum. Conserved restricted solid-on-solid models belonging to the VLDS class were simulated in d=1 and d=2. The tail in d=1 has the form exp(-x2) and, in d=2, has a simple exponential decay, but is quantitatively different from the distribution of the linear fourth-order (Mullins-Herring) theory. It is not possible to fit any of the above distributions to those of 1/f(alpha) noise interfaces, in contrast with recently studied models with depinning transitions.

Generating function, path integral representation, and equivalence for stochastic exclusive particle systems

We present the path integral representation of the generating function for classical exclusive particle systems. By introducing hard-core bosonic creation and annihilation operators and appropriate commutation relations, we construct the Fock space structure. Using the state vector, the generating function is defined and the master equation of the system is transformed into the equation for the generating function. Finally, the solution of the linear equation for the generating function is derived in the form of the path integral. Applying the formalism, the equivalence of reaction-diffusion processes of single species and two species is described.

What is the connection between ballistic deposition and the Kardar-Parisi-Zhang equation?

Ballistic deposition (BD) is believed to belong to the Kardar-Parisi-Zhang (KPZ) universality class. In this paper we study the validity of this belief by rigorously deriving a continuum equation from the BD microscopic rules, which deviates from the KPZ equation. We show that in one dimension and in the presence of noise the deviation is not important. This is not the case in the absence of noise. In more than one dimension and in the presence of noise we obtain an equation that superficially seems to be a continuum equation but in which the symmetry under rotations around the growth direction is broken.

Numerical study of discrete models in the class of the nonlinear molecular beam epitaxy equation

We study numerically some discrete growth models belonging to the class of the nonlinear molecular beam epitaxy equation, or the Villain-Lai-Das Sarma (VLDS) equation. The conserved restricted solid-on-solid model (CRSOS) with maximum height differences Delta H(max)=1 and Delta H(max)=2 was analyzed in substrate dimensions d=1 and d=2 . The Das Sarma and Tamborenea (DT) model and a competitive model involving random deposition and CRSOS deposition were studied in d=1. For the CRSOS model with Delta H(max)=1, we obtain the more accurate estimates of scaling exponents in d=1:roughness exponent alpha=0.94+/-0.02 and dynamical exponent z=2.88+/-0.04. These estimates are significantly below the values of one-loop renormalization for the VLDS theory, which confirms Janssen's proposal of the existence of higher-order corrections. The roughness exponent in d=2 is very near the one-loop result alpha=2/3, in agreement with previous works. The moments W(n) of orders n=2 , 3, 4 of the height distribution were calculated for all models, and the skewness S triple bond W3/W(3/2)(2) and the kurtosis Q triple bond W4/W(2)2-3 were estimated. At the steady states, the CRSOS models and the competitive model have nearly the same values of S and Q in d=1, which suggests that these amplitude ratios are universal in the VLDS class. The estimates for the DT model are different, possibly due to their typically long crossover to asymptotic values. Results for the CRSOS models in d=2 also suggest that those quantities are universal.

Probing universality classes in solid-on-solid deposition

We consider several stochastic processes corresponding to the same physical solid-on-solid deposition problem. Simplified models presenting the same (conditional) mean and variance for each process are also introduced as well as generalizations in terms of the deposition of blobs and probabilistic deposition rules. We compare the evolution of the roughness as a function of time for a three-parameter family that includes as limit cases the Family model and the Edwards-Wilkinson equation, showing that in all cases the derived models with the same mean and variance are indistinguishable from the originating models in terms of the evolution of the roughness. Finally, we show that although all the models studied belong to the same universality class, some relevant features such as the final surface roughness are reproduced only for models within a restricted class determined by sharing the same (conditional) mean and variance.

Edwards-Wilkinson equation from lattice transition rules

Continuum equations of motion for the height fluctuations of lattice growth models are derived from their transition rules by regularizing and coarse-graining the associated discrete Langevin equations. For models with random deposition followed by instantaneous relaxation to a neighboring site based on identifying the local height minimum, our methodology yields the Edwards-Wilkinson equation. The application of this procedure to other growth models is discussed.

Universality classification of restricted solid-on-solid type surface growth models

We consider the restricted solid-on-solid (RSOS) type surface growth models and classify them into dynamic universality classes according to their symmetry and conservation law. Four groups of RSOS-type microscopic models--asymmetric (A), asymmetric-conserved (AC), symmetric (S), and symmetric-conserved (SC) groups--are introduced and the corresponding stochastic differential equations (SDEs) are derived. Analyzing these SDEs using dynamic renormalization group theory, we confirm the previous results that A-RSOS, AC-RSOS, and S-RSOS groups belong to the Kardar-Parisi-Zhang class, the Villain-Lai-Das Sarma class, and the Edwards-Wilkinson class, respectively. We also find that SC-RSOS group belongs to a new universality class featuring the conserved-cubic nonlinearity.