Ehrlich-Schwoebel instability in molecular-beam epitaxy: A minimal model

Epitaxial growth in one dimension

The final structure and properties of layers grown by epitaxy techniques are determined in the very early stage of the process. This review describes one-dimensional models for epitaxial growth, emphasizing the basic theoretical concepts employed to analyze nucleation and aggregation phenomena in the submonolayer regime. The main findings regarding the evolution of quantities that define the properties of the system, such as monomer and island densities, and the associated island size, gap length, and capture zone distributions are discussed, as well as the analytical tools used to evaluate them. This review provides a concise overview of the most widely used algorithms for simulating growth processes, discusses relevant experimental results, and establishes connections with existing theoretical studies.

Colloidal model for nucleation and aggregation in one dimension: Accessing the interaction parameters

Through a one-dimensional colloidal model for epitaxial growth, we characterize the nucleation and aggregation processes occurring in a gap between adjacent islands. The timescales associated with deposition, diffusion, aggregation, and nucleation inside the gap are studied in terms of the parameters defining the interaction between colloidal particles. Numerical results from molecular-dynamics (MD) simulations are compared with analytical models and good agreement is found between both data sets. The results for the timescales are used to calculate the associated rates to generate kinetic Monte Carlo (KMC) simulations, which allow exploring larger systems and longer timescales in comparison with MD simulations. The KMC simulations reproduce the global behavior of the densities of islands and monomers as well as the gap length distribution between adjacent islands.

Island density scaling in reversible submonolayer growth with attachment barriers

Island nucleation and growth play an important role in thin-film growth. One quantity of particular interest is the exponent χ, which describes the dependence of the saturation island density N_{sat}∼(D_{h}/F)^{-χ} on the ratio D_{h}/F of the monomer hopping rate D_{h} to the deposition rate F. While standard rate equation (RE) theory predicts that χ=i/(i+2) (where i is the critical island size), more recently it has been predicted that in the presence of a strong barrier to the attachment of monomers to islands, a significantly larger value χ=2i/(i+3) may be observed. While this prediction has recently been tested using kinetic Monte Carlo simulations for the case of irreversible growth corresponding to i=1, it has not been tested for the case of reversible island growth corresponding to i>1. Here we present a mean-field self-consistent RE method which we have used to study the dependence of the effective value of χ on D_{h}/F and barrier-strength for i=1,2,3, and 6. Both the no-nucleation-barrier case in which there exists a barrier for monomers to attach to islands larger than the critical island size (but not to smaller islands) and the nucleation-barrier case in which there is a barrier for monomers to attach to islands of all sizes are studied. In all cases, we find that the existence of attachment barriers significantly increases the effective value of χ for a given barrier strength. In addition, for i=1 we find good agreement between our extrapolated asymptotic value of χ and the theoretical strong-barrier prediction both with and without a nucleation barrier. In contrast, for i>1 the value of χ is significantly larger in the presence of a nucleation barrier than in its absence. In particular, while an asymptotic analysis of our results for i>1 also leads to excellent agreement with the strong barrier prediction in the presence of a nucleation barrier, in the absence of a nucleation barrier the asymptotic values are significantly lower.

Interface collisions with diffusive mass transport

We report on a linear Langevin model that describes the evolution of the roughness of two interfaces that move towards each other and are coupled by a diffusion field. This model aims at describing the closing of the gap between two 2D material domains during growth, and the subsequent formation of a rough grain boundary. We assume that deposition occurs in the gap between the two domains and that the growth units diffuse and may attach to the edges of the domains. These units can also detach from edges, diffuse, and reattach elsewhere. For slow growth, the edge roughness increases monotonously and then saturates at some equilibrium value. For fast growth, the roughness exhibits a maximum just before the collision between the two interfaces, which is followed by a minimum. The peak of the roughness can be dominated by statistical fluctuations or by edge instabilities. A phase diagram with three regimes is obtained: Slow growth without peak, peak dominated by statistical fluctuations, and peak dominated by instabilities. These results reproduce the main features observed in kinetic Monte Carlo simulations.

Plasmonic metal nanostructures with extremely small features: new effects, fabrication and applications

