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.     

System of Mathematical Models for the Analysis of Industrial FZ-Si-Crystal Growth Processes

A system of coupled mathematical models and the corresponding program package is developed to study the interface shape, heat transfer, thermal stresses, fluid flow as well as the transient dopant segregation in the floating zone (FZ) growth of large silicon crystals (diameter more than 100mm) grown by the needle-eye technique. The floating zone method with needle-eye technique is used to produce high-purity silicon single crystals for semiconductor devices to overcome the problems resulting from the use of crucibles. The high frequency electric current induced by the pancake induction coil, the temperature gradients and the feed/crystal rotation determine the free surface shape of the molten zone and cause the fluid motion. The quality of the growing crystal depends on the shape of the growth interface, the temperature gradients and corresponding thermal stresses in the single crystal, the fluid flow, and especially on the dopant segregation near the growth interface. From the calculated transient dopant concentration fields in the molten zone the macroscopic and microscopic resistivity distribution in the single crystal is derived. The numerical results of the resistivity distributions are compared with the resistivity distributions measured in the grown crystal.     

Numerical simulation and validation of the Peltier pulse marking of solid/liquid interfaces

We describe a numerical approach of the solidification of binary alloys to study the motion of a crystal/melt interface submitted to current pulses involving a modification of the dopant concentration field. For the thermal aspect, the Thomson effect, the Peltier effect and Joule heating have been included in the heat flow. For the solutal segregation, our model is based on mass transports which occur in the liquid phase, namely diffusion and convection. Numerical computations are validated by comparison with experimental data and thus could find applications in the prediction of the effects of Peltier pulse marking in crystal growth.     

Optimization of control parameters of cadmium zinc telluride Bridgman single crystal growth

The temperature gradient within a furnace chamber and the crucible pull rate are the key control parameters for cadmium zinc telluride Bridgman single crystal growth. Their effects on the heat and mass transfer in front of the solid‐liquid interface and the solute segregation in the grown crystal were investigated with numerical modeling. With an increase of the temperature gradient, the convection intensity in the melt in front of the solid‐liquid interface increases almost proportionally to the temperature gradient. The interface concavity decreases rapidly at faster crucible pull rates, while it increases at slow pull rates. Moreover, the solute concentration gradient in the melt in front of the solid‐liquid interface decreases significantly, as does the radial solute segregation in the grown crystal. In general, a decrease of the pull rate leads to a strong decrease of the concavity of the solid‐liquid interface and of the radial solute segregation in the grown crystal, while the axial solute segregation in the grown crystal increases slightly. A combination of a low crucible pull rate with a medium temperature gradient within the furnace chamber will make the radial solute segregation of the grown crystal vanish. (© 2007 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.     

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)     

Distribution coefficient of boron in Si crystal ingots grown in cusp-magnetic Czochralski process

Silicon single crystals are grown by the Czochralski method with various growing conditions. Effective segregation coefficient of boron is found to depend on the magnetic field in cusp-magnetic Cz method. Effects of zero-Gauss plane (ZGP), ZGP shape and magnetic intensity (MI) on the dopant concentration and its distribution in the crystal are experimentally investigated. The shape of ZGP is not only flat but also parabolic due to the magnetic ratio (MR), which is the ratio of the lower to upper electric-current densities in the configurations of the cusp-magnetic field. Equilibrium distribution coefficient of boron calculated by BPS model is 0.698. With the crystal rotation (w) of 16 rpm and the crucible rotation of ?0.5 rpm, the effective distribution coefficient (k_e) is 0.728 in zero magnetic intensity and increases up to 0.8093 in the parabolic ZGP shape. Although the magnetic strength near the crystal–melt interface decreases with increasing MR, it increases in the bulk melt, and hence k_e increases. Flow stability in the bulk melt influences k_e. At the magnetic field and growing conditions, k_e increases with increasing initial charge size of the silicon melt. There is no significant influence of ZGP on the radial distribution of the boron concentration. Simulation results of melt flow in the presence of a parabolic ZGP are outlined, and the segregation results in the experiments are compared with published experimental data.     

The numerical simulation of heat and mass transfer during growth of a binary system by the travelling heater method

