Growth equations for the Wolf-Villain and Das Sarma-Tamborenea models of molecular-beam epitaxy

Bifractality in the one-dimensional Wolf-Villain model

We introduce a multifractal optimal detrended fluctuation analysis to study the scaling properties of the one-dimensional Wolf-Villain (WV) model for surface growth. This model produces coarsened surface morphologies for long timescales (up to 10^{9} monolayers) and its universality class remains an open problem. Our results for the multifractal exponent τ(q) reveal an effective local roughness exponent consistent with a transient given by the molecular beam epitaxy (MBE) growth regime and Edwards-Wilkinson (EW) universality class for negative and positive q values, respectively. Therefore, although the results corroborate that long-wavelength fluctuations belong to the EW class in the hydrodynamic limit, as conjectured in the recent literature, a bifractal signature of the WV model with an MBE regime at short wavelengths was observed.

Effects of a kinetic barrier on limited-mobility interface growth models

The role played by a kinetic barrier originated by out-of-plane step edge diffusion, introduced by Leal et al. [J. Phys.: Condens. Matter 23, 292201 (2011)JCOMEL0953-898410.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 β=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.

Temperature behaviour of the average size of nanoparticle lattices co-deposited with an amorphous matrix. Analysis of Ge + Al2O3 and Ni + Al2O3 thin films

We theoretically interpret the thermal behaviour of the average radius versus substrate temperature of regular quantum dot/nanocluster arrays formed by sputtering semiconductor/metal atoms with oxide molecules. The analysis relies on a continuum theory for amorphous films with given surface quantities, perturbed by a nanoparticle lattice. An account of the basic thermodynamic contributions is given in terms of force-flux phenomenological coefficients of each phase (Ge, Ni, Al₂O₃). Average radii turn out to be expressible by a characteristic length scale and a dimensionless parameter, which mainly depend upon temperature through diffusion lengths, film pressures and finite-size corrections to interfacial tensions. The numerical agreement is good in both Ge ([Formula: see text]) and Ni ([Formula: see text]) lattices grown at temperatures [Formula: see text]800 K, despite the lower temperature behaviour of quantum dots seeming to suggest further driving forces taking part in such processes.

Asymptotic dynamic scaling behavior of the (1+1)-dimensional Wolf-Villain model

Extensive kinetic Monte Carlo simulations are presented for the Wolf-Villain model in (1+1) dimensions. Asymptotic dynamic scaling is found for lattice sizes L≥2048. The exponents obtained from our simulations, α=0.50±0.02 and β=0.25±0.02, are in excellent agreement with the exact values α=1/2 and β=1/4 for the one-dimensional Edwards-Wilkinson equation. Our findings explain the widespread discrepancies of previous reports for exponents of the Wolf-Villain model in (1+1) dimensions, and the results are also consistent with the theoretical predictions of López et al. [J. M. López, M. Castro, and R. Gallego, Phys. Rev. Lett. 94, 166103 (2005)].

Extensive studies on linear growth processes with spatiotemporally correlated noise in arbitrary substrate dimensions

An extensive analytical and numerical study on a class of growth processes with spatiotemporally correlated noise in arbitrary dimension is undertaken. In addition to the conventional investigation on the interface morphology and interfacial widths, we pay special attention to exploring the characteristics of the slope-slope correlation function S(r,t) and the [Q]-th degree residual local interfacial width w[Q](l,t), whose importance has been somewhat overlooked in the literature. Based on the above analysis, we give a plausible theoretical explanation about the various experimental observations of kinetically and thermodynamically unstable surface growth. Furthermore, through explicit examples, we show that the statistical methods of calculating the exponents (including the dynamic exponent z, the global roughness exponent α, and the local roughness exponent α(loc)), based on the scaling of S(r,t) and w[Q](l,t), are very reliable and rarely influenced by the finite time and/or finite-size effects. Another important issue we focus on in this paper is related to numerical calculation. For the specific class of growth processes discussed in this paper, we develop a very efficient and accurate algorithm for numerical calculation of the dynamics of interface configuration, the structure factor, the various correlation functions, the interfacial width and its variants in arbitrary dimensions, even with very large system size and very late time. The proposed systematical algorithm can be easily generalized to other linear processes and some special nonlinear processes.

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.

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.

Crossover and universality in the Wolf-Villain model

The transition rules of the Wolf-Villain model for the deposition and instantaneous relaxation of particles on a lattice are expressed as a Langevin equation for the height fluctuations at each site. A coarse-graining transformation of this equation shows directly that this model belongs to the Edwards-Wilkinson universality class, in agreement with kinetic Monte Carlo simulations. The crossover from the Mullins-Herring equation is explained by the transformation under coarse graining of the coefficients in the equation of motion.

Theoretical continuous equation derived from the microscopic dynamics for growing interfaces in quenched media

We present an analytical continuous equation for the Tang and Leschhorn model [Phys. Rev. A 45, R8309 (1992)] derived from their microscopic rules using a regularization procedure. As well in this approach, the nonlinear term (nablah)(2) arises naturally from the microscopic dynamics even if the continuous equation is not the Kardar-Parisi-Zhang equation [Phys. Rev. Lett. 56, 889 (1986)] with quenched noise (QKPZ). Our equation is similar to a QKPZ equation but with multiplicative quenched and thermal noise. The numerical integration of our equation reproduces all the scaling exponents of the directed percolation depinning model.

Effect of symmetry on volume conserving surface models

We study the effect of symmetry on volume conserving models without deposition and evaporation. By using the master equation approach, we identify two types of stochastic continuum equation with a conservative noise, depending on the symmetry of hopping rate in diffusion rules. In the model with symmetric hopping rate, a Laplacian term is essentially absent from the continuum equation. The dynamic scaling of this model is thus determined by the nonlinear fourth order equation with a conservative noise. When the symmetry is broken, a Laplacian term may be present, so the asymptotic scaling behavior is governed by the Laplacian term with nonzero coefficient. We verify this result by investigating a simple discrete model analytically.