Control of Large-size SiC Growth Process Using Modeling tool

Large size bulk silicon carbide (SiC) crystals are commonly grown by the physical vapor transport (PVT) method. The PVT growth of SiC crystals involves sublimation and condensation, chemical reactions, stoichiometry, mass transport, induced thermal stress, as well as defect and micropipes generation and propagation. The quality and polytype of as-grown SiC crystals are related to the temperature distribution inside the growth chamber during the growth process, it is critical to predict the temperature distribution from the measured temperatures outside the crucible by pyrometers. A radio-frequency induction-heating furnace was used for the growth of large-size SiC crystals by the PVT method in the present study. Modeling and simulation have been used to develop the SiC growth process and to improve the SiC crystal quality. Parameters such as the temperature measured at the top of crucible, temperature measured at the bottom of the crucible, and inert gas pressure are used to control the SiC growth process. By measuring the temperatures at the top and bottom of the crucible, the temperatures inside the crucible were predicted with the help of modeling tool. SiC crystals of 6H polytype were obtained and characterized by the Raman scattering spectroscopy and SEM, and crystals of few millimeter size grown inside the crucible were found without micropipes. Expansion of the crystals were also performed with the help of modeling and simulation.

3D Finite Volume Scheme for Czochralski Crystal Growth

A general three-dimensional model is developed for simulation of the growth process of silicon single crystals by Czochralski technique. The numerical scheme is based on the curvilinear non-orthogonal finite volume discretization. Numerical solutions show that the flow and temperature fields in the melt are asymmetric and unsteady for 8’’ silicon growth. The effects of rotation of crystal on the flow structure are studied. The rotation of crystal forms the Ekman layer in which the temperature gradient along solid/melt surface is small.

3D Finite Volume Scheme for Czochralski Single Crystal Growth

Czochralski (Cz) technique, which is used for growing single crystals, has dominated the production of single crystals for electronic applications. The Cz growth process involves multiple phases, moving interface and three-dimensional behavior. Much has been done to study these phenomena by means of numerical methods as well as experimental observations. A three-dimensional curvilinear finite volume based algorithm has been developed to model the Cz process. A body-fitted transformation based approach is adopted in conjunction with a multizone adaptive grid generation (MAGG) technique to accurately handle the three-dimensional problems of phase-change in irregular geometries with free and moving surfaces. The multizone adaptive model is used to perform a three-dimensional simulation of the Cz growth of silicon single crystals.Since the phase change interface are irregular in shape and they move in response to the solution, accurate treatment of these interfaces is important from numerical accuracy point of view. The multizone adaptive grid generation (MAGG) is the appropriate scheme for this purpose. Another challenge encountered is the moving and periodic boundary conditions, which is essential to the numerical solution of the governing equations. Special treatments are implemented to impose the periodic boundary condition in a particular direction and to determine the internal boundary position and shape varying with the combination of ambient physicochemical transport process and interfacial dynamics. As indicated above that the applications and processes characterized by multi-phase, moving interfaces and irregular shape render the associated physical phenomena three-dimensional and unsteady. Therefore a generalized 3D model rather than a 2D simulation, in which the governing equations are solved in a general non-orthogonal coordinate system, is constructed to describe and capture the features of the growth process. All this has been implemented and validated by using it to model the low pressure Cz growth of silicon. Accuracy of this scheme is demonstrated by agreement of simulation data with available experimental data. Using the quasi-steady state approximation, it is shown that the flow and temperature fields in the melt under certain operating conditions become asymmetric and unsteady even in the absence of extrinsic sources of asymmetry. Asymmetry in the flow and temperature fields, caused by high shear initiated phenomena, affects the interface shape in the azimuthal direction thus results in the thermal stress distribution in the vicinity, which has serious implications from crystal quality point of view.