Surface plasmons in metals promise many fascinating properties and applications in optics, sensing, photonics and nonlinear fields. Plasmonic nanostructures with extremely small features especially demonstrate amazing new effects as the feature sizes scale down to the sub-nanometer scale, such as quantum size effects, quantum tunneling, spill-out of electrons and nonlocal states etc. The unusual physical, optical and photo-electronic properties observed in metallic structures with extreme feature sizes enable their unique applications in electromagnetic field focusing, spectra enhancing, imaging, quantum photonics, etc. In this review, we focus on the new effects, fabrication and applications of plasmonic metal nanostructures with extremely small features. For simplicity and consistency, we will focus our topic on the plasmonic metal nanostructures with feature sizes of sub-nanometers. Subsequently, we discussed four main and typical plasmonic metal nanostructures with extremely small features, including: (1) ultra-sharp plasmonic metal nanotips; (2) ultra-thin plasmonic metal films; (3) ultra-small plasmonic metal particles and (4) ultra-small plasmonic metal nanogaps. Additionally, the corresponding fascinating new effects (quantum nonlinear, non-locality, quantum size effect and quantum tunneling), applications (spectral enhancement, high-order harmonic wave generation, sensing and terahertz wave detection) and reliable fabrication methods will also be discussed. We end the discussion with a brief summary and outlook of the main challenges and possible breakthroughs in the field. We hope our discussion can inspire the broader design, fabrication and application of plasmonic metal nanostructures with extremely small feature sizes in the future.

Two-step unconventional protocol for epitaxial growth in one dimension with hindered reactions

We study the effect of hindered aggregation and/or nucleation on the island formation process in a two-step growth protocol. In the proposed model, the attachment of monomers to islands and/or other monomers is hindered by additional energy barriers which decrease the hopping rate of the monomers to the occupied sites of the lattice. For zero and weak barriers, the attachment is limited by diffusion while for strong barriers it is limited by reaction. We describe the time evolution of the system in terms of the monomer and island densities, N_{1} and N. We also calculate the gap length, the capture zone and the island distributions. For all the sets of barriers considered, the results given by the proposed analytical model are compared with those from kinetic Monte Carlo simulations. We found that the behavior of the system depends on the ratio of the nucleation barrier to the aggregation barrier. The two-step growth protocol allows more control and understanding on the island formation mechanism because it intrinsically separates the nucleation and aggregation processes in different time regimes.

Coherent X-ray measurement of step-flow propagation during growth on polycrystalline thin film surfaces

The properties of artificially grown thin films are strongly affected by surface processes during growth. Coherent X-rays provide an approach to better understand such processes and fluctuations far from equilibrium. Here we report results for vacuum deposition of C₆₀ on a graphene-coated surface investigated with X-ray Photon Correlation Spectroscopy in surface-sensitive conditions. Step-flow is observed through measurement of the step-edge velocity in the late stages of growth after crystalline mounds have formed. We show that the step-edge velocity is coupled to the terrace length, and that there is a variation in the velocity from larger step spacing at the center of crystalline mounds to closely-spaced, more slowly propagating steps at their edges. The results extend theories of surface growth, since the behavior is consistent with surface evolution driven by processes that include surface diffusion, the motion of step-edges, and attachment at step edges with significant step-edge barriers.

Designing Nanostructured Ti6Al4V Bioactive Interfaces with Directed Irradiation Synthesis toward Cell Stimulation to Promote Host-Tissue-Implant Integration

A new generation of biomaterials are evolving from being biologically inert toward bioactive surfaces, which can further interact with biological components at the nanoscale. Here, we present directed irradiation synthesis (DIS) as a novel technology to selectively apply plasma ions to bombard any type of biomaterial and tailor the nanofeatures needed for in vitro growth stimulation. In this work, we demonstrate for the first time, the influence of physiochemical cues (e.g., self-organized topography at nanoscale) of medical grade Ti₆Al₄V results in control of cell shape, adhesion, and proliferation of human aortic smooth muscle stem cells. The control of surface nanostructures was found to be correlated to ion-beam incidence angle linked to a surface diffusive regime during irradiation synthesis with argon ions at energies below 1 keV and a fluence of 2.5 × 10¹⁷ cm^-2. Cell viability and cytoskeleton morphology were evaluated at 24 h, observing an advance cell attachment state on post-DIS surfaces. These modified surfaces showed 84% of cell biocompatibility and an increase in cytoplasmatic protusions ensuring a higher cell adhesion state. Filopodia density was promoted by a 3-fold change for oblique incidence angle DIS treatment compared to controls (e.g., no patterning) and lamellipodia structures were increased more than a factor of 2, which are indicators of cell attachment stimulation due to DIS modification. In addition, the morphology of the nanofeatures were tailored, with high fidelity control of the main DIS parameters that control diffusive and erosive regimes of self-organization. We have correlated the morphology and the influence in cell behavior, where nanoripple formation is the most active morphology for cell stimulation.

