1The production of self-organized nanostructures by ion beam sputtering has been of keen interest to researchers for many decades. Despite numerous experimental and theoretical efforts to understand ion-induced nanostructures, there are still many basic questions open to discussion, such as the role of erosion or curvature-dependent sputtering. In this work, a hybrid MD/kMC (Molecular Dynamics/kinetic Monte Carlo) multiscale atomistic model is developed to investigate these knowledge gaps, and its predictive ability is validated across the experimental parameter space. This model uses crater functions, which were obtained from MD simulations, to model the prompt mass redistribution due to single-ion impacts. Defect migration, which is missing from previous models that use crater functions, is treated by a kMC Arrhenius method. Using this model, a systematic study was performed for silicon bombarded by Ar + ions of various energies (100 eV, 250 eV, 500 eV, 700 eV and Finally, simulations of dual-beam sputtering were performed, with the resulting patterns determined by the flux ratio of the two beams and the pattern types resulting from single-beam sputtering under the same conditions. 2