Numerical Study on Flow Field and Temperature Distribution in Growth Process of 200 mm Czochralski Silicon Crystals

The melt flow and temperature distribution in a 200 mm silicon Czochralski furnace with a cusp magnetic field was modeled and simulated by using a finite-volume based FLUTRAPP ( Fluid Flow and Transport Phenomena Program) code. The melt flow in the crucible was focused, which is a result of the competition of buoyancy, the centrifugal forces caused by the rotations of the crucible and crystal, the thermocapillary force on the free surfaces and the Lorentz force induced by the cusp magnetic field. The zonal method for radiative heat transfer was used in the growth chamber, which was confined by the crystal surface, melt surface, crucible, heat shield, and pull chamber. It was found that the cusp magnetic field could strength the dominant counter-rotating swirling flow cell in the crucible and reduce the flow oscillation and the pulling-rate fluctuation. The fluctuation of dopant and oxygen concentration in the growing crystal could thus be smoothed.

Effects of Thermocapillary Convection on the Growth Process in Industrial 8-inch Silicon Czochralski Growth

Czochralski (CZ) crystal growth process is a widely used technique in manufacturing of silicon crystals and other semiconductor materials. The ultimate goal of the IC industry is to have the highest quality substrates, which are free of point defect, impurities and micro defect clusters. The scale up of silicon wafer size from 200 mm to 300 mm requires large crucible size and more heat power. Transport phenomena in crystal growth processes are quite complex due to melt and gas flows that may be oscillatory and/or turbulent, coupled convection and radiation, impurities and dopant distributions, unsteady kinetics of the growth process, melt crystal interface dynamics, free surface and meniscus, stoichiometry in the case of compound materials. A global model has been developed to simulate the temperature distribution and melt flow in an 8-inch system. The present program features the fluid convection, magnetohydrodynamics, and radiation models. A multi-zone method is used to divide the Cz system into different zones, e.g., the melt, the crystal and the hot zone. For calculation of temperature distribution, the whole system inside the stainless chamber is considered. For the convective flow, only the melt is considered. The widely used zonal method divides the surface of the radiation enclosure into a number of zones, which has a uniform distribution of temperature, radiative properties and composition. The integro-differential equations for the radiative heat transfer are solved using the matrix inversion technique. The zonal method for radiative heat transfer is used in the growth chamber, which is confined by crystal surface, melt surface, heat shield, and pull chamber. Free surface and crystal/melt interface are tracked using adaptive grid generation. The competition between the thermocapillary convection induced by non-uniform temperature distributions on the free surface and the forced convection by the rotation of the crystal determines the interface shape, dopant distribution, and striation pattern. The temperature gradients on the free surface are influenced by the effects of the thermocapillary force on the free surface and the rotation of the crystal and the crucible.

Progress In Modeling Of Fluid Flows In Crystal Growth Processes

Modeling of fluid flows in crystal growth processes has become an important research area in theoretical and applied mechanics. Most crystal growth processes involve fluid flows, such as flows in the melt, solution or vapor. Theoretical modeling has played an important role in developing technologies used for growing semiconductor crystals for high performance electronic and optoelectronic devices. The application of devices requires large diameter crystals with a high degree of crystallographic perfection, low defect density and uniform dopant distribution. In this article, the flow models developed in modeling of the crystal growth processes such as Czochralski, ammonothermal and physical vapor transport methods are reviewed. In the Czochralski growth modeling, the flow models for thermocapillary flow, turbulent flow and MHD flow have been developed. In the ammonothermal growth modeling, the buoyancy and porous media flow models have been developed based on a single-domain and continuum approach for the composite fluid-porous layer systems. In the physical vapor transport growth modeling, the Stefan flow model has been proposed based on the flow-kinetics theory for the vapor growth. In addition, perspectives for future studies on crystal growth modeling are proposed. (c) 2008 National Natural Science Foundation of China and Chinese Academy of Sciences. Published by Elsevier Limited and Science in China Press. All rights reserved.

Effects and numerical analysis of argon gas flow on the oxygen concentration in Czochralski silicon single crystal growth

Argon gas, as a protective environment and carrier of latent heat, has an important effect on the temperature distribution in crystals and melts. Numeric simulation is a potent tool for solving engineering problems. In this paper, the relationship between argon gas flow and oxygen concentration in silicon crystals was studied systematically. A flowing stream of argon gas is described by numeric simulation for the first time. Therefore, the results of experiments can be explained, and the optimum argon flow with the lowest oxygen concentration can be achieved. (C) 2002 Elsevier Science B.V. All rights reserved.

Effects and numerical analysis of argon gas flow on the oxygen concentration in Czochralski silicon single crystal growth

Argon gas, as a protective environment and carrier of latent heat, has an important effect on the temperature distribution in crystals and melts. Numeric simulation is a potent tool for solving engineering problems. In this paper, the relationship between argon gas flow and oxygen concentration in silicon crystals was studied systematically. A flowing stream of argon gas is described by numeric simulation for the first time. Therefore, the results of experiments can be explained, and the optimum argon flow with the lowest oxygen concentration can be achieved. (C) 2002 Elsevier Science B.V. All rights reserved.

Void-like defects in annealed Czochralski silicon

Void-like defects of octahedron structure having {111} facets were observed in annealed Czochralski silicon. The amorphous coverage of SiOx and SiCx on the inner surface of the defects was identified using transmission electron microscopy and electron energy-loss spectroscopy. It is suggested that these defects are a kind of amorphous precipitate origin. A mechanism for the generation of these defects and the previously reported solid amorphous precipitates is proposed. (C) 1998 American Institute of Physics. [S0003-6951(98)02842-3].

PROTON-IMPLANTED NEW TYPE SILICON MATERIAL SUBJECTED TO 2-STEP FURNACE ANNEALING

A high-resistivity defect layer buried beneath the silicon surface layer by using proton implantation and two-step conventional furnace annealing is described. During the first annealing step (600-degrees-C), implanted hydrogen atoms move towards the damage region and then coalesce into hydrogen gas bubbles at the residual defect layer. During the second annealing step (1180-degrees-C) these bubbles do not move due to their large volume. Structural defects are formed around the bubbles at a depth of approximately 0.5-mu-m. The defect layer results in a high resistivity value. Experiments show that the quality of the surface layer has been improved because the surface Hall mobility increased by 20%. The sample was investigated by transmission electron microscopy.

Dislocation movement in nitrogen-doped Czochralski silicon

Dislocation movement in N-doped Czochralski silicon (Cz-Si) was surveyed by four point bend method. Dislocation movement velocities in Cz-Si doped with nitrogen, with both nitrogen and antimony, and with only antimony were investigated. The order of measured dislocation movement velocities, at 700 degrees C less than or equal to T less than or equal to 800 degrees C and under resolved stress sigma=4.1 kg/mm(2), was V-Sb.O > V-n.Sb.O>V-N.O. The experiments showed that nigtrogen doping could retard the movement of dislocations.

COMPARISON OF PROPERTIES OF SOLID-PHASE EPITAXIAL SILICON-ON-SAPPHIRE FILMS RECRYSTALLIZED BY RAPID THERMAL ANNEALING AND FURNACE ANNEALING

Chemically vapour deposited silicon on sapphire (SOS) films 0.25 mu m thick were implanted with Si-28(+) and recrystallized in solid phase by furnace annealing (FA) and IR rapid thermal annealing (RTA) in our laboratory. An improvement in crystalline quality can be obtained using both annealing procedures. After FA, it is hard to retain the intrinsic high resistivity value(10(4)-10(5) Ohm cm) observed in as-grown SOS films, so the improvement process cannot be put to practical use effectively. However, it is demonstrated that by properly adjusting the implantation and RTA conditions, significant improvements in both film quality and film autodoping can be accomplished. This work describes a modified double solid phase epitaxy process in which the intrinsic high resistivities of the as grown SOS films are retained. The mechanism of suppression of Al autodoping is discussed.

Oxygen precipitation in Czochralski grown silicon heat treated at 525 degrees C

The diffusion of interstitial oxygen In silicon at 525 degrees C is studied using time-of-flight small-angle neutron scattering (SANS) to separate the elastic scattering from oxygen-containing aggregates from the inelastic scattering from neutron-phonon interactions. The growth of oxygen-containing aggregates as a function of time gives a diffusion coefficient, D, calculated from Ham's theory, that is I factor of similar to 3.8 +/- 1.4 times higher than that expected by extrapolation of higher and lower temperature data (D = 0.13 exp(-2.53 eV kT(-1)) cm(2) s(-1)). This result confirms previous observations of enhanced diffusion at intermediate temperatures (400 degrees C-650 degrees C) although the magnitude of the enhancement we find is Much smaller than that reported by some others.

LSSA large area silicon sheet task continuous Czochralski process development : quarterly report no. 1, October-December, 1977 /

Includes bibliographical references.

Position-sensitive silicon strip detector characterization using particle beams

Silicon strip detectors are fast, cost-effective and have an excellent spatial resolution. They are widely used in many high-energy physics experiments. Modern high energy physics experiments impose harsh operation conditions on the detectors, e.g., of LHC experiments. The high radiation doses cause the detectors to eventually fail as a result of excessive radiation damage. This has led to a need to study radiation tolerance using various techniques. At the same time, a need to operate sensors approaching the end their lifetimes has arisen. The goal of this work is to demonstrate that novel detectors can survive the environment that is foreseen for future high-energy physics experiments. To reach this goal, measurement apparatuses are built. The devices are then used to measure the properties of irradiated detectors. The measurement data are analyzed, and conclusions are drawn. Three measurement apparatuses built as a part of this work are described: two telescopes measuring the tracks of the beam of a particle accelerator and one telescope measuring the tracks of cosmic particles. The telescopes comprise layers of reference detectors providing the reference track, slots for the devices under test, the supporting mechanics, electronics, software, and the trigger system. All three devices work. The differences between these devices are discussed. The reconstruction of the reference tracks and analysis of the device under test are presented. Traditionally, silicon detectors have produced a very clear response to the particles being measured. In the case of detectors nearing the end of their lifefimes, this is no longer true. A new method benefitting from the reference tracks to form clusters is presented. The method provides less biased results compared to the traditional analysis, especially when studying the response of heavily irradiated detectors. Means to avoid false results in demonstrating the particle-finding capabilities of a detector are also discussed. The devices and analysis methods are primarily used to study strip detectors made of Magnetic Czochralski silicon. The detectors studied were irradiated to various fluences prior to measurement. The results show that Magnetic Czochralski silicon has a good radiation tolerance and is suitable for future high-energy physics experiments.

Study of formation of silicon carbide in the Acheson process

Formation of silicon carbide in the Acheson process was studied using a mass transfer model which has been developed in this study. The century old Acheson process is still used for the mass production of silicon carbide. A heat resistance furnace is used in the Acheson process which uses sand and petroleum coke as major raw materials.: It is a highly energy intensive process. No mass transfer model is available for this process. Therefore, a mass transfer model has been developed to study the mass transfer aspects of the process along with heat transfer. The reaction kinetics of silicon carbide formation has been taken from the literature. It has been shown that reaction kinetics has a reasonable influence on the process efficiency. The effect of various parameters on the process such as total gas pressure, presence of silicon carbide in the initial charge, etc. has been studied. A graphical user interface has also been developed for the Acheson process to make the computer code user friendly.