Combined effects of crucible geometry and Marangoni convection on silicon Czochralski crystal growth

In order to understand the influence of crucible geometry combined with natural convection and Marangoni convection on melt flow pattern, temperature and pressure fields in silicon Czochralski crystal growth process, a set of numerical simulations was conducted. We carry out calculation enable us to determine temperature, pressure and velocity fields in function of Grashof and Marangoni numbers. The essential results show that the hemispherical geometry of crucible seems to be adapted for the growth of a good quality crystal and the pressure field is strongly affected by natural and Marangoni convection and it is more sensitive than temperature. (© 2009 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Three dimensional numerical study of the effect of large Grashof number on HEM crystal growth

Heat transfer and fluid flow in HEM crystal growth of silicon in cylindrical cavity is studied numerically. The walls of the crucible are heated to a fixed temperature. The exchanger that causes and induces natural convection is seated at the middle‐bottom of the crucible. The finite‐volume method is employed to solve the governing equations with proper boundary conditions. The effects of transport mechanism on the temperature distribution, melt flow, pressure and stream function are presented. We focus our work on the pressure field which has not yet been studied in HEM crucible. Also, we extend our work on a wide range Grashof number and for large numbers until 10¹² not yet studied in HEM furnace. It is found that the onset of flow fluctuations appears at Gr = 10¹⁰. Uniform temperature is observed in the entire melt at high Grashof number with development of a thermal boundary layer close to the exchanger. The thermal boundary layer thickness is calculated for strong buoyancy regime. Besides, for very high Gr number, buoyancy has less effect on temperature and then on melt‐crystal interface shape. During enlarging Gr, pressure evolution is related to temperature variation more than flow pattern.     

勾形磁场下提拉法生产单晶硅的数值模拟

本文给出了提拉单晶硅时,勾形磁场强度的计算公式,并对单晶硅在有无勾形磁场情况下熔体内流场和氧的浓度分布进行了数值模拟,计算出磁场作用下磁场强度和洛伦兹力及有无磁场时流函数、垂直截面处的速度场和氧的浓度分布.通过分析表明,勾形磁场能使流动更为平稳,能有效地降低熔体内及生长界面氧的浓度,并对产生这一现象的机理作了理论分析.     

大尺寸电阻加热式碳化硅晶体生长热场设计与优化

大尺寸低缺陷碳化硅(SiC)单晶体是功率器件和射频(RF)器件的重要基础材料,物理气相传输(physical vapor transport, PVT)法是目前生长大尺寸SiC单晶体的主要方法。获得大尺寸高品质晶体的核心是通过调节组分、温度、压力实现气相组分在晶体生长界面均匀定向结晶,同时尽可能减小晶体的热应力。本文对电阻加热式8英寸(1英寸=2.54 cm)碳化硅大尺寸晶体生长系统展开热场设计研究。首先建立描述碳化硅原料受热分解热质输运及其多孔结构演变、系统热输运的物理和数学模型,进而使用数值模拟方法研究加热器位置、加热器功率和辐射孔径对温度分布的影响及其规律,并优化热场结构。数值模拟结果显示,通过优化散热孔形状、保温棉的结构等设计参数,电阻加热式大尺寸晶体生长系统在晶锭厚度变化、多孔介质原料消耗的情况下均能达到较低的晶体横向温度梯度和较高的纵向温度梯度。     

Melt stirring effect of a weak magnetic field on crystal growth

Melt stirring effect of a weak magnetic field for the natural convection of liquid metal in an electrically adiabatic cubic enclosure heated from one vertical wall and cooled from an opposing wall was studied by a fully transient three-dimensional numerical analyses and the reasoning for melt stirring effect was clarified from the numerical results. Similar techniques were applied for the melt convection in a cylindrical Czochralski crystal growing crucible with an application of a vertical magnetic field. In a static crucible, central fluid column rotated in a magnetic field and in a rotating crucible, central fluid column did not rotate in a magnetic field. These peculiar characteristics could have been explained due to the Lorentz force.     

Experiments on the oscillatory convection of low Prandtl number liquid in Czochralski configuration for crystal growth with cusp magnetic field

This paper presents results of experiments on the oscillatory convection of mercury in a Czochralski configuration with cusp magnetic field. Temperature fluctuation measurements are carried out to determine the critical Rayleigh number for the onset of time dependent natural convection. The effects of a cusp magnetic field on the supercritical natural convection coupled with rotation of crystal disk are investigated. In the presence of a rotating flow it is found that a cusp magnetic field can induce a new long wave instability and can amplify the temperature fluctuation depending on the magnitude of the relevant flow similarity parameters and the melt aspect ratios. A flow regime diagram for the amplification and damping of the temperature fluctuations is presented to provide an experimental data base for finding optimum growth conditions in the cusp magnetic field Czochralski process.     

Three dimensional simulation of melt flow in Czochralski crystal growth with steady magnetic fields

Three-dimensional transient numerical simulations were carried out to investigate the melt convection and temperature fluctuations within an industrial Czochralski crucible. To study the magnetic damping effects on the growth process, a vertical magnetic field and a cusp magnetic field were considered. Due to our special interest in the melt convection, only local simulation was conducted. The melt flow was calculated by large-eddy simulation (LES) and the magnetic forces were implemented in the CFD code by solving a set of user-defined scalar (UDS) functions. In the absence of magnetic fields, the numerical results show that the buoyant plumes rise from the crucible to the free surface and the crystal–melt interface, which indicates that the heat and mass transfer phenomena in Si melt can be characterized by the turbulent flow patterns. In the presence of a vertical magnetic field, the temperature fluctuations in the melt are significantly damped, with the buoyant plumes forming regular cylindrical geometries. The cusp magnetic field could also markedly reduce the temperature fluctuations, but the buoyant plumes would break into smaller vortical structures, which gather around the crystal as well as in the center of the crucible bottom. With the present crucible configurations, it is found that the vertical magnetic field with an intensity of 128 mT can damp the temperature fluctuations more effectively than the 40 mT cusp magnetic field, especially in the region near the growing crystal.     

Direct observation and numerical simulation of molten silicon flow during crystal growth under magnetic fields by x-ray radiography and large-scale computation

Control of melt flow during Czochralski (CZ) crystal growth by application of magnetic fields is an important technique for large-diameter (>300 mm) silicon single crystals. Melt convection under magnetic fields is an interesting problem for electromagnetic-hydrodynamics. This paper reviews the effects of a vertical magnetic field and a cusp-shaped magnetic field on melt flow during CZ crystal growth. Melt flow in vertical magnetic fields or cusp-shaped magnetic fields was investigated by the direct observation method based on X-ray radiography and by numerical simulation. The first part of this review shows the result of direct observation of molten silicon flow under magnetic fields. It also compares the results of experimental and numerical simulation. The second part shows the details of the numerical simulation of the behavior of molten silicon in magnetic fields.     

Static magnetic fields in semiconductor floating-zone growth

Heat and mass transfer in semiconductor float-zone processing are strongly influenced by convective flows in the zone, originating from sources such as buoyancy convection, thermocapillary (Marangoni) convection, differential rotation, or radio frequency heating. Because semiconductor melts are conducting, flows can be damped by the use of static magnetic fields to influence the interface shape and the segregation of dopants and impurities. An important objective is often the suppression of time-dependent flows and the ensuing dopant striations. In RF-heated Si-FZ-crystals, fields up to 0.5Tesla show some flattening of the interface curvature and a reduction of striation amplitudes. In radiation-heated (small-scale) Si-FZ crystals, fields of 0.2–0.5Tesla already suppress the majority of the dopant striations. The uniformity of the radial segregation is often compromised by using a magnetic field, due to the directional nature of the damping. Transverse fields lead to an asymmetric interface shape and thus require crystal rotation (resulting in rotational dopant striations) to achieve a radially symmetric interface, whereas axial fields introduce a coring effect. A complete suppression of dopant striations and a reduction of the coring to insignificant values, combined with a shift of the axial segregation profile towards a more diffusionlimited case, are possible with axial static fields in excess of 1Tesla. Strong static magnetic fields, however, can also lead to the appearance of thermoelectromagnetic convection, caused by the interaction of thermoelectric currents with the magnetic field.     

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.     

旋转磁场对空间浮区内流场及浓度场的影响

采用全浮区模型数值研究了旋转磁场作用下熔区内热毛细对流流动特性,分析了磁场强度对流场及浓度场的影响.研究发现,无磁场时,熔体内杂质浓度场和流场呈现三涡胞对称振荡特征;温度场主要由扩散作用决定,呈对称分布.旋转磁场作用下,Ma数基本保持不变.当磁场强度B0≤1 mT时,熔体内杂质浓度场和流场与无磁场时结构类似,但旋转磁场的搅拌作用使得熔体内周期性振荡提前出现,且当旋转磁场产生的洛伦兹力相对较大时,表面张力产生的三维振荡对流得到很好地抑制.B0=5 mT时,周向波动被完全抑制,熔区内流场和浓度场呈二维轴对称分布.旋转磁场对熔体流动产生的轴向抑制作用和周向搅拌作用,都有助于熔体流动的稳定性、浓度分布以及温度分布的均匀性,从而有利于高质量晶体的生长.     

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)     