The influence of convection and heat and mass transfer on the shape and position of melt/solid interfaces and on radial composition segregation is analysed numerically for the travelling heater method growth of a binary alloy in a vertical transparent ampoule. Results are presented for crystal and melt with thermophysical properties similar to Cd_xHg_1−xTe with the assumption that the pseudobinary CdTe-HgTe phase diagram is true. The two-dimensional axisymmetric heat transfer equation, hydrodynamical equation and convective diffusion equation are included in the mathematical model. The rates of crystal growth and dissolution are supposed to be proportional to the compositional supercooling in the melt near the interfaces. It is shown for the conditions when convection is absent that the interfaces are asymmetrically positioned respectively to the heater centre line. Intensive convection makes their position more symmetrical but the length of the liquid zone greater. The flow pattern in the melt appears to be greatly influenced by solutal gravitational convection. The nonlinear dependence of the melt density on the temperature and composition are used in the model. The cases when speed of the heater is antiparallel (stable density stratification) or parallel (unstable stratification) to the vector of gravitational acceleration are considered.     

蓝宝石晶体热性能的各向异性对SAPMAC法晶体生长的影响

采用有限元法对冷心放肩微量提拉法蓝宝石晶体生长过程中晶体内的温度、应力分布进行了模拟计算,结合实验结果讨论了蓝宝石晶体热性能的各向异性对晶体生长的影响.研究结果表明,对于冷心放肩微量提拉蓝宝石晶体生长系统,较大的轴向热导率有利于提高晶体的生长速率和界面稳定性,而稍大的径向热导率则有利于保持微凸的生长界面.晶体内的热应力受径向热膨胀系数的影响显著,随着径向热膨胀系数的增大而增大,最大热应力总是出现在籽晶与新生晶体的界面区域.在实验中选α轴为结晶取向,成功生长出了直径达230mm、高质量蓝宝石晶体.     

半导体单晶生长过程中的位错研究

阐述了现有的半导体单晶位错模型,即临界切应力模型和粘塑性模型的基本理论及应用状况.分析了熔体法单晶生长过程中影响位错产生、增殖的各种因素,以及抑制位错增殖的措施.与熔体不润湿、与晶体热膨胀系数相近的坩埚材料,低位错密度的籽晶可有效地抑制生长晶体的位错密度;固液界面的形状及晶体内的温度梯度是降低位错密度的关键控制因素,而两因素又受到炉膛温度梯度、长晶速率、气体和熔体对流等晶体生长工艺参数的影响.最后,对熔体单晶生长过程的位错研究进行了展望.     

Simulation of temperature and flow fields in an inductively heated melt growth system

The goal of the research presented here is to apply a global analysis of an inductively heated Czochralski furnace for a real sapphire crystal growth system and predict the characteristics of the temperature and flow fields in the system. To do it, for the beginning stage of a sapphire growth process, influence of melt and gas convection combined with radiative heat transfer on the temperature field of the system and the crystal‐melt interface have been studied numerically using the steady state two‐dimensional finite element method. For radiative heat transfer, internal radiation through the grown crystal and surface to surface radiation for the exposed surfaces have been taken into account. The numerical results demonstrate that there are a powerful vortex which arises from the natural convection in the melt and a strong and large vortex that flows upwards along the afterheater side wall and downwards along the seed and crystal sides in the gas part. In addition, a wavy shape has been observed for the crystal‐melt interface with a deflection towards the melt. (© 2010 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Analysis of the dopant segregation effects at the floating zone growth of large silicon crystals

A computer simulation is carried out to study the dopant concentration fields in the molten zone and in the growing crystal for the floating zone (FZ) growth of large (> 100 mm) Si crystals with the needle-eye technique and with feed/crystal rotation. The mathematical model developed in the previous work is used to calculate the shape of the molten zone and the velocity field in the melt. The influence of melt convection on the dopant concentration field is considered. The significance of the rotation scheme of the feed rod and crystal on the dopant distribution is investigated. The calculated dopant concentration directly at the growth interface is used to determine the normalized lateral resistivity distribution in the single crystal. The calculated resistivity distributions are compared with lateral spreading resistivity measurements in the single crystal.     

热屏优化对大直径单晶硅生长影响的数值模拟

通过对28英寸热场生长300 mm硅单晶过程中结晶速率、固液界面形状、晶体中热应力及晶体中氧含量的数值计算提出了在该热场条件下热屏的优化方案。数值计算结果表明:对热屏底端与晶体表面和熔体自由液面的距离以及热屏材料(优化前热屏使用单一石墨材料,优化后采用辐射率较高的内壁材料结合反射率较高的外壁材料组成复合式热屏)的优化可以减少主加热器对晶体的热辐射使得固液界面更加平坦,藉此增加结晶速率,减小晶体内热应力和熔体中氧含量。     

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)     

