Due to the complexity of the binary component system and the difficulty of tracing individual particles in experiments, it is highly desirable to develop simulation tools and models to further reveal the microscopic nucleation behavior of Si and C atoms. In this article, self-consistent rate equations (SCRE) theory combined with the Kinetic Monte Carlo (KMC) model are used to study the nucleation mechanism in the early stage of SiC(0001) surface epitaxial growth under a constant deposition flux. A set of rate equations describing the time evolution of the density of monomers, dimers, and islands are established. By introducing the effective absorption length, the rate equations can be solved self-consistently. In the KMC model, a set of crystal lattices of SiC, including the coordinates of individual Si and C particles and the bond indication, are established. In this model, deposition, adatom diffusion, attachment to and detachment from the clusters of Si or C, and their diffusion along the edge of clusters are considered, and the Hoshen–Kopelman algorithm is implemented to identify and label the clusters. The results show that the time evolution of the density of monomers, dimers, and islands and their dependences on the deposition flux, which are obtained from the SCRE theory, are consistent with the results of the KMC model. The nucleation rate of the dimer increases with the increase in the deposition flux until the number of islands becomes saturated, which leads to a higher density of stable clusters at higher flux. Two regimes governing the absorption length of the monomer are presented. First, before the appearance of the nucleus, the absorption length of the monomer is mainly dominated by monomer density and the diffusion coefficient of adatom. Second, with the growth process, the total capture rate of stable islands gradually plays a dominant role, which is responsible for the further reduction in absorption length and the equal absorption length of Si and C.
In order to investigate the microscopic evolution of the step flow growth process and reveal the microscopic origins of crystalline anisotropy during the epitaxial growth of 3C-SiC (0001) vicinal surface, a three-dimensional Kinetic Monte Carlo model is established, in which Si and C are considered individually. The helicoidal boundary condition is applied to the direction perpendicular to the step, and the periodic boundary condition is used in the direction along the step. Then, the effects of crystalline anisotropy on lateral growth rate, morphologies of step patterns, and growth mode are studied. The results show that the lateral growth rate in [1−210] is larger than that in [Formula: see text], and the zigzag and meandering patterns of step are constructed in [1−210] and [Formula: see text] directions, respectively, which is consistent with the experimental observations. Two possible origins of anisotropy are also revealed: one is the higher concentration of the edge sites of the step and the larger bonding energy in the [1−210] direction and another is the adatom diffusion along the edge of the step. Finally, a larger area of pure step flow growth mode is obtained in the [1−210] direction, which is good for lowering the cost for 3C-SiC epitaxial layers.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.