Molecular Dynamics Calculation for Low-Energy Ion Implantation Process with Dynamic Annealing Effect

Abstract

In this paper, we present a molecular dynamics (MD) study on a low-energy ion implantation process for nanoscale CMOS (Complementary Metal Oxide Semiconductor) processes. To model the profiles of interstitials and vacancies, the recoil interaction approximation (RIA) was employed, while the kinetic Monte Carlo (KMC) approach was used for modeling the dynamic annealing effect between cascades. The simulation results performed for as-implanted boron profiles were compared with the results of the binary collision approximation (BCA) calculation by UT-MARLOWE as well as with the experimental SIMS data. The simulation revealed that the dynamic annealing effect between cascades is essential for the accurate estimation of defect distribution as well as as-implanted ion distribution. The dynamic annealing effect was carefully investigated for a case of boron implantation with an ion implantation energy of 2 keV, doses of 1×1014 ions/cm2 and 1×1015 ions/cm2, and a dose rate of 1×1012 ions/cm2·s.