Surface growth

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:

${\displaystyle A_{in}=A_{0,in}\exp \left(-{\frac {E_{a,in}}{kT}}\right)}$

where A is thermal frequency of molecular vibration, k is Boltzmann constant.

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 described exactly as in the absorption process, with the exception of a different activation energy barrier.

${\displaystyle A_{out}=A_{0,out}\exp \left(-{\frac {E_{a,out}}{kT}}\right)}$

For example, if all positions on the surface of the crystal are energy equivalent, the rate of growth can be calculated from Turnbull formula:

${\displaystyle V_{c}=hC_{0}(A_{out}-A_{0,out})=hC_{0}\exp \left(-{\frac {E_{a,in}}{kT}}\right)\cdot \left(1-\exp \left(-{\frac {\Delta G}{kT}}\right)\right)}$

where, ∆G = E_in – E_out, A_out, A_{o 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, C_o concentration of the molecules in direct distance from the surface.

Diffusion process can also be calculated with Arrhenius equation:

${\displaystyle D=D_{0}\exp \left(-{\frac {E_{d}}{kT}}\right)}$

where, D is diffusion coefficient, E_d 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.

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.

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.

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.

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.

Conclusion: low impingement rate helps creating smoother surface.

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:

• Das Sarma, S.; Tamborenea, P. (21 January 1991). "A new universality class for kinetic growth: One-dimensional molecular-beam epitaxy". Physical Review Letters. American Physical Society (APS). 66 (3): 325–328. doi:10.1103/physrevlett.66.325. ISSN 0031-9007. PMID 10043777.
• Levi, Andrea C; Kotrla, Miroslav (13 January 1997). "Theory and simulation of crystal growth". Journal of Physics: Condensed Matter. IOP Publishing. 9 (2): 299–344. doi:10.1088/0953-8984/9/2/001. ISSN 0953-8984.
• Meng, B.; Weinberg, W.H. (1996). "Dynamical Monte Carlo studies of molecular beam epitaxial growth models: interfacial scaling and morphology". Surface Science. Elsevier BV. 364 (2): 151–163. doi:10.1016/0039-6028(96)00597-3. ISSN 0039-6028.
• Wadley, H.N.G; Zhou, X; Johnson, R.A; Neurock, M (2001). "Mechanisms, models and methods of vapor deposition". Progress in Materials Science. Elsevier BV. 46 (3–4): 329–377. doi:10.1016/s0079-6425(00)00009-8. ISSN 0079-6425.
• Wolf, D. E; Villain, J (1 October 1990). "Growth with Surface Diffusion". Europhysics Letters (EPL). IOP Publishing. 13 (5): 389–394. doi:10.1209/0295-5075/13/5/002. ISSN 0295-5075.
• Xiao, Rong-Fu; Alexander, J. Iwan D.; Rosenberger, Franz (1 February 1991). "Growth morphologies of crystal surfaces". Physical Review A. American Physical Society (APS). 43 (6): 2977–2992. doi:10.1103/physreva.43.2977. ISSN 1050-2947.
• Lars Röntzsch. "Vicinal surface diffusion". Retrieved 23 May 2019.