The Role of Marangoni Convection for the FZ-Growth of Silicon

Fluctuations of the electrical resistivity due to inhomogeneous dopant distribution are still a serious problem for the industrial processing yield of doped silicon crystals. In the case of silicon floating-zone growth, the main sources of these inhomogeneities are time- dependent flows in the liquid phase during the growth process. Excluding radio frequency (RF) induced convection, buoyancy and thermocapillary (Marangoni) convection are the two natural reasons for fluid flow. Both originate from temperature/concentration gradients in the melt, buoyancy convection through thermal/concentrational volume expansion, and thermocapillary convection through the temperature/concentration dependence of the surface tension. To improve the properties of grown crystals, knowledge of the strength, the characteristic, and the relation of these two flow mechanisms is essential. By the use of microgravity, the effect and the strength of buoyancy (gravity dependent) and thermocapillary (gravity independent) convection can be separated and clarified. Applying magnetic fields, both convective modes can be influenced: fluid flow can either be damped (static magnetic fields) or overlaid by a regular flow regime (rotating magnetic fields). Two complementary approaches have been pursued: Silicon full zones (experiments on the German sounding rockets TEXUS 7, 12, 22, 29, and 36) with the maximum temperature at half of the zone height and silicon half zones (experiments on the Japanese sounding rockets TR-IA4 and 6) with the maximum temperature at the top of the melt. With the full zone arrangement, the intensity and the frequency of the dopant striations could be determined and the critical Marangoni number could be identified. The half zone configuration is suited to classify the flow pattern and to measure the amplitude and the frequency of temperature fluctuations in the melt by inserting thermocouples or temperature sensors into the melt. All experiments have been carried out in monoellipsoid mirror furnaces. Typical zone geometries are approx. 8 to 14 mm in diameter and height. The crystals grown under microgravity are compared to crystals grown in static axial magnetic fields (B<5 tesla) and in transversal rotating magnetic fields (B<7.5 mT / f=50 Hz). Experimental results are completed by 3D numerical simulations: the obtained temperature and concentration distribution in the melt confirm the damping effect of rotating magnetic fields and explain the change in the radial segregation under static magnetic fields.