Advances in Crystallization Processes 612 microvoids as a result of aggregation of vacancies (Goethem et al., 2008;Kulkarni, 2008a). In this physical model, the interaction between the impurities and intrinsic point defects is not considered (Kulkarni et. al., 2004). R e c e n t v e r s i o n s o f t h i s m o d e l h a v e s u g g e s t e d t h a t p a r t o f t h e v a c a n c i e s ( v ) i n t h e temperature range 1683 ... 1373 K, due to the interaction with oxygen (O) and nitrogen (N) impurities, are bound into complexes of the vO, vO 2 , and vN types (Kulkarni 2007;2008b). After the formation of microvoids, the aforementioned complexes grow and take up vacancies. This model has ignored the growth of the complexes by means of the injection of intrinsic interstitial silicon atoms and the interaction of an impurity with intrinsic interstitial silicon atoms (Kulkarni 2007;2008b).In the general case recombination-diffusion model assumes that the process of defect formation in dislocation-free silicon single crystals occurs in four stages: (i) fast recombination of intrinsic point defects near the crystallization front; (ii) the formation in the narrow temperature range 1423...1223 K depending on the value of V g /G microvoids or interstitial dislocation loops; (iii) the formation of oxygen clusters in the temperature range 1223...1023 K; (iv) growth of precipitates as a result of subsequent heat treatments.Recombination-diffusion model is the physical basis for models of the dynamics of point defects. The mathematical model of point defect dynamics in silicon quantitatively explains the homogeneous mechanism of formation of microvoids and dislocation loops. It should be noted that, in the general case, the model of point defect dynamics includes three approximations: rigorous, simplified, and discrete-continuum approaches (Sinno, 1999;Dornberger et. al., 2001;Wang & Brown, 2001;Kulkarni et. al., 2004;Kulkarni, 2005;Prostomolotov & Verezub, 2009). The rigorous model requires the solution to integrodifferential equations for point defect concentration fields, and the distribution of grown-in microdefects in this model is a function of the coordinates, the time, and the time of evolution of the size distribution of microdefects. A high consumption of time and cost for the performance of calculations required the development of a simplified model in which the average defect radius is approximated by the square root of the average defect area. This approximation is taken into account in the additional variable, which is proportional to the total area of the defect surface. The simplified model is effective for calculating the twodimensional distribution of grown-in microdefects. Both models use the classical nucleation theory and suggest the calculation of the formation of stable nuclei and the kinetics of diffusion-limited growth of defects. The discrete-continuum approximation suggests a complex approach: the solution to discrete equations for the smallest defects and the solution to the Fokker-Planck equation for...