紊流模型模拟分析旋转对提拉大直径单晶硅的影响

本文采用紊流模型对提拉大直径单晶硅时,对晶体旋转、坩埚旋转及二者共同作用三种情况下,熔体内的流线、等温线、氧的浓度分布、紊流粘性系数、紊动能等作了数值模拟,发现晶体的旋转能提高氧的径向均匀性,紊流粘性系数和紊动能随着坩埚转速的提高先增加后下降.晶体坩埚同时旋转时并不能有效降低紊流粘性系数,但能使子午面上的流动受到抑制,等温线更为平坦,有利于晶体生长.     

Modulating dopant segregation in floating-zone silicon growth in magnetic fields using rotation

The feasibility of modulating dopant segregation using rotation for floating-zone silicon growth in axisymmetric magnetic fields is investigated through computer simulation. In the model, heat and mass transfer, fluid flow, magnetic fields, melt/solid interfaces, and the free surface are solved globally by a robust finte-volume/Newton's method. Different rotation modes, single- and counter-rotations, are applied to the growth under both axial and cusp magnetic fields. Under the magnetic fields, it is observed that dopant mixing is poor in the quiescent core region of the molten zone, and the weak convection there is responsible for the segregation. Under an axial magnetic field, moderate counter-rotation or crystal rotation improves dopant uniformity. However, excess counter-rotation or feed rotation alone results in more complicated flow structures, and thus induces larger radial segregation. For the cusp fields, rotation can enhance more easily the dopant mixing in the core melt and thus improve dopant uniformity.     