Fragmentation approach to the point-island model with hindered aggregation: Accessing the barrier energy

We study the effect of hindered aggregation on the island formation process in a one- (1D) and two-dimensional (2D) point-island model for epitaxial growth with arbitrary critical nucleus size i. In our model, the attachment of monomers to preexisting islands is hindered by an additional attachment barrier, characterized by length l_{a}. For l_{a}=0 the islands behave as perfect sinks while for l_{a}→∞ they behave as reflecting boundaries. For intermediate values of l_{a}, the system exhibits a crossover between two different kinds of processes, diffusion-limited aggregation and attachment-limited aggregation. We calculate the growth exponents of the density of islands and monomers for the low coverage and aggregation regimes. The capture-zone (CZ) distributions are also calculated for different values of i and l_{a}. In order to obtain a good spatial description of the nucleation process, we propose a fragmentation model, which is based on an approximate description of nucleation inside of the gaps for 1D and the CZs for 2D. In both cases, the nucleation is described by using two different physically rooted probabilities, which are related with the microscopic parameters of the model (i and l_{a}). We test our analytical model with extensive numerical simulations and previously established results. The proposed model describes excellently the statistical behavior of the system for arbitrary values of l_{a} and i=1, 2, and 3.

A kinetic Monte Carlo simulation method of van der Waals epitaxy for atomistic nucleation-growth processes of transition metal dichalcogenides

Controlled growth of crystalline solids is critical for device applications, and atomistic modeling methods have been developed for bulk crystalline solids. Kinetic Monte Carlo (KMC) simulation method provides detailed atomic scale processes during a solid growth over realistic time scales, but its application to the growth modeling of van der Waals (vdW) heterostructures has not yet been developed. Specifically, the growth of single-layered transition metal dichalcogenides (TMDs) is currently facing tremendous challenges, and a detailed understanding based on KMC simulations would provide critical guidance to enable controlled growth of vdW heterostructures. In this work, a KMC simulation method is developed for the growth modeling on the vdW epitaxy of TMDs. The KMC method has introduced full material parameters for TMDs in bottom-up synthesis: metal and chalcogen adsorption/desorption/diffusion on substrate and grown TMD surface, TMD stacking sequence, chalcogen/metal ratio, flake edge diffusion and vacancy diffusion. The KMC processes result in multiple kinetic behaviors associated with various growth behaviors observed in experiments. Different phenomena observed during vdW epitaxy process are analysed in terms of complex competitions among multiple kinetic processes. The KMC method is used in the investigation and prediction of growth mechanisms, which provide qualitative suggestions to guide experimental study.

Recent developments in the synthesis of nanostructured chalcopyrite materials and their applications: a review

Nanostructured chalcopyrites: synthesis and applications.

The Ehrlich-Schwoebel barrier on an oxide surface: a combined Monte-Carlo and in situ scanning tunneling microscopy approach

The controlled growth of epitaxial films of complex oxides requires an atomistic understanding of key parameters determining final film morphology, such as termination dependence on adatom diffusion, and height of the Ehrlich-Schwoebel (ES) barrier. Here, through an in situ scanning tunneling microscopy study of mixed-terminated La5/8Ca3/8MnO3 (LCMO) films, we image adatoms and observe pile-up at island edges. Image analysis allows determination of the population of adatoms at the edge of islands and fractions on A-site and B-site terminations. A simple Monte-Carlo model, simulating the random walk of adatoms on a sinusoidal potential landscape using Boltzmann statistics is used to reproduce the experimental data, and provides an estimate of the ES barrier as ∼0.18 ± 0.04 eV at T = 1023 K, similar to those of metal adatoms on metallic surfaces. These studies highlight the utility of in situ imaging, in combination with basic Monte-Carlo methods, in elucidating the factors which control the final film growth in complex oxides.

Directed Kinetic Self-Assembly of Mounds on Patterned GaAs (001): Tunable Arrangement, Pattern Amplification and Self-Limiting Growth

We present results demonstrating directed self-assembly of nanometer-scale mounds during molecular beam epitaxial growth on patterned GaAs (001) surfaces. The mound arrangement is tunable via the growth temperature, with an inverse spacing or spatial frequency which can exceed that of the features of the template. We find that the range of film thickness over which particular mound arrangements persist is finite, due to an evolution of the shape of the mounds which causes their growth to self-limit. A difference in the film thickness at which mounds at different sites self-limit provides a means by which different arrangements can be produced.

Impact of nucleation on step-meandering instabilities during step-flow growth on vicinal surfaces

Step-meandering instabilities can manifest during step-flow growth on vicinal surfaces [Bales and Zangwill, Phys. Rev. B 41, 5500 (1990); Pierre-Louis, D'Orsogna, and Einstein, Phys. Rev. Lett. 82, 3661 (1999)]. A phase diagram based on the various growth regimes of a vicinal surface allows us to study the impact of nucleation on these meanders and to predict a meandering instability caused by the nucleation and the coalescence of both islands and steps. Using an accelerated kinetic Monte Carlo method, we find that the coalescence of islands with steps produces large protrusions and deep ripples and that the resulting meandering instability is reinforced by the growth of the islands at almost the same positions from one monolayer to the other. A coarsening phenomenon occurs for the instability wavelength until mounds appear, favored by a large Ehrlich-Schwoebel barrier. Such a meandering instability could be exploited for periodic self-assembly.

Coarsening dynamics in one dimension: the phase diffusion equation and its numerical implementation

Many nonlinear partial differential equations (PDEs) display a coarsening dynamics, i.e., an emerging pattern whose typical length scale L increases with time. The so-called coarsening exponent n characterizes the time dependence of the scale of the pattern, L(t)≈t(n), and coarsening dynamics can be described by a diffusion equation for the phase of the pattern. By means of a multiscale analysis we are able to find the analytical expression of such diffusion equations. Here, we propose a recipe to implement numerically the determination of D(λ), the phase diffusion coefficient, as a function of the wavelength λ of the base steady state u(0)(x). D carries all information about coarsening dynamics and, through the relation |D(L)|=/~L(2)/t, it allows us to determine the coarsening exponent. The main conceptual message is that the coarsening exponent is determined without solving a time-dependent equation, but only by inspecting the periodic steady-state solutions. This provides a much faster strategy than a orward time-dependent calculation. We discuss our method for several different PDEs, both conserved and not conserved.

Nonlinear wavelength selection in surface faceting under electromigration

We report on the control of the faceting of crystal surfaces by means of surface electromigration. When electromigration reinforces the faceting instability, we find perpetual coarsening with a wavelength increasing as t(1/2). For strongly stabilizing electromigration, the surface is stable. For weakly stabilizing electromigration, a cellular pattern is obtained, with a nonlinearly selected wavelength. The selection mechanism is not caused by an instability of steady states, as suggested by previous works in the literature. Instead, the dynamics is found to exhibit coarsening before reaching a continuous family of stable nonequilibrium steady states.

Geometric evolution law for modeling strongly anisotropic thin-film morphology

The morphology of the solid-vapor interface of a nanoscale thin crystalline film is influenced by many factors including surface diffusion, attachment-detachment, deposition, and interface kinetics. Using a high-order accurate and efficient numerical method, we investigate the dynamics of two dimensional thin films when all of these effects are considered. The observed morphologies consist of facets of constant slope separated by narrow transition intervals: kinks (valleys) and antikinks (hills). The number of kinks and antikinks decreases as the system coarsens in time. Our numerical results confirm that when deposition is present, the only possible coarsening event is the kink-ternary where two kinks meet and annihilate an antikink. We characterize the total amount of coarsening, the time over which the coarsening occurs and the associated coarsening scaling laws when all effects are considered. As found in previous work that considered only attachment-detachment, or surface diffusion, there are three distinct coarsening regimes associated with increasing magnitudes of the deposition flux-fast coarsening, a regime in which periodic structures form with little or no subsequent coarsening, and a regime in which the film surface evolves chaotically. We find that the inclusion of attachment-detachment leads to additional coarsening compared to the dynamics that result from driven surface diffusion alone. When deposition and interface kinetics are both considered, the slowdown of evolution caused by the kinetic effects necessitates a decrease in the deposition flux in order to produce a nonchaotic coarsening regime. Together, these provide testable predictions for experiments of thin-film dynamics.

Continuum imit of a tep Flow odel of Epitaxial Growth

We examine a class of step ow models of epitaxial growth obtained from a Burton-Cabrera-Frank (BCF) type approach in one space dimension. ur goal is to derive a consistent contin uummodel for the ev olutionof the lm surface. Away from peaks and valleys, the surface height solves a Hamilton- acobi equation (H E). he peaks are free boundaries for this H E. heir evolution must be speci ed by boundary conditions re ecting the microscopic physics of nucleation. e investigate this boundary condition by numerical simulation of the step ow dynamics using a simple n ucleationlaw. ur results rev ealthe presence of sp ecial structures in the pro le near a peak; we discuss the relationship between these structures and the contin uumequation. e further address the importance of ev aporationfor matching the local behavior near the peak to the solution of the contin uum equation.

Growth instability due to lattice-induced topological currents in limited-mobility epitaxial growth models

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.

Demixing processes in AgPd superlattices

The present scanning tunneling microscopy (STM) study describes the growth of silver-palladium heterostructures at room temperature, with ab initio simulations of ordered AgPd phases supporting the interpretation of STM images. First, the growth of Pd on an Ag(111) surface proceeds in a multilayer mode, leading to the formation of a columnar structure. Then, upon Ag deposition on this structure, Ag and Pd partially mix and form a two-dimensional AgPd alloy on top of the columns. Finally, an atomically flat Ag(111) surface is restored, and two-dimensional growth continues. An interpretation of this peculiar growth mode including interfacial alloying is proposed based on thermodynamic and kinetic arguments.

Formation of complex wedding-cake morphologies during homoepitaxial film growth of Ag on Ag(111): atomistic, step-dynamics, and continuum modeling

An atomistic lattice-gas model is developed which successfully describes all key features of the complex mounded morphologies which develop during deposition of Ag films on Ag(111) surfaces. We focus on this homoepitaxial thin film growth process below 200 K. The unstable multilayer growth mode derives from the presence of a large Ehrlich-Schwoebel step-edge barrier, for which we characterize both the step-orientation dependence and the magnitude. Step-dynamics modeling is applied to further characterize and elucidate the evolution of the vertical profiles of these wedding-cake-like mounds. Suitable coarse-graining of these step-dynamics equations leads to instructive continuum formulations for mound evolution.

Scale decomposition of molecular beam epitaxy

In this work, a study of epitaxial growth was carried out by means of the wavelets formalism. We showed the existence of a dynamic scaling form in a wavelet discriminated linear molecular beam epitaxy (MBE) equation where diffusion and noise are the dominant effects. We determined simple and exact scaling functions involving the scale of the wavelets when the system size is set to infinity. Exponents were determined for both correlated and uncorrelated noise. The wavelet methodology was applied to a computer model simulating linear epitaxial growth; the results showed very good agreement with analytical formulation. We also considered epitaxial growth with the additional Ehrlich-Schwoebel effect. We characterized the coarsening of mounds formed on the surface during the nonlinear phase using the wavelet power spectrum. The latter has an advantage over other methods, in the sense that one can track the coarsening in both frequency (or scale) space and real space simultaneously. Wavelets analysis also provides a quantitative tool for the characterization of the mounded surfaces through its concise scale discrimination. We showed that the averaged wavelet power spectrum (also called scalegram) over all the positions on the surface profile identified the existence of a dominant scale a(*), which increases with time following a power law relation of the form a(*)∼t(n), where [Formula: see text].

Theoretical analysis of mound slope selection during unstable multilayer growth

A "step dynamics" model is developed for mound formation during multilayer homoepitaxy. Downward funneling of atoms deposited at step edges is incorporated and controls mound slope selection. Behavior of the selected slope differs from that predicted by phenomenological continuum treatments where the lateral mass current vanishes identically. Instead, this current is shown to vary periodically and vanish only on average. An exact coarse-grained continuum formulation with appropriate boundary conditions is derived and recovers step dynamics results.

Kinetics of step bunching during growth: a minimal model

We study a minimal stochastic model of step bunching during growth on a one-dimensional vicinal surface. The formation of bunches is controlled by the preferential attachment of atoms to descending steps (inverse Ehrlich-Schwoebel effect) and the ratio d of the attachment rate to the terrace diffusion coefficient. For generic parameters (d>0) the model exhibits a very slow crossover to a nontrivial asymptotic coarsening exponent beta approximately 0.38 . In the limit of infinitely fast terrace diffusion (d=0) linear coarsening (beta=1) is observed instead. The different coarsening behaviors are related to the fact that bunches attain a finite speed in the limit of large size when d=0 , whereas the speed vanishes with increasing size when d>0 . For d=0 an analytic description of the speed and profile of stationary bunches is developed, and a connection to the problem of front propagation into an unstable state is pointed out.

Interplay of instabilities in mounded surface growth

We numerically study a one-dimensional conserved growth equation with competing linear (Ehrlich-Schwoebel) and nonlinear instabilities. As a control parameter is varied, this model exhibits a nonequilibrium phase transition between two mounded states, one of which exhibits slope selection and the other does not. The coarsening behavior of the mounds in these two phases is studied in detail. In the absence of noise, the steady-state configuration depends crucially on which of the two instabilities dominates the early time behavior.

Kinetic roughening of ion-sputtered Pd(001) surface: beyond the Kuramoto-Sivashinsky model

We investigate the kinetic roughening of Ar+ ion-sputtered Pd(001) surface both experimentally and theoretically. In situ real-time x-ray reflectivity and in situ scanning tunneling microscopy show that nanoscale adatom islands form and grow with increasing sputter time t. Surface roughness W(t) and lateral correlation length xi(t) follow the scaling laws W(t) approximately t(beta) and xi(t) approximately t(1/z) with the exponents beta approximately 0.20 and 1/z approximately 0.20, for an ion beam energy epsilon=0.5 keV, which is inconsistent with the prediction of the Kuramoto-Sivashinsky (KS) model. We thereby extend the KS model by applying the coarse-grained continuum approach of the Sigmund theory to the order of O(inverted Delta(4),h(2)), where h is the surface height, and derive a new term of the form inverted Delta(2)(inverted Delta h)(2) which plays a decisive role in describing the observed morphological evolution of the sputtered surface.

When does coarsening occur in the dynamics of one-dimensional fronts?

Dynamics of a one-dimensional growing front with an unstable straight profile are analyzed. We argue that a coarsening process occurs if and only if the period lambda of the steady-state solution is an increasing function of its amplitude A. This statement is rigorously proved for two important classes of conserved and nonconserved models by investigating the phase diffusion equation of the steady pattern. We further provide clear numerical evidence for the growth equation of a stepped crystal surface.

Mound formation and coarsening from a nonlinear instability in surface growth

We study spatially discretized versions of a class of one-dimensional, nonequilibrium, conserved growth equations for both nonconserved and conserved noise using numerical integration. An atomistic version of these growth equations is also studied using stochastic simulation. The models with nonconserved noise are found to exhibit mound formation and power-law coarsening with slope selection for a range of values of the model parameters. Unlike previously proposed models of mound formation, the Ehrlich-Schwoebel step-edge barrier, usually modeled as a linear instability in growth equations, is absent in our models. Mound formation in our models occurs due to a nonlinear instability in which the height (depth) of spontaneously generated pillars (grooves) increases rapidly if the initial height (depth) is sufficiently large. When this instability is controlled by the introduction of a nonlinear control function, the system exhibits a first-order dynamical phase transition from a rough self-affine phase to a mounded one as the value of the parameter that measures the effectiveness of control is decreased. We define an "order parameter" that may be used to distinguish between these two phases. In the mounded phase, the system exhibits power-law coarsening of the mounds in which a selected slope is retained at all times. The coarsening exponents for the spatially discretized continuum equation and the atomistic model are found to be different. An explanation of this difference is proposed and verified by simulations. In the spatially discretized growth equation with conserved noise, we find the curious result that the kinetically rough and mounded phases are both locally stable in a region of parameter space. In this region, the initial configuration of the system determines its steady-state behavior.

Surface kinetics and generation of different terms in a conservative growth equation

A method based on the kinetics of adatoms on a growing surface under epitaxial growth at low temperature in (1+1) dimensions is proposed to obtain a closed form of the local growth equation. It can be generalized to any growth problem where surface morphology is governed by adatom diffusion. The method can be easily extended to higher dimensions. The kinetic processes contributing to various terms in the growth equation are identified from the analysis of in-plane and downward hops. In particular, processes corresponding to the term that breaks h-->-h symmetry and the curvature dependent term are discussed. Effects of these terms on the stable to unstable transition in (1+1) dimension are analyzed. In (2+1) dimensions, it is shown that an additional asymmetric term is generated due to the in-plane curvature associated with mound-like structures. This term is independent of any diffusion barrier differences between in-plane and out-of-plane migration. It is shown that terms generated in the presence of downward hops are the relevant terms in a growth equation. A growth equation in closed form is obtained for various growth models introduced to capture most of the processes in experimental molecular beam epitaxial growth. The effect of dissociation is also considered and is seen to have a stabilizing effect on growth. It is shown that for uphill current the growth equation approach fails to describe the growth since a given single equation does not apply over the entire substrate.

Mounding instability and incoherent surface kinetics

Mounding instability in a conserved growth from vapor is analyzed within the framework of adatom kinetics on the growing surface. The analysis shows that depending on the local structure of the surface, kinetics of adatoms may vary, leading to disjoint regions in the sense of a continuum description. This is manifested particularly under the conditions of instability. Mounds grow on these disjoint regions and their lateral growth is governed by the flux of adatoms hopping across the steps in the downward direction. Asymptotically ln t dependence is expected in (1+1) dimensions. Simulation results confirm the prediction. Growth in (2+1) dimensions is also discussed.

Layer-by-layer epitaxy in limited mobility nonequilibrium models of surface growth

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.

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.

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: on 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.

Oscillatory driving of crystal surfaces: a route to controlled pattern formation

We show that the oscillatory driving of crystal surfaces can induce pattern formation or smoothening. Depending on driving conditions, step bunching and meandering, mound formation, or surface smoothening may be seen in the presence of a kinetic asymmetry at the steps or kinks. We employ a step model to calculate the induced mass flux along misoriented surfaces, which accounts for surface dynamics and stability. Slope selection, surface metastability, and frequency-dependent surface stability are found. Quantitative predictions for pattern formation on metal surfaces in an electrolyte are provided.

Instability and wavelength selection during step flow growth of metal surfaces vicinal to fcc(001)

We study the onset and development of ledge instabilities during growth of vicinal metal surfaces using kinetic Monte Carlo simulations. We observe the formation of periodic patterns at [110] close packed step edges on surfaces vicinal to fcc(001) under realistic molecular beam epitaxy conditions. The corresponding wavelength and its temperature dependence are studied in detail. Simulations suggest that the ledge instability on fcc(1,1,m) vicinal surfaces is controlled by the strong kink Ehrlich-Schwoebel barrier, with the wavelength determined by dimer nucleation at the step edge. Our results are in agreement with recent continuum theoretical predictions, and experiments on Cu(1,1,17) vicinal surfaces.

(2+1)-dimensional stochastic growth model and its application to some experimental observations

A (2+1)-dimensional stochastic growth model is proposed in order to understand various experimental observations at low temperature growth. The model includes step attached adatom hopping in the upward direction to provide a source of additional particle current. The flexibility obtained due to additional current control allows one to express most of the experimental observations within the framework of this model. It is argued that the time evolution exponent beta for surface roughness is >0.5 only if upward hops are included. Simulations using this model are shown to be qualitatively consistent with the observations for homoepitaxial growth on Cu(100) and Ge(100) at different temperatures.