A Python code for simulations of RHEED intensity oscillations within the one-dimensional dynamical approximation

We present a Python-based implementation of 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 crystal potential for heteroepitaxial structures including the possible existence of various diffuse scattering models through the layer parallel to the surface. The calculations are based on the use of a one-dimensional dynamical diffraction theory. Although this theory has some limitations, in practice it is useful under so-called one-beam condition. Computation performance has been improved by using Numba as an open source, NumPy-aware optimising compiler for Python.

The previous version of this program (AETW_v1_0) may be found at https://doi.org/10.1016/j.cpc.2014.07.003.