The diffusion equation for non-equilibrium interstitial impurity atoms in the form convenient for numerical solution was obtained. The proposed equation takes into account all different charge states of interstitial atoms and drift of all mobile charged species in the built-in electric field, although only the concentration of the neutral impurity interstitials must be obtained to solve the equation. The absorption of impurity interstitials by immobile sinks, mobile vacancies, and due to the passivation of dopant atoms was also taken into account. The derived effective coefficients that describe diffusion and trapping of impurity interstitials have the following characteristic features: (i) their dependences on the concentration of substitutionally dissolved dopant atoms were smooth and monotone; (ii) these concentration dependences were obtained in the form traditionally used for the approximation of effective diffusivity of substitutionally dissolved dopant atoms in processing-simulation codes. To illustrate the usability and efficiency of the derived equation for simulation of diffusion of non-equilibrium interstitial impurity atoms, the migration of deuterium interstitials in the silicon substrate doped with boron was simulated. The calculated profile of deuterium concentration in the passivated region agrees well with the experimental data.

Interstitial Diffusion under Conditions of Trapping of Interstitial Impurity Atoms. O.I.Velichko, Y.P.Shaman: Materials Science in Semiconductor Processing, 2010, 13[1], 13-20