Influence of a high vertical magnetic field on Te dopant segregation in InSb grown by the vertical gradient freeze method

A vertical gradient freeze apparatus was set up to investigate the influence of a vertical magnetc field on Te dopant segregation in InSb. Te-doped InSb crystals were grown in the presence and absence of an 80.0 kG magnetic field. The axial profile of the Te concentration in the crystal grown in the magnetic field was observed to be more uniform than that grown without magnetic field, which was attributed to the influence of the high magnetic field on Te dopant segregation by reducing convection in the melt.     

采用延时修正法实施QUICK格式模拟提拉单晶体的生长

为了数值模拟提拉(又名Czochralski)法获得单晶体的生长过程,本文采用有限容积法离散控制方程,采用非均匀的交错网格避免不合理的振荡压力场,采用三阶精度QUICK (Quadratic Upwind Interpolation of Convective Kinematics)格式离散对流项,采用延时修正来实施QUICK格式获得满足主对角占优的代数方程组,采用SIMPLE(Semi-implicit Method for Pressure Linked Equations)算法耦合压力和速度场,给出了基于上述方法的方程、算法,并发展了程序,计算了Wheeler标准问题,计算结果与文献相当一致,同时本算法能模拟计算高葛拉晓夫数时的流动,显示出非均匀网格QUICK格式模拟晶体生长的优越性;另外本文将这一算法运用到单晶硅的数值模拟中,计算结果令人满意.     

Numerical modeling of frequency influence on the electromagnetic stirring of semiconductor melts

Alternating magnetic fields can be used in order to increase the level of convection and to mix the doped semiconductor alloys. A numerical analysis of the electromagnetic induced convection in GaInSb semiconductor melts is performed by using the software package CrysVUn. The magnetic field parameters are varied in order to obtain a maximum efficiency of the induced convection with a minimum quantity of the heat released in the melt. The influence of the electrical current frequency on the convection intensity is analyzed for samples with various radii (R = 0.5 – 3cm). Numerical procedure is validated by comparing the numerical results obtained in mercury samples with the experimental data given from the literature, which show a maximum stirring for a magnetic skin depth δ = 0.2R , in the case of a mercury sample with the radius R = 10 cm. This maximum corresponds to a shielding parameter R _ω = 40. Our numerical results show that the value of the shielding parameter for which the convection intensity reaches the maximum depends on the sample radius and increases when the sample radius increases. The results of this analysis are important in the case of samples with small radius, when a good mixing of the melt can be obtained for frequencies much lower than those corresponding to a shielding parameter R_ω = 40. (© 2006 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)     

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)     

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)