Persistence probabilities of the interface height in ( 1+1 ) - and ( 2+1 ) -dimensional atomistic, solid-on-solid, stochastic models of surface growth are studied using kinetic Monte Carlo simulations, with emphasis on models that belong to the molecular beam epitaxy (MBE) universality class. Both the initial transient and the long-time steady-state regimes are investigated. We show that for growth models in the MBE universality class, the nonlinearity of the underlying dynamical equation is clearly reflected in the difference between the measured values of the positive and negative persistence exponents in both transient and steady-state regimes. For the MBE universality class, the positive and negative persistence exponents in the steady-state are found to be theta(S)(+) =0.66+/-0.02 and theta(S)(-) =0.78+/-0.02, respectively, in ( 1+1 ) dimensions, and theta(S)(+) =0.76+/-0.02 and theta(S)(-) =0.85+/-0.02, respectively, in ( 2+1 ) dimensions. The noise reduction technique is applied on some of the ( 1+1 ) -dimensional models in order to obtain accurate values of the persistence exponents. We show analytically that a relation between the steady-state persistence exponent and the dynamic growth exponent, found earlier to be valid for linear models, should be satisfied by the smaller of the two steady-state persistence exponents in the nonlinear models. Our numerical results for the persistence exponents are consistent with this prediction. We also find that the steady-state persistence exponents can be obtained from simulations over times that are much shorter than that required for the interface to reach the steady state. The dependence of the persistence probability on the system size and the sampling time is shown to be described by a simple scaling form.

We investigate, using the noise reduction technique, the asymptotic universality class of the well-studied nonequilibrium limited mobility atomistic solid-on-solid surface growth models introduced by Wolf and Villain (WV) and Das Sarma and Tamborenea (DT) in the context of kinetic surface roughening in ideal molecular beam epitaxy. We find essentially all the earlier conclusions regarding the universality class of DT and WV models to be severely hampered by slow crossover and extremely long-lived transient effects. We identify the correct asymptotic universality class(es) that differs from earlier conclusions in several instances.

We study, through large-scale stochastic simulations using the noise reduction technique, surface growth via vapor deposition, e.g., molecular beam epitaxy ͑MBE͒, for simple nonequilibrium limited mobility solid-onsolid growth models, such as the Family model, the Das Sarma-Tamborenea model, the Wolf-Villain ͑WV͒ model, the larger-curvature ͑LC͒ model, and other related models. We find that (dϭ2ϩ1)-dimensional surface growth in several noise reduced models ͑most notably the WV and LC models͒ exhibits spectacular quasiregular mound formation with slope selection in their dynamical surface morphology in contrast to the standard statistically scale-invariant kinetically rough surface growth expected ͑and earlier reported in the literature͒ for such growth models. The mounding instability in these epitaxial growth models does not involve the Ehrlich-Schwoebel step-edge diffusion barrier. The mounded morphology in these growth models arises from the interplay between the line tension along step edges in the plane parallel to the average surface and the suppression of noise and island nucleation. The line tension tends to stabilize some of the step orientations that coincide with in-plane high-symmetry crystalline directions, and thus the mounds that are formed assume quasiregular structures. The noise reduction technique developed originally for Eden-type models can be used to control the stochastic noise and enhance diffusion along the step edge, which ultimately leads to the formation of quasiregular mounds during growth. We show that by increasing the diffusion surface length together with supression of nucleation and deposition noise, one can obtain a self-organization of the pyramids in quasiregular patterns. The mounding instability in these simple epitaxial growth models is closely related to the cluster-edge diffusion ͑as opposed to step-edge barrier͒ driven mounding in MBE growth, which has been recently discussed in the literature. The epitaxial mound formation studied here is a kinetic-topological instability ͓which can happen only in (dϭ2ϩ1)-dimensional, or higher dimensional, growth, but not in (dϭ1 ϩ1)-dimensional growth because no cluster diffusion around a closed surface loop is possible in ''onedimensional'' surfaces͔, which is likely to be quite generic in real MBE-type surface growth. Our extensive numerical simulations produce mounded ͑and slope-selected͒ surface growth morphologies which are strikingly visually similar to many recently reported experimental MBE growth morphologies.

We study, using noise reduction techniques, layer by layer epitaxial growth in limited mobility solid-on-solid nonequilibrium surface growth models, which have been introduced in the context of kinetic surface roughening in ideal molecular beam epitaxy. Multiple hit noise reduction and long surface diffusion length lead to qualitatively similar layer by layer epitaxy in 1+1 and 2+1 dimensional limited mobility growth simulations. We discuss the dynamic scaling characteristics connecting the transient layer by layer growth regime with the asymptotic kinetically rough growth regime.

Wolf-Villain (WV) model is a simple model used to study ideal molecular beam epitaxy (MBE) growth by using computer simulations. In this model, an adatom diffuses instantaneously within a finite diffusion length to maximize its coordination number. We study statistical properties of thin films grown by this model. The morphology of the WV model is found to be kinetically rough with a downhill particle diffusion current. In real MBE growth, however, there are additional factors such as the existence of a potential barrier that is known as the Ehrlich-Schwoebel (ES) barrier. The ES barrier is an additional barrier for an adatom that diffuses over a step edge from the upper to a lower terrace which is known to induce an uphill particle current. We found that with the addition of the ES barrier, the WV-ES model morphology is rough with mound formation on the surface when the barrier is strong enough. To confirm these results, the correlation function is also studied. We find no oscillation in the correlation function in the WV model. For the WV-ES model, the correlation function oscillates. These results confirm that a strong enough ES barrier can cause mound formation on the WV surface in our study.

The effects of the initial height on the temporal persistence probability of steady-state height fluctuations in up-down symmetric linear models of surface growth are investigated. We study the (1+1)-dimensional Family model and the (1+1)- and (2+1)-dimensional larger curvature (LC) model. Both the Family and LC models have up-down symmetry, so the positive and negative persistence probabilities in the steady state, averaged over all values of the initial height h(0), are equal to each other. However, these two probabilities are not equal if one considers a fixed nonzero value of h(0). Plots of the positive persistence probability for negative initial height versus time exhibit power-law behavior if the magnitude of the initial height is larger than the interface width at saturation. By symmetry, the negative persistence probability for positive initial height also exhibits the same behavior. The persistence exponent that describes this power-law decay decreases as the magnitude of the initial height is increased. The dependence of the persistence probability on the initial height, the system size, and the discrete sampling time is found to exhibit scaling behavior.

No abstract

The energetically driven Ehrlich-Schwoebel barrier had been generally accepted as the primary cause of the growth instability in the form of quasiregular moundlike structures observed on the surface of thin film grown via molecular-beam epitaxy (MBE) technique. Recently the second mechanism of mound formation was proposed in terms of a topologically induced flux of particles originating from the line tension of the step edges which form the contour lines around a mound. Through large-scale simulations of MBE growth on a variety of crystalline lattice planes using limited-mobility, solid-on-solid models introduced by Wolf-Villain and Das Sarma-Tamborenea in 2+1 dimensions, we show that there exists a topological uphill particle current with strong dependence on specific lattice crystalline structure. Without any energetically induced barriers, our simulations produce spectacular mounds very similar, in some cases, to what have been observed in many recent MBE experiments. On a lattice where these currents cease to exist, the surface appears to be scale invariant, statistically rough as predicted by the conventional continuum growth equation.