The role of inner and internal radiation on the melt growth of sapphire crystal

In this paper, for an inductively heated Czochralski furnace used to grow sapphire single crystal, influence of the inner (wall‐to‐wall) and crystal internal (bulk) radiation on the characteristics of the growth process such as temperature and flow fields, structure of heat transfer and crystal‐melt interface has been studied numerically using the 2D quasi‐steady state finite element method. The obtained results of global analysis demonstrate a strong dependence of thermal field, heat transport structure and crystal‐melt interface on both types of radiative heat transfer within the growth furnace.     

Global analysis of heat transfer in CZ crystal growth of oxide taking into account three-dimensional unsteady melt convection: Effect of meniscus shape

A global model of heat transfer analysis for the Czochralski crystal growth of oxides, in which a three-dimensional unsteady melt flow was taken into account, was developed in our recent study. In the model, however, the focal point was the methodology of formulating the model by coupling a conventional global model of heat transfer, which is based on a pseudosteady axisymmetric assumption, with a model of a three-dimensional, unsteady melt flow via two interface models. Therefore, for simplicity, the shape of the melt free surface was assumed to be flat. In this study, the global model was further developed by considering the meniscus of the melt free surface. It was found that the meniscus of the melt free surface caused the melt flow to be more unstable and shifted the critical Reynolds number at which the melt/crystal interface inversion occurs toward a much lower value.     

Dopant transport during semiconductor crystal growth: Axial versus transverse magnetic fields

The present study investigates the effects of magnetic field orientation, magnetic field strength and growth rate on the dopant segregation in semiconductor crystals, and presents results of dopant composition in the crystal and in the melt at several different times during growth for several combinations of process parameters. The crystal's lateral segregation depends on the magnetic field's orientation and strength while the axial segregation depends on the magnetic field's strength and the growth rate. If either convective or diffusive transport truly dominates, then the crystal's dopant distribution is laterally uniform. The axial distribution in the crystal approaches the well-mixed limit if the melt motion is strong and the growth rate is slow, and the distribution approaches the diffusion-controlled limit if the melt motion is slow and the growth rate is fast. The deviations of the dopant distribution in the crystal from lateral uniformity and from the classical limits are quantified for several combinations of process parameters.     

Effect of the growth rate on thermal conditions during the growth of single crystals with melt feeding

An attempt is made to establish a correlation between the radial and axial growth rates and the change in the conditions of heat transfer from a growing crystal to the atmosphere of the water-cooled vacuum furnace for the growth of large alkali halide single crystals. It is found experimentally that an increase in the growth rate leads to an increase in the automatic compensation of the melt temperature by the main heater. In this case, the thickness of the layer of melt condensate on the end face and the lateral surface of the crystal decreases. It is revealed that the possibility of growing infinitely long ingots in the presence of intense melt evaporation is restricted by the possibilities of the heat transport through the boundary between the furnace atmosphere and the cooled furnace walls, onto which melt condensate deposits.     

Radial segregation due to weak convection in a floating zone

The influence of weak convection, caused by surface tension forces, on radial dopant segregation occurring in crystals grown under microgravity conditions is studied numerically. The geometry considered corresponds to a floating-zone configuration with partially coated melt surfaces consisting of small evenly distributed spots of free surfaces. In order to distinguish dopant distribution due to weak convection clearly from distribution due to diffusion the spots only cover one quarter of the periphery. Thus, surface tension-driven convection is allowed only over one quarter of the floating-zone configuration resulting in an asymmetric dopant distribution. The percentage of free surfaces present is varied in order to alter the Marangoni flow rates. The maximum dopant concentration due to radial segregation is plotted as a function of a certain convection level. The results of the present numerical study are supposed to be used to design corresponding space experiments launched at the end of the year 2000.     

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.