Modelling of the isothermal melt flow due to rotating magnetic fields in crystal growth

The use of a rotating magnetic field promises the feature of a contactless flow control in crystal growth especially for configurations where an increase of the material transport in a definite way is desired. This paper gives the comparison of numerically calculated and experimentally obtained results on the flow due to a rotating magnetic field as well as numerical results on the influence of the field parameters (frequency, amplitude) on the fluid flow in the melt.     

Applicability of LES turbulence modeling for CZ silicon crystal growth systems with traveling magnetic field

To examine the applicability of LES turbulence modeling for CZ silicon crystal growth systems with traveling magnetic fields, LES calculations with Smagorinsky–Lilly turbulence model and van Driest damping at the solid walls are carried out. The program package for the calculations was developed on the basis of the open-source code library OpenFOAM^®${\mathrm{OpenFOAM}}^{\circledR }$. A previously published laboratory model with low temperature melt InGaSn, a 20” crucible, and process parameters corresponding to industrial Czochralski silicon systems is considered. Flow regimes with two crystal and crucible rotation rates and with different strengths of the traveling magnetic field “down” are analyzed. The calculated distributions of averaged temperature and standard temperature deviations are compared with measured ones in the laboratory system, and a relatively good agreement between them is shown. The influence of chosen time steps and grid sizes is analyzed by comparing Fourier spectra of temperature time-autocorrelation functions and temperature spatial distributions, and it is shown that the used moderate meshes of few hundred thousand cells can be applied for practical calculations.     

Effect of magnetic field on melt flow and crystal growth of oxide crystals

The flow in an oxide melt such as LiNbO, and TiO₂ in a high magnetic field was observed by using magnetic-field-applied Czochralski equipment for oxide crystals. It was found that the flows in oxides melts were very much different from these in a semiconductor melt. The single crystals of TiO₂ were grown in a magnetic field by using this equipment.     

Analysis of SiC crystal sublimation growth by fully coupled compressible multi-phase flow simulation

A fully coupled compressible multi-phase flow solver was developed to effectively design a large furnace for producing large-size SiC crystals. Compressible effect, convection and buoyancy effects, flow coupling between argon gas and species, and the Stefan effect are included. A small and experimental furnace is used to validate the solver. First, the essentiality of 2D flow calculation and the significance of incorporating buoyancy effect and gas convection, the Stefan effect, and flow interaction between argon gas and species were investigated by numerical results. Then the effects of argon gas on deposition rate, growth rate, graphitization on the powder source, and supersaturation and stoichiometry on the seed were analyzed. Finally, the advantages of an extra chamber design were explained, and improvement of growth rate was validated by the present solver.     

Silicon transport under rotating and combined magnetic fields in liquid phase diffusion growth of SiGe

The effect of applied rotating and combined (rotating and static) magnetic fields on silicon transport during the liquid phase diffusion growth of SiGe was experimentally studied. 72‐hour growth periods produced some single crystal sections. Single and polycrystalline sections of the processed samples were examined for silicon composition. Results show that the application of a rotating magnetic field enhances silicon transport in the melt. It also has a slight positive effect on flattening the initial growth interface. For comparison, growth experiments were also conducted under combined (rotating and static) magnetic fields. The processed samples revealed that the addition of static field altered the thermal characteristics of the system significantly and led to a complete melt back of the germanium seed. Silicon transport in the melt was also enhanced under combined fields compared with experiments with no magnetic field. (© 2010 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Numerical investigation of carbon and silicon carbide contamination during the melting process of the Czochralski silicon crystal growth

Carbon contamination in single crystalline silicon is detrimental to the minority carrier lifetime, one of the critical parameters for electronic wafers. In order to study the generation and accumulation of carbon contamination, transient global modeling of heat and mass transport was performed for the melting process of the Czochralski silicon crystal growth. Carbon contamination, caused by the presence of carbon monoxide in argon gas and silicon carbide in the silicon feedstock, was predicted by the fully coupled chemical model; the model included six reactions taking place in the chamber. A simplified model for silicon carbide generation by the reaction between carbon monoxide and solid silicon was proposed using the closest packing assumption for the blocky silicon feedstock. The accumulation of carbon in the melted silicon feedstock during the melting and stabilization stages was predicted. Owing to this initial carbon content in the melt, controlling carbon contamination before the growth stage becomes crucial for reducing the carbon incorporation in a growing crystal.     

Numerical simulation of thermal and mass transport during Czochralski crystal growth of sapphire

The thermal and flow transport in an inductively heated Czochralski crystal growth furnace during a crystal growth process is investigated numerically. The temperature and flow fields inside the furnace, coupled with the heat generation in the iridium crucible induced by the electromagnetic field generated by the RF coil, are computed. The results indicate that for an RF coil fixed in position during the growth process, although the maximum value of the magnetic, temperature and velocity fields decrease, the convexity of the crystal‐melt interface increases for longer crystal growth lengths. The convexity of the crystal‐melt interface and the power consumption can be reduced by adjusting the relative position between the crucible and the induction coil during growth. (© 2010 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Global simulation of coupled carbon and oxygen transport in a Czochralski furnace for silicon crystal growth

For accurate prediction of carbon and oxygen impurities in a single crystal produced by the Czochralski method, global simulation of coupled oxygen and carbon transport in the whole furnace was implemented. Both gas-phase transportation and liquid-phase transportation of oxygen and carbon were considered. With five chemical reactions considered, SiO and CO concentrations in gas and C and O atom concentrations in silicon melt were solved simultaneously. The simulation results show good agreement with experimental data.     

Time-dependent simulation of the growth of large silicon crystals by the Czochralski technique using a turbulent model for melt convection

A numerical model is developed to perform the dynamic and global simulation of Czochralski growth. The effect of melt convection is taken into account by means of an eddy viscosity flow model, which can represent the mixing effect of flow oscillations on the heat transfer. Our method is used to investigate the dynamics of the growth of a 40 cm silicon crystal.     

Influence of Coriolis force on thermal convection and impurity segregation during crystal growth under microgravity

In contrast with the generally accepted viewpoint, it is shown that the Coriolis force caused by rotation of an orbital station can appreciably affect natural convection and impurity distribution during the growth of crystals from a melt in orbital flight conditions. 2D and 3D steady and oscillatory convection in a rectangular enclosure is considered. The resonance phenomenon arising due to the interaction of the Coriolis force and harmonic oscillations of the gravity force is demonstrated. It is shown that for moderate values of the Ekman number the Coriolis force suppresses convection in one direction and amplifies it in the other, which in turn results in deformation of the impurity distribution over the cross-section of the crystal.