Elsevier

Computer Physics Communications

Volume 207, October 2016, Pages 536-538
Computer Physics Communications

RHEED intensities from two-dimensional heteroepitaxial nanoscale systems of GaN on a Si surface

https://doi.org/10.1016/j.cpc.2016.05.029Get rights and content

Abstract

This paper presents a computer program, which facilitates the calculation of changes to the intensity of RHEED oscillations from the heteroepitaxial structures of (0001)GaN films nucleated on a Si surface. The calculations are based on the use of a dynamical diffraction theory and different models of scattering crystal potential.

New version program summary

Title of program: RHEED_DIFF_W

Program Files doi:10.17632/hn3pt6ytky.1

Licensing provisions: GNU General Public License 3

Programming language used: C++

Journal reference of previous version: Computer Physics Communications 185 (2014) 3001–3009

Does the new version supersede the previous version?: No. It does supersede version AETW_v1_0.

Reasons for the new version: Responding to the users’ feedback we presented a practical procedure of construction of simulation program, which facilitates the calculation of changes to the intensity of RHEED oscillations in the function of the glancing angle of incidence of the electron beam, employing various models of scattering crystal potential for heteroepitaxial structures of hexagonal (0001)GaN films nucleated on a Si surface. GaN has a wurtzite crystal structure with alternating planes of Ga and N atoms and lattice constants a=3.189Å and c=5.185Å [1]. Choosing the top surface to be z=0, the N atoms are located at z=nd where d=5.185/2Å and n is an integer, and Ga atoms are at z=ndu, where u=3c/8 (as shown in Fig. 1).

The nature of physical problem: Since the discovery of molecular beam epitaxy, which allows formation of layers (crystals) with a controlled composition of subsequent monolayers, the research attempts of many laboratories have been focused on studying the dynamics of formation of both monolayers and specific heterojunction structures, that is the structures with monocrystalline interfaces created by materials of various chemical compositions (different physical and chemical properties). The fundamental scientific problem in such research is to specify both interface type and growth mechanism for subsequent layers. Researchers and technologists manufacturing two-dimensional nanoscale systems frequently use RHEED rocking curves (the specular beam intensities versus the glancing angle) to control growth of films at the atomistic level of accuracy. The structural properties of the GaN/Si interface have been meticulously studied by many [3–5], and crucial problem of the direct growth of GaN on Si is the formation of amorphous silicon nitride (SiN) at the GaN/Si interface.

Typical running time: The typical running time is machine and user-parameters dependent, i.e. the number of layers for both the substrate and the growing layers included in the calculations.

The method of solution: RHEED intensities are calculated within the general framework described in Ref. [2] with the model of the scattering potential for heterostructures: (1)UCombined(z)=nUsubstrate(z)+n(Ulayer(z)+Uaddlayer(z)), where the Uaddlayer represents the diffuse scattering on the topmost layer.

Summary of revisions: The presented version of computer program implements original algorithms of self-consistent calculations for scattering potentials GaN(0001)/Si(111) and solving a time-independent Schrödinger equation for high-energy electrons (Fig. 2).

  1. Download : Download high-res image (238KB)
  2. Download : Download full-size image

Fig. 1. A plot of the real part of the one-dimensional potential of GaN(0001) for 12 layers thickness.

  1. Download : Download high-res image (134KB)
  2. Download : Download full-size image

Fig. 2. The specular beam intensities versus the glancing angle for the GaN(0001) corresponds to the data shown in Fig. 1. Electron energy equals 10 keV.

  1. Download : Download high-res image (175KB)
  2. Download : Download full-size image

Fig. 3. The one-dimensional (real part) potential of GaN(0001)/Si(111) for 10 layer thickness of Si and one of GaN corresponds to the data shown in Fig. 4.

  1. Download : Download high-res image (144KB)
  2. Download : Download full-size image

Fig. 4. A set of rocking curves calculated during the growth of a single GaN layer on the Si(111) substrate for two different values θGaN of the coverage of the topmost layer. For the calculations the MODEL3 [2] of the scattering potential was used.

During the numerical calculations of the changes in the intensity of the specular beam in the function of the glancing angle, it was assumed that the electron energy equals 10 keV, the glancing angle was increased from 0.5° to 7°, and the value of α parameter =0.1 and β=0.5 for the MODEL3 of the scattering potential [2]. Figs. 34 present the situation where on the crystalline Si(111) substrate a single GaN layer was placed, where the GaN growth on Si is having a thin amorphous layer. It can be observed that for a very thin layer of GaN, a clear minimal amplitude of RHEED oscillation occurs with the coverage of the GaN layer equal to about 0.5.

References

  • [1]

    A. Ichimiya, P.I. Cohen, Reflection High-Energy Electron Diffraction, Cambridge University Press, Cambridge, 2004.

  • [2]

    A. Daniluk, Comput. Phys. Commun. 185 (2014) 3001.

  • [3]

    M. Tamura, M. López-López, Superficies y Vacío 13 (2001) 80.

  • [4]

    J. Liang, S.-K. Hong, N. Kouklin, R. Beresford, and J.M. Xu, Appl. Phys. Lett. 83(9) (2003) 1752.

  • [5]

    M. Sobanska, K. Klosek, Z.R. Zytkiewicz, J. Borysiuk, B.S. Witkowski, E. Lusakowska, A. Reszka, and R. Jakiela, Cryst. Res. Technol. 47(3) (2012) 307.

References (0)

Cited by (0)

View full text