Characterization of a Czochralski grown silicon crystal doped with 1020 cm‐3 germanium

A dislocation‐free silicon single crystal doped with 10²⁰ cm^‐3 germanium (Ge) has been grown using the Czochralski (CZ) growth technique. The Ge concentration in the seed‐end and tang‐end of the crystal was 8×10¹⁹cm^‐3and 1.6×10²⁰ cm^‐3, respectively. The effective segregation coefficient of Ge, the distribution of flow pattern defects (FPDs) and the wafer warpage have been characterized. Both the effective segregation coefficient and the equilibrium segregation coefficient of Ge in silicon were evaluated. Then, the density of FPDs was traced from seed‐end to tang‐end of the ingot, a suppression of FPDs by Ge doping was shown. That is probably because the Ge atoms consume free vacancies and thus a higher density of smaller voids is formed. Furthermore, the mechanical strength of wafers has also been characterized by batch warpage analysis. The warpage in the seed‐end was larger than that in the tang‐end of the ingot, showing that the mechanical strength of wafers is enhanced by Ge doping. Such improvement is interpreted by an enhanced dislocation pinning effect associated with the enhanced nucleation of grown‐in oxygen precipitates in the Ge‐doped silicon wafers. (© 2011 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Time-dependent simulation of Czochralski silicon crystal growth

We have developed a detailed mathematical model and numerical simulation tools based on the streamline upwind/Petrov-Galerkin (SUPG) finite element formulation for the Czochralski silicon crystal growth. In this paper we consider the mathematical modeling and numerical simulation of the time-dependent melt flow and temperature field in a rotationally symmetric crystal growth environment. Heat inside the Czochralski furnace is transferred by conduction, convection and radiation, Radiating surfaces are assumed to be opaque, diffuse and gray. Hence the radiative heat exchange can be modeled with a non-local boundary condition on the radiating part of the surface. The position of the crystal-melt interface is solved by the enthalpy method. The melt flow is assumed to be laminar and governed by the cylindrically symmetric and incompressible Navier-Stokes equations coupled with the calculation of temperature.     

Global simulation of a Czochralski furnace for silicon crystal growth against the assumed thermophysical properties

In order to understand the effects of the thermophysical properties of the melt on the transport phenomena in the Czochralski (Cz) furnace for the single crystal growth of silicon, a set of global analyses of momentum, heat and mass transfer in small Cz furnace (crucible diameter: 7.2 cm, crystal diameter: 3.5 cm, operated in a 10 Torr argon flow environment) was carried out using the finite‐element method. The global analysis assumed a pseudosteady axisymmetric state with laminar flow. The results show that different thermophysical properties will bring different variations of the heater power, the deflection of the melt/crystal interface, the axial temperature gradient in the crystal on the center of the melt/crystal interface and the average oxygen concentration along the melt/crystal interface. The application of the axial magnetic field is insensitive to this effect. This analysis reveals the importance of the determination of the thermophysical property in numerical simulation. (© 2006 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Numerical investigation of heat transport and fluid flow during the seeding process of oxide Czochralski crystal growth Part 1: non‐rotating seed

For the seeding process of oxide Czochralski crystal growth, the flow and temperature field of the system as well as the seed‐melt interface shape have been studied numerically using the finite element method. The configuration usually used initially in a real Czochralski crystal growth process consists of a crucible, active afterheater, induction coil with two parts, insulation, melt, gas and non‐rotating seed crystal. At first the volumetric distribution of heat inside the metal crucible and afterheater inducted by the RF coil was calculated. Using this heat source the fluid flow and temperature field were determined in the whole system. We have considered two cases with respect to the seed position: (1) before and (2) after seed touch with the melt. It was observed that in the case of no seed rotation (ω_seed = 0), the flow pattern in the bulk melt consists of a single circulation of a slow moving fluid. In the gas domain, there are different types of flow motion related to different positions of the seed crystal. In the case of touched seed, the seed‐melt interface has a deep conic shape towards the melt. It was shown that an active afterheater and its location with respect to the crucible, influences markedly the temperature and flow field of the gas phase in the system and partly in the melt. (© 2007 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Measurement, modelling and simulation of defects in as-grown Czochralski silicon

Defects in as-grown and heat-treated 150 and 200 mm Czochralski silicon crystals are investigated for different crystal pulling conditions and thermal treatments. First results are presented using noncontact carrier recombination imaging for detection of electrically active defects. The defect densities and sizes are compared with the results of computer modelling, suggesting that the observed defects in as-grown material are most likely large voids, possibly partially filled with silicon oxide. In contrast, the defects observed after prolonged heat treatments are silicon oxide precipitates with a density which is several orders of magnitude larger. The voids nucleate at temperatures above 1100°C by a homogeneous nucleation process and grow further to the observed size during further cooling of the crystal. They are responsible for the midfield breakdown of 10–50 nm gate oxides.     

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.     

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)     

Three‐dimensional study of the pressure field and advantages of hemispherical crucible in silicon Czochralski crystal growth

The effects of several growth parameters in cylindrical and spherical Czochralski crystal process are studied numerically and particularly, we focus on the influence of the pressure field. We present a set of three‐dimensional computational simulations using the finite volume package Fluent in two different geometries, a new geometry as cylindro‐spherical and the traditional configuration as cylindro‐cylindrical. We found that the evolution of pressure which is has not been studied before; this important function is strongly related to the vorticity in the bulk flow, the free surface and the growth interface. It seems that the pressure is more sensitive to the breaking of symmetry than the other properties that characterize the crystal growth as temperature or velocity fields. (© 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.     

Three‐dimensional numerical investigation of the effects of an open window on the convective heat transfer of an oxide Czochralski system

A three‐dimensional numerical analysis was carried out for a real Czochralski crystal growth furnace containing only gas and without any melt and crystal in order to investigate the effects of a small observation window on the temperature and flow field of the system. For this approach, the induction heating equations, the Navier‐Stokes equation with Boussinesq approximation, the continuity and energy equations have been solved in cylindrical coordinates using the finite element method. It has been found that the flow and thermal fields in the system are obviously three‐dimensional and non‐axisymmetric. The gas enters the system through the window is directed towards the opposite side wall where it is divided into two parts of vertical direction as well as expands in horizontal direction. Consequently, there is a spiral gas flow in the crucible and afterheater which rotates upwards in azimuthal direction along the walls. (© 2007 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Numerical investigation of heat transport and fluid flow during the seeding process of oxide Czochralski crystal growth Part 2: rotating seed

In this paper, the role of seed rotation on the characteristics of the two‐dimensional temperature and flow field in the oxide Czochralski crystal growth system has been studied numerically for the seeding process. Based on the finite element method, a set of two‐dimensional quasi‐steady state numerical simulations were carried out to analyze the seed‐melt interface shape and heat transfer mechanism in a Czochralski furnace with different seed rotation rates: ω_seed = 5‐30 rpm. The results presented here demonstrate the important role played by the seed rotation for influencing the shape of the seed‐melt interface during the seeding process. The seed‐melt interface shape is quite sensitive to the convective heat transfer in the melt and gaseous domain. When the local flow close to the seed‐melt interface is formed mainly due to the natural convection and the Marangoni effect, the interface becomes convex towards the melt. When the local flow under the seed‐melt interface is of forced convection flow type (seed rotation), the interface becomes more concave towards the melt as the seed rotation rate (ω_seed) is increased. A linear variation of the interface deflection with respect to the seed rotation rate has been found, too. (© 2007 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Influence of the crucible bottom shape on the heat transport and fluid flow during the seeding process of oxide Czochralski crystal growth

For the seeding process of oxide Czochralski crystal growth, influence of the crucible bottom shape on the heat generation, temperature and flow field of the system and the seed‐melt interface shape have been studied numerically using the finite element method. The configuration usually used in a real Czochralski crystal growth process consists of a crucible, active afterheater, induction coil with two parts, insulation, melt, gas and seed crystal. At first, the volumetric distribution of heat inside the metal crucible and afterheater inducted by the RF‐coil was calculated. Using this heat generation in the crucible wall as a source the fluid flow and temperature field of the entire system as well as the seed‐melt interface shape were determined. We have considered two cases, flat and rounded crucible bottom shape. It was observed that using a crucible with a rounded bottom has several advantages such as: (i) The position of the heat generation maximum at the crucible side wall moves upwards, compared to the flat bottom shape. (ii) The location of the temperature maximum at the crucible side wall rises and as a result the temperature gradient along the melt surface increases. (iii) The streamlines of the melt flow are parallel to the crucible bottom and have a curved shape which is similar to the rounded bottom shape. These important features lead to increasing thermal convection in the system and influence the velocity field in the melt and gas domain which help preventing some serious growth problems such as spiral growth. (© 2007 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Thermocapillary‐buoyancy flow of silicon melt in a shallow annular pool

In order to understand the nature of surface spoke patterns on silicon melt in industrial Czochralski furnaces, a series of unsteady three‐dimensional numerical simulations were conducted for thermocapillary‐buoyancy flow of silicon melt in annular pool (inner radius r_i = 15 mm, outer radius r_o = 50 mm, depth d = 3 mm). The pool is heated from the outer cylindrical wall and cooled at the inner wall. Bottom and top surfaces either are adiabatic or allow heat transfer in the vertical direction. Results show that a small temperature difference in the radial direction generates steady roll‐cell thermocapillary‐buoyancy flow. With large temperature difference, the simulation can predict three‐dimensional oscillatory flow, which is characterized by spoke patterns traveling in the azimuthal direction. The small vertical heat flux (3 W/cm²) does not have significant effects on the characteristics of this oscillatory flow. Details of the flow and temperature disturbances are discussed and the critical conditions for the onset of the oscillatory flow are determined. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Global heat loss and thermal stress analysis in Czochralski crystal growth

Based on our invention of an energy‐efficient Czochalski crystal growth furnace, a 2D‐axisymmetric numerical simulation model of LiNbO₃ crystal growth is developed. The heat transfer, melt and gas flow, radiation and the interface deflection have been examined. Heat losses in the furnace and the insulator, as well as the heating power and thermal stress distribution at three stages of crystal growth are calculated in detail. It is found that a large proportion of heat dissipates through the water‐cooling system, and at the steel shell of the furnace, gas convection heat transfer is the major cooling mechanism. Less heat dissipation by radiation and more heat flux by gas convection to the crystal sidewall results in a larger concentrated thermal stress, which may induce large crystal cracks in the growth process. The simulation results of heating power are in coincidence with the actual power of our furnace, which verifies the feasibility of our model. The detailed information with respect to the device obtained from simulation can help to optimize the energy‐saving design and growth process.     

Numerical investigation of crystal growth process of bulk Si and nitrides – a review

Heat and mass transfer during crystal growth of bulk Si and nitrides by using numerical analysis was studied. A three‐dimensional analysis was carried out to investigate temperature distribution and solid‐liquid interface shape of silicon for large‐scale integrated circuits and photovoltaic silicon. The analysis enables prediction of the solid‐liquid interface shape of silicon crystals. The result shows that the interface shape became bevel like structure in the case without crystal rotation. We also carried out analysis of nitrogen transfer in gallium melt during crystal growth of gallium nitride using liquid‐phase epitaxy. The result shows that the growth rate at the center was smaller than that at the center. (© 2007 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Numerical investigation of induction heating and heat transfer in a SiC growth system

A global simulation model is applied for a silicon carbide growth system heated by induction coils. A finite‐volume method (FVM) and a global model are applied to solve the equations for electromagnetic field, conductive and radiative heat transfer. The growth rate is predicted by Hertz‐Knudsen equation and one‐dimensional mass transfer equation. Further, simulations for five different coil positions are carried out to investigate the effects of coil position on temperature distribution in the furnace. The numerical results reveal that the variation of temperature in the radial direction along the substrate surface and the temperature difference between the powder and substrate are greatly affected by the coil position. The predicted growth rate along the substrate surface for five coil positions is also studied. Finally, a reasonable range of coil positions maintaining a balance between large‐diameter crystal, high growth rate, temperature limitation of material and lower electrical power consumption is obtained. (© 2007 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Influence of internal radiation on the heat transfer during growth of YAG single crystals by the Czochralski method

Heat and mass transfer taking place during growth of Y₃Al₅O₁₂ (YAG) crystals by the Czochralski method, including inner radiation, is analyzed numerically using a Finite Element Method. For inner radiative heat transfer through the crystal the band approximation model and real transmission characteristics, measured from obtained crystals, are used. The results reveal significant differences in temperature and melt flow for YAG crystals doped with different dopands influencing the optical properties of the crystals. When radiative heat transport through the crystal is taken into account the melt‐crystal interface shape is different from that when the radiative transport is not included. Its deflection remains constant over a wide range of crystal rotation rates until it finally rapidly changes in a narrow range of rotation rates. (© 2003 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

ACRT forced convection and its effects on solute segregation and heat and mass transfer during single crystal growth

A numerical simulation study was carried out for CdZnTe vertical Bridgman method crystal growth with the accelerated crucible rotation technique (ACRT). The convection, heat and mass transfer in front of the solid‐liquid interface, and their effects on the solute segregation of the grown crystal can be characterized with the following. ACRT brings about a periodic forced convection in the melt, of which the intensity and the incidence are far above the ones of the natural convection without ACRT. This forced convection is of multiformity due to the changes of the ACRT parameters. It can result in the increases of both the solid‐liquid interface concavity and the temperature gradient of the melt in front of the solid‐liquid interface, of which magnitudes vary from a little to many times as the ACRT wave parameters change. It also enhances the mass transfer in the melt in a great deal, almost results in the complete uniformity of the solute distribution in the melt. With suitable wave parameters, ACRT forced convection decreases the radial solute segregation of the crystal in a great deal, even makes it disappear completely. However, it increases both the axial solute segregation and the radial one notably with bad wave parameters. An excellent single crystal could be gotten, of which the most part is with no segregation, by adjusting both the ACRT wave parameters and the crystal growth control parameters, e.g. the initial temperature of the melt, the temperature gradient, and the crucible withdrawal rate. (© 2008 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Numerical investigation of magnetic field effect on pressure in cylindrical and hemispherical silicon CZ crystal growth

The effect of axial magnetic field of different intensities on pressure in silicon Czochralski crystal growth is investigated in cylindrical and hemispherical geometries with rotating crystal and crucible and thermocapillary convection. As one important thermodynamic variable, the pressure is found to be more sensitive than temperature to magnetic field with strong dependence upon the vorticity field. The pressure at the triple point is proposed as a convenient parameter to control the homogeneity of the grown crystal. With a gradual increase of the magnetic field intensity the convection effect can be reduced without thermal fluctuations in the silicon melt. An evaluation of the magnetic interaction parameter critical value corresponding to flow, pressure and temperature homogenization leads to the important result that a relatively low axial magnetic field is required for the spherical system comparatively to the cylindrical one.     

Boron segregation control in silicon crystal ingots grown in Czochralski process

Boron-doped silicon single crystals of 207 mm diameter with various growing conditions are grown from a large amount of the melt in the cusp-magnetic Czochralski method, and the effects of growing parameters on dopant concentrations in the crystals are experimentally investigated. Equilibrium distribution coefficient of boron calculated by BPS model is 0.716. With the crystal rotation (ω) of 13 rpm and the crucible rotation of , the effective distribution coefficient (k_e) is 0.751 in zero magnetic strength and increases up to 0.78 in the magnetic strength of 640 G. For , there is no significant influence of ω on k_e. With , k_e is almost unity. The experimental results are compared with theory.