Fast kinetic Monte Carlo simulation of strained heteroepitaxy in three dimensions

Accelerated algorithms for simulating the morphological evolution of strained heteroeptiaxy based on a ball and spring lattice model in three dimensions are explained. We derive exact Green's function formalisms for boundary values in the associated lattice elasticity problems. The computational efficiency is further enhanced by using a superparticle surface coarsening approximation. Atomic hoppings simulating surface diffusion are sampled using a multi-step acceptance-rejection algorithm. It utilizes quick estimates of the atomic elastic energies from extensively tabulated values modulated by the local strain. A parameter controls the compromise between accuracy and efficiency of the acceptance-rejection algorithm.Comment: 10 pages, 4 figures, submitted to Proceedings of Barrett Lectures 2007, Journal of Scientific Computin

Off-lattice Kinetic Monte Carlo simulations of strained heteroepitaxial growth

An off-lattice, continuous space Kinetic Monte Carlo (KMC) algorithm is discussed and applied in the investigation of strained heteroepitaxial crystal growth. As a starting point, we study a simplifying (1+1)-dimensional situation with inter-atomic interactions given by simple pair-potentials. The model exhibits the appearance of strain-induced misfit dislocations at a characteristic film thickness. In our simulations we observe a power law dependence of this critical thickness on the lattice misfit, which is in agreement with experimental results for semiconductor compounds. We furthermore investigate the emergence of strain induced multilayer islands or "Dots" upon an adsorbate wetting layer in the so-called Stranski-Krastanov (SK) growth mode. At a characteristic kinetic film thickness, a transition from monolayer to multilayer islands occurs. We discuss the microscopic causes of the SK-transition and its dependence on the model parameters, i.e. lattice misfit, growth rate, and substrate temperature.Comment: 17 pages, 6 figures Invited talk presented at the MFO Workshop "Multiscale modeling in epitaxial growth" (Oberwolfach, Jan. 2004). Proceedings to be published in "International Series in Numerical Mathematics" (Birkhaeuser

Computation of Strained Epitaxial Growth in Three Dimensions by Kinetic Monte Carlo

A numerical method for computation of heteroepitaxial growth in the presence of strain is presented. The model used is based on a solid-on-solid model with a cubic lattice. Elastic effects are incorporated using a ball and spring type model. The growing film is evolved using Kinetic Monte Carlo (KMC) and it is assumed that the film is in mechanical equilibrium. The strain field in the substrate is computed by an exact solution which is efficiently evaluated using the fast Fourier transform. The strain field in the growing film is computed directly. The resulting coupled system is solved iteratively using the conjugate gradient method. Finally we introduce various approximations in the implementation of KMC to improve the computation speed. Numerical results show that layer-by-layer growth is unstable if the misfit is large enough resulting in the formation of three dimensional islands

Kinetic Monte Carlo simulation of faceted islands in heteroepitaxy using multi-state lattice model

A solid-on-solid model is generalized to study the formation of Ge pyramid islands bounded by (105) facets on Si(100) substrates in two dimensions. Each atomic column is not only characterized by the local surface height but also by two deformation state variables dictating the local surface tilt and vertical extension. These deformations phenomenologically model surface reconstructions in (105) facets and enable the formation of islands which better resemble faceted pyramids. We demonstrate the model by application to a kinetic limited growth regime. We observe significantly reduced growth rates after faceting and a continuous nucleation of new islands until overcrowding occurs.Comment: 7 pages, 5 figure

Multiscale Kinetic Monte Carlo Simulation of Self-Organized Growth of GaN/AlN Quantum Dots

A three-dimensional kinetic Monte Carlo methodology is developed to study the strained epitaxial growth of wurtzite GaN/AlN quantum dots. It describes the kinetics of effective GaN adatoms on an hexagonal lattice. The elastic strain energy is evaluated by a purposely devised procedure: first, we take advantage of the fact that the deformation in a lattice-mismatched heterostructure is equivalent to that obtained by assuming that one of the regions of the system is subjected to a properly chosen uniform stress (Eshelby inclusion concept), and then the strain is obtained by applying the Green’s function method. The standard Monte Carlo method has been modified to implement a multiscale algorithm that allows the isolated adatoms to perform long diffusion jumps. With these state-of-the art modifications, it is possible to perform efficiently simulations over large areas and long elapsed times. We have taylored the model to the conditions of molecular beam epitaxy under N-rich conditions. The corresponding simulations reproduce the different stages of the Stranski–Krastanov transition, showing quantitative agreement with the experimental findings concerning the critical deposition, and island size and density. The influence of growth parameters, such as the relative fluxes of Ga and N and the substrate temperature, is also studied and found to be consistent with the experimental observations. In addition, the growth of stacked layers of quantum dots is also simulated and the conditions for their vertical alignment and homogenization are illustrated. In summary, the developed methodology allows one to reproduce the main features of the self-organized quantum dot growth and to understand the microscopic mechanisms at play

Competing roughening mechanisms in strained heteroepitaxy: a fast kinetic Monte Carlo study

We study the morphological evolution of strained heteroepitaxial films using kinetic Monte Carlo simulations in two dimensions. A novel Green's function approach, analogous to boundary integral methods, is used to calculate elastic energies efficiently. We observe island formation at low lattice misfit and high temperature that is consistent with the Asaro-Tiller-Grinfeld instability theory. At high misfit and low temperature, islands or pits form according to the nucleation theory of Tersoff and LeGoues.Comment: 4 pages, 4 figures, ReVTe

Kinetic Monte Carlo Simulation of Strained Heteroepitaxial Growth with Intermixing

An efficient method for the simulation of strained heteroepitaxial growth with intermixing using kinetic Monte Carlo is presented. The model used is based on a solid-on-solid bond counting formulation in which elastic effects are incorporated using a ball and spring model. While idealized, this model nevertheless captures many aspects of heteroepitaxial growth, including nucleation, surface diffusion, and long range effects due elastic interaction. The algorithm combines a fast evaluation of the elastic displacement field with an efficient implementation of a rejection-reduced kinetic Monte Carlo based on using upper bounds for the rates. The former is achieved by using a multigrid method for global updates of the displacement field and an expanding box method for local updates. The simulations show the importance of intermixing on the growth of a strained film. Further the method is used to simulate the growth of self-assembled stacked quantum dots