An efficient program with a user-friendly graphical interface has been developed for simulating reflection high-energy electron diffraction patterns obtained on crystal surfaces. The calculations are based on a kinematical approach which considers single electron scattering events at the surface. This time-efficient approach is in most cases sensitive enough for distinguishing different structural models, even if they differ very subtly. The results are presented in realistic two-dimensional density plots which can be directly compared to experimental observations. This program provides a useful tool in studying different structural models for crystal surfaces. Program summaryProgram title: RHEEDsimCatalogue identifier: AEJC_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEJC_v1_0.htmlProgram obtainable from: CPC Program Library, Queenʼs University, Belfast, N. IrelandLicensing provisions: Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.htmlNo. of lines in distributed program, including test data, etc.: 2752No. of bytes in distributed program, including test data, etc.: 367 205Distribution format: tar.gzProgramming language: Matlab (version > 7.6.0)Computer: Personal ComputersOperating system: Windows with Matlab environmentRAM: Greater than 1 MBClassification: 7.4Nature of problem: Reflection high-energy electron diffraction (RHEED) has been widely used in studying crystal surface structures all over the world, especially in combination with ultra-high vacuum molecular beam epitaxy systems. In addition to determining the surface smoothness, RHEED is also a very useful tool in studying surface reconstructions, which are often encountered at the growth surfaces of semiconductors and alloys. While the positions of the fractional streaks can be used to determine the basic information of the surface supercells, the intensity modulations on the fractional streaks provide further insights on the details within these unit cells. Kinematic approach is an efficient method for simulating the RHEED patterns based on various surface structural models, which can help unraveling the surface atomic structure.Solution method: The kinematic approach utilized here assumes single scattering events. Furthermore, the Ewald sphere is approximated into a planar surface for computing the streak intensities (which are most relevant to real experiments). Structure factors are calculated based on a given input of atomic species and their coordinates, with user modifiable form factors. In addition to the intensity modulation within the surface plane, additional modulations extending into the z direction is also taken into consideration, resulting in realistic density maps of the RHEED streaks, which can be directly compared to experimental observations.Unusual features: The main program RHEEDsim.m calls several local subprograms for certain computational tasks. As a result, all programs should be extracted into a single folder, and that folder should be set as the main directory in Matlab.Running time: The computing time is computer and user parameter dependent, but typically ranges from few seconds to few minutes.