A computer simulation model for oxygen precipitation during crystal growth and wafer annealing of Czochralski silicon is constructed. In this model, the precipitate morphology is determined by minimizing the excess free energy. This model is connected with the formation models of grown-in defects, such as crystal-originated particles ͑COPs͒ in the vacancy-rich region or extrinsic stacking faults in the self-interstitials rich region. That is, the simulation of oxygen precipitation starts at 1000 or 900°C with the input data of residual point defect concentrations after COP or stacking-fault formation. The model well simulates oxygen precipitation behavior. For example, the simulated results of isothermal annealing at 700 and 800°C show that (i) the precipitate morphology is initially spherical, then changes to oblate spheroidal with an increase in annealing time, and (ii) the aspect ratio  is almost constant at (2.0-6.0) ϫ 10 Ϫ2 for longer than 200 h. The simulated  agrees approximately with experimentally determined  ϭ (0.8-2.5) ϫ 10 Ϫ2 of platelet precipitates reported in the literature. The values of both the density and the growth rate of precipitates agree for the simulations and experiments. The model also shows that oxygen precipitation during crystal growth strongly depends on the residual point defect concentrations after grown-in defect formation.Further advances in device processing require double side polished silicon wafers to achieve higher flatness and low particle generation. For such wafers, implementation of external gettering at the back surface ͑e.g., poly or sandblasting͒ is difficult to realize. Therefore, all candidates for 300 mm diam wafers, such as crystaloriginated particles ͑COPs͒ less polished wafers, high temperature annealed wafers, and epitaxial ͑epi͒ wafers require internal gettering ͑IG͒ by oxide precipitates. Because oxygen precipitation behavior depends not only on oxygen concentration but on the thermal history of the crystal ͑i.e., the pulling rate of the crystal, the crystal portion, etc.͒, the control of oxygen precipitation becomes far more difficult in comparison with that of 200 mm diam wafers.A computer simulation technique of oxygen precipitation behavior is a powerful tool applicable to various Czochralski ͑CZ͒ silicon wafers which undergo various device fabrication processes. 1-7 Hartzell et al. first reported the model, which describes the kinetics of oxygen precipitation by rate equations associated by a FokkerPlanck equation. 2 The continuity equations for point defect distribution and the strain relief by the emission of self-interstitials were incorporated in the model by Schrems 3 and Esfandyari et al. 4 Senkader et al. expanded the model to deal with oxygen precipitation associated with the formation of stacking faults. 5 However, the first model that included the nucleation process and took into account the effects of point defects was reported by Kobayashi. 6 The model showed that the precipitate nucleation occurs in a growing crystal in a temperature ran...