The process of irreversible nucleation in multilayer growth. II. Exact results in one and two dimensions

We study irreversible dimer nucleation on top of terraces during epitaxial growth in one and two dimensions, for all values of the step-edge barrier. The problem is solved exactly by transforming it into a first passage problem for a random walker in a higher-dimensional space. The spatial distribution of nucleation events is shown to differ markedly from the mean-field estimate except in the limit of very weak step-edge barriers. The nucleation rate is computed exactly, including numerical prefactors.Comment: 22 pages, 10 figures. To appear in Phys. Rev.

Determination of step--edge barriers to interlayer transport from surface morphology during the initial stages of homoepitaxial growth

We use analytic formulae obtained from a simple model of crystal growth by molecular--beam epitaxy to determine step--edge barriers to interlayer transport. The method is based on information about the surface morphology at the onset of nucleation on top of first--layer islands in the submonolayer coverage regime of homoepitaxial growth. The formulae are tested using kinetic Monte Carlo simulations of a solid--on--solid model and applied to estimate step--edge barriers from scanning--tunneling microscopy data on initial stages of Fe(001), Pt(111), and Ag(111) homoepitaxy.Comment: 4 pages, a Postscript file, uuencoded and compressed. Physical Review B, Rapid Communications, in press

Relaxation kinetics in two-dimensional structures

We have studied the approach to equilibrium of islands and pores in two dimensions. The two-regime scenario observed when islands evolve according to a set of particular rules, namely relaxation by steps at low temperature and smooth at high temperature, is generalized to a wide class of kinetic models and the two kinds of structures. Scaling laws for equilibration times are analytically derived and confirmed by kinetic Monte Carlo simulations.Comment: 6 pages, 7 figures, 1 tabl

Diffusion on a stepped substrate

We present results for collective diffusion of adatoms on a stepped substrate with a submonolayer coverage. We study the combined effect of the additional binding at step edge, the Schwoebel barrier, the enhanced diffusion along step edges, and the finite coverage on diffusion as a function of step density. In particular, we examine the crossover from step--dominated diffusion at high step density to terrace-dominated behavior at low step density in a lattice-gas model using analytical Green's function techniques and Monte Carlo simulations. The influence of steps on diffusion is shown to be more pronounced than previously anticipated.Comment: 4 pages, RevTeX, 3 Postscript figure

Coarsening of Surface Structures in Unstable Epitaxial Growth

We study unstable epitaxy on singular surfaces using continuum equations with a prescribed slope-dependent surface current. We derive scaling relations for the late stage of growth, where power law coarsening of the mound morphology is observed. For the lateral size of mounds we obtain $\xi \sim t^{1/z}$ with $z \geq 4$. An analytic treatment within a self-consistent mean-field approximation predicts multiscaling of the height-height correlation function, while the direct numerical solution of the continuum equation shows conventional scaling with z=4, independent of the shape of the surface current.Comment: 15 pages, Latex. Submitted to PR

The process of irreversible nucleation in multilayer growth. I. Failure of the mean-field approach

The formation of stable dimers on top of terraces during epitaxial growth is investigated in detail. In this paper we focus on mean-field theory, the standard approach to study nucleation. Such theory is shown to be unsuitable for the present problem, because it is equivalent to considering adatoms as independent diffusing particles. This leads to an overestimate of the correct nucleation rate by a factor N, which has a direct physical meaning: in average, a visited lattice site is visited N times by a diffusing adatom. The dependence of N on the size of the terrace and on the strength of step-edge barriers is derived from well known results for random walks. The spatial distribution of nucleation events is shown to be different from the mean-field prediction, for the same physical reason. In the following paper we develop an exact treatment of the problem.Comment: 19 pages, 3 figures. To appear in Phys. Rev.

Changing shapes in the nanoworld

What are the mechanisms leading to the shape relaxation of three dimensional crystallites ? Kinetic Monte Carlo simulations of fcc clusters show that the usual theories of equilibration, via atomic surface diffusion driven by curvature, are verified only at high temperatures. Below the roughening temperature, the relaxation is much slower, kinetics being governed by the nucleation of a critical germ on a facet. We show that the energy barrier for this step linearly increases with the size of the crystallite, leading to an exponential dependence of the relaxation time.Comment: 4 pages, 5 figures. Accepted by Phys Rev Let

Lattice Effects in Crystal Evaporation

We study the dynamics of a stepped crystal surface during evaporation, using the classical model of Burton, Cabrera and Frank, in which the dynamics of the surface is represented as a motion of parallel, monoatomic steps. The validity of the continuum approximation treated by Frank is checked against numerical calculations and simple, qualitative arguments. The continuum approximation is found to suffer from limitations related, in particular, to the existence of angular points. These limitations are often related to an adatom detachment rate of adatoms which is higher on the lower side of each step than on the upper side ("Schwoebel effect").Comment: DRFMC/SPSMS/MDN, Centre d'Etudes Nucleaires de Grenoble, 25 pages, LaTex, revtex style. 8 Figures, available upon request, report# UBFF30119

Growth of Patterned Surfaces

During epitaxial crystal growth a pattern that has initially been imprinted on a surface approximately reproduces itself after the deposition of an integer number of monolayers. Computer simulations of the one-dimensional case show that the quality of reproduction decays exponentially with a characteristic time which is linear in the activation energy of surface diffusion. We argue that this life time of a pattern is optimized, if the characteristic feature size of the pattern is larger than $(D/F)^{1/(d+2)}$, where $D$ is the surface diffusion constant, $F$ the deposition rate and $d$ the surface dimension.Comment: 4 pages, 4 figures, uses psfig; to appear in Phys. Rev. Let