Numerical Simulation of Heat Transfer and Kinetics in the Bulk Growth of Silicon Carbide

The growth process of 2-inch silicon carbide (SiC) single crystals by the physical vapor transport method (or modified Lely method) has been modeled and simulated. The comprehensive process model incorporates the calculations of radio frequency (RF) induction heating, heat and mass transfer and growth kinetics. The transport equations for electromagnetic field, heat transfer, and species transport are solved using a finite volume-based numerical scheme called MASTRAPP (Multizone Adaptive Scheme for Transport and Phase Change Process). Temperature distribution for a 2-inch growth system is calculated, and the effects of induction heating frequency and current on the temperature distribution and growth rate are investigated. The predicted results have been compared with the experimental data.

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.

Numerical Simulation of the Flow Field and Concentration Distribution in the Bulk Growth of Silicon Carbide Crystals

The physical vapor transport (PVT) method is being widely used to grow large-size single SiC crystals. The growth process is associated with heat and mass transport in the growth chamber, chemical reactions among multiple species as well as phase change at the crystal/gas interface. The current paper aims at studying and verifying the transport mechanism and growth kinetics model by demonstrating the flow field and species concentration distribution in the growth system. We have developed a coupled model, which takes into account the mass transport and growth kinetics. Numerical simulation is carried out by employing an in-house developed software based on finite volume method. The results calculated are in good agreement with the experimental observation.

Growth of Silicon Carbide Bulk Crystals by Physical Vapor Transport Method and Modeling Efforts in the Process Optimization

Silicon carbide bulk crystals were grown in an induction-heating furnace using the physical vapor transport method. Crystal growth modeling was performed to obtain the required inert gas pressure and temperatures for sufficiently large growth rates. The SiC crystals were expanded by designing a growth chamber having a positive temperature gradient along the growth interface. The obtained 6H-SiC crystals were cut into wafers and characterized by Raman scattering spectroscopy and X-ray diffraction, and the results showed that most parts of the crystals had good crystallographic structures.

Formation of dendritic nanostructures in Pyrex glass anodically bonded to silicon coated with an aluminum thin film

Anodic bonding with thin films of metal or alloy as an intermediate layer, finds increasing applications in micro/nanoelectromechanical systems. At the bonding temperature of 350 degrees C, voltage of 400 V, and 30 min duration, the anodic bonding is completed between Pyrex glass and crystalline silicon coated with an aluminum thin film with a thickness comprised between 50 and 230 nm. Sodium-depleted layers and dendritic nanostructures were observed in Pyrex 7740 glass adjacent to the bonding interface. The sodium depletion width does not increase remarkably with the thickness of aluminum film. The dendritic nanostructures result from aluminum diffusion into the Pyrex glass. This experimental research is expected to enhance the understanding of how the depletion layer and dendritic nanostructures affect the quality of anodic bonding. (C) 2007 Elsevier B.V. All rights reserved.

Crystal structure and phase transition studies in perovskite-type oxides using powder-diffraction techniques and symmetry-mode analysis : SrLnMRuO6 (Ln=La,Pr,Nd; M=Zn,Co,Mg,Ni,Fe) and ALn2CuTi2O9 (A=Ca,Ba; Ln=La,Pr,Nd,Sm)

La tesis se ha centrado en la síntesis y caracterización estructural de materiales tipo perovskita: SrLnMRuO6 (Ln=La,Pr,Nd; M=Zn,Co,Mg,Ni,Fe) y ALn2CuTi2O9 (A=Ca,Ba; Ln=La,Pr,Nd,Sm). El estudio de las estructuras de los materiales se ha realizado mediante el análisis de los patrones de difracción en polvo de rayos-X, sincrotrón y/o neutrones. En el refinamiento por el método de Rietveld de las estructuras se han sustituido las coordenadas atómicas (el método más común), por coordenadas colectivas: las amplitudes de los modos que describen la distorsión de la fase prototipo. Los resultados generales para la serie SrLnMRuO6 (Ln=La,Pr,Nd; M=Zn,Co,Mg,Ni) a temperatura ambiente se ha recogido en un diagrama en el que se han indicado las amplitudes de los modos que transforman de acuerdo a las irreps en función del factor de tolerancia, ya que todos ellos cristalizan en la misma fase monoclínica (P21/n); y a temperaturas altas se ha construido un diagrama de fase. Los materiales SrLnFeRuO6 ( Ln=La,Pr,Nd) y CaLn2CuTi2O9 cristalizan en la fase ortorrómbica Pbnm a temperatura ambiente; mientras que BaLn2CuTi2O9 tienen una estructura más simétrica, I4/mcm. A altas temperaturas se han identificado las transiciones de fase inducidas por el cambio de temperatura.A temperaturas bajas se han analizado las estructuras magnéticas de algunos de los compuestos mediante difracción de neutrones.