Samples of epitaxially grown $n$-type silicon have been implanted at room temperature with low doses $({10}^{6}--{10}^{9}\text{ }{\text{cm}}^{\ensuremath{-}2})$ of He, C, Si, and I ions using energies from 2.75 to 46 MeV. Deep level transient spectroscopy studies reveal that the generation of divacancy $({V}_{2})$ and vacancy-oxygen (VO) pairs has a distinct ion mass dependence. Especially, the doubly negative charge state of the divacancy, ${V}_{2}(=/\ensuremath{-})$, decreases in intensity with increasing ion mass compared to that of the singly negative charge state of the divacancy, ${V}_{2}(\ensuremath{-}/0)$. In addition, the measurements show also a decrease in the intensity of the level assigned to VO compared to that of ${V}_{2}(\ensuremath{-}/0)$ with increasing ion mass. Carrier capture cross-section measurements demonstrate a reduction in the electron capture rate with increasing ion mass for all the three levels ${V}_{2}(\ensuremath{-}/0)$, ${V}_{2}(=/\ensuremath{-})$, and VO; but a gradual recovery occurs with annealing. Concurrently, the strength of the ${V}_{2}(\ensuremath{-}/0)$ level decreases in a wide temperature range starting from below $200\text{ }\ifmmode^\circ\else\textdegree\fi{}\text{C}$, accompanied by an increase in the amplitudes of both the VO and ${V}_{2}(=/\ensuremath{-})$ peaks. In order to account for these results a model is introduced where local carrier compensation is a key feature and where two modes of ${V}_{2}$ are considered: (1) ${V}_{2}$ centers located in regions with a high defect density around the ion track $({V}_{2}^{\text{dense}})$ and (2) ${V}_{2}$ centers located in regions with a low defect density $({V}_{2}^{\text{dilute}})$. The ${V}_{2}^{\text{dense}}$ fraction does not give any contribution to the ${V}_{2}(=/\ensuremath{-})$ signal due to local carrier compensation, and the amplitude of the ${V}_{2}(=/\ensuremath{-})$ level is determined by the ${V}_{2}^{\text{dilute}}$ fraction only. The spatial distributions of defects generated by single-ion impacts were simulated by Monte Carlo calculations in the binary collision approximation, and to distinguish between the regions with ${V}_{2}^{\text{dense}}$ and ${V}_{2}^{\text{dilute}}$ a threshold for the defect generation rate was introduced. The model is shown to give good quantitative agreement with the experimentally observed ion mass dependence for the ratio between the amplitudes of the ${V}_{2}(=/\ensuremath{-})$ and ${V}_{2}(\ensuremath{-}/0)$ peaks. In particular, the threshold value for the defect generation rate remains constant $(\ensuremath{\sim}1.2\text{ }\text{vacancies}/\text{ion}/\text{\AA{}})$ irrespective of the type of ion used, which provides strong evidence for the validity of the model. Annealing at temperatures above $\ensuremath{\sim}300\text{ }\ifmmode^\circ\else\textdegree\fi{}\text{C}$ is found to reduce the spatial localization of the defects and migration of ${V}_{2}$ occurs with subsequent trapping by interstitial oxygen atoms and formation of divacancy-oxygen pairs.
Read full abstract