Kinetic Monte Carlo surface growth method
![]() | This article has multiple issues. Please help improve it or discuss these issues on the talk page. (Learn how and when to remove these template messages)
|
Kinetic Monte Carlo (KMC) is a form of computer simulation in which atoms and molecules are allowed to interact at given rate that could be controlled based on known physics. This simulation method is typically used in the micro-electrical industry to study crystal surface growth, and it can provide accurate models surface morphology in different growth conditions on a time scales typically ranging from micro-seconds to hours. Experimental methods such as Scanning Electron Microscopy (SEM), X-ray diffraction, and Transmission Electron Microscopy (TEM), and other computer simulation methods such as Molecular Dynamics (MD), and Monte Carlo (MC) simulation are widely used.
How KMC surface growth works
1. Absorption process
First, the model tries to predict where an atom would land on a surface and its rate at particular environmental conditions, such as temperature and vapor pressure. In order to land on a surface, atoms have to overcome the so-called activation energy barrier. The frequency of passing through the activation barrier can by calculated by the Arrhenius equation:
Ain = Ao inexp (-Ea in/kT)
where A is thermal frequency of molecule vibration, k is Boltzmann constant.
2. Desorption process
When atoms land on a surface, there are two possibilities. First, they would diffuse on the surface and find other atoms to make a cluster, which will be discussed below. Second, they could come off of the surface or so-called desorption process. The desorption is describe exactly as in the absorption process, with the exception of a different activation energy barrier.
Aout = Ao outexp (-Ea out/kT)
For example, if all positions on the surface of the crystal are energy equivalent, the rate of growth can be calculated from Turnbull formula:
Vc = h*Co (Aout - Ao out) = h*Coexp (-Ea in/kT) * (1 - exp(-∆G/kT))
where, ∆G = Ein – Eout, Aout, Ao out are frequencies to go in or out of crystal for any given molecule on the surface, h - height of the molecule in the growth direction, Co concentration of the molecules in direct distance from the surface.
3. Diffusion process on surface
Diffusion process can also be calculated with Arrhenius equation:
D = Do exp(-Ed/kT)
where, D is diffusion coefficient, Ed is diffusion activation energy.
All three processes strongly depend on surface morphology at a certain time. For example, atoms tend to lend at the edges of a group of connected atoms, the so-called island, rather than on a flat surface, this reduces the total energy. When atoms diffuse and connect to an island, each atom tends to diffuse no further, because activation energy to detach itself out of the island is much higher. Moreover, if an atom landed on top of an island, it would not diffuse fast enough, and the atom would tend to move down the steps and enlarge it.
Simulation Methods
Because of limited computing power, specialized simulation models have been developed for various purposes depending on the time scale:
a) Electronic scale simulations (Density Function Theory, Ab-initio Molecular Dynamics): sub-atomic length scale in femto-second time scale
b) Atomic scale simulations (MD): nano to micro-meter length scale in nano-second time scale
c) Film scale simulation (KMC): micro-meter length scale in micro to hour time scale.
d) Reactor scale simulation (phase field model): meter length scale in year time scale.
multiscale modeling techniques have also been developed to deal with overlapping time scales.
How to use growth conditions in KMC
The interest of growing a smooth and defect-free surface requires a combination set of physical conditions throughout the process. Such conditions are bond strength, temperature, surface-diffusion limited and supersaturation (or impingement) rate. Using KMC surface growth method, following pictures describe final surface structure at different conditions.
1. Bond strength and temperature
Bond strength and temperature certainly play important roles in the crystal grow process. For high bond strength, when atoms land on a surface, they tend to be closed to atomic surface clusters, which reduce total energy. This behavior results in many isolated cluster formations with a variety of size yielding a rough surface. Temperature, on the other hand, controls the high of the energy barrier.
Conclusion: high bond strength and low temperature is preferred to grow a smoothed surface.
2. Surface and bulk diffusion effect
Thermodynamically, a smooth surface is the lowest ever configuration, which has the smallest surface area. However, it requires a kinetic process such as surface and bulk diffusion to create a perfectly flat surface.
Conclusion: enhancing surface and bulk diffusion will help create a smoother surface.
3. Supersaturation level
Conclusion: low impingement rate helps creating smoother surface.
4. Morphology at different combination of conditions
With the control of all growth conditions such as temperature, bond strength, diffusion, and saturation level, desired morphology could be formed by choosing the right parameters. Following is the demonstration how to obtain some interesting surface features:
References
- S. Das Sarma and P. Tamborenea, "A new universality class for kinetic growth: one-dimensional molecular-beam epitaxy", Phys. Rev. Lett. 66 (1991), 325
- Levi A.C. and Kotrla M., J. Phys.: Condens. Matter 9 (1997), 299-344
- Meng B. and Weinberg W.H., Surface Science 364 (1996), 151-163.
- Wadley H.N.C., Zhou X., Johnson R.A., and Neurock M., Progress in Material Science 46 (2001), 329-377
- Wolf D.E. and Villain J., "Growth with surface diffusion", Europhys. Lett. 13 (1990), 389-394
- Xiao R.F., Alexander J.I.D., and Rosenberger F., Phys Rev. A. 43 (1991), 2977
- http://www.roentzsch.org/SurfDiff/