Search for content and authors
 

Simulation of ultrasound influence on melt convection for the growth of striation-free GaxIn1-xSb and Si single crystals by the Czochralski method

Gennadiy N. Kozhemyakin 1Lyudmila Nemets Ann Bulankina 1

1. Volodymyr Dahl East Ukrainian National University, Luhansk 91050, Ukraine

Abstract

It has been shown previously that the presence of ultrasound parallel to the pulling axis in the melt decreased growth striations in semiconductor single crystals grown by the Czochralski method. The effect of ultrasound on striations, which decrease in number and disappear in the central part of InSb, GaAs, BixSb1-x and GaxIn1-xSb solid solutions pulled single crystals, was studied and reported in several publications of our group [1-4]. We showed experimentally that ultrasonic waves at frequencies from 0.15 to 5 MHz can be used for damping of convection in the melt during crystal growth. The damping of convection is attributed to forming of ultrasonic standing waves under the sold-liquid (S/L) interface. However, the behavior of convective flow under the effect of ultrasound was not studied.

In the present study, we modeled convection at the growth conditions of small GaxIn1-xSb and large Si single crystals. The flow simulation for GaxIn1-xSb and Si melts was conducted for quasi-steady conditions, when Raleigh number was less than 105. For these conditions applied to Ga0.03In0.97Sb and Si single crystals, the growth parameters have the values listed in Table 1.

Table 1. Parameters of the growth crystals

Variable Symbol Crystal
Ga0.03In0.97Sb Si
Crystal radius, mm  r 5 75
Crucible radius, mm Rc 15 150
Melt height, mm h 23 40
Crystal rotation rate, rpm Ω 1 ­– 20 5 – 10
Maximum velocity of flow, cm/s V 0.047–0.9 12.6 – 28
Particle radius, nm rp 5 5
Particle mass, kg Pp 3×10-21 1,3×10-21
To formulate the mixed convection problem these melts, a viscous incompressible flow in a cylindrical crucible was considered. The growing crystals were pulled very slowly from the center of the stationary crucible, and the vertical motion of the crystals was assumed to be negligible. The temperature of crucible walls and bottom surfaces was fixed at a constant value. The Navier-Stokes equation was used for the calculations of convective flow velocity in GaxIn1-xSb and silicon melts for quasi-steady conditions. The flow was assumed to remain laminar and incompressible, with constant fluid properties.

The radial and axial distributions of flow in GaxIn1-xSb and silicon melts were calculated. The radial velocity distribution of flow increases from the center to Rc/2 in both crucibles, and then decreases. The maximum velocity of axial flow was under the S/L interface at the axis of the crucibles. However, summary values of the radial and axial flow velocities had maximum values under the S/L interface near periphery of the crystals. The maximum velocities were 0.9 cm/s and 28 cm/s for GaxIn1‑xSb and silicon melts, respectively.

The maximum calculated flow velocity was used for the calculation of the force Fc acting on Si particles in the melt. Particles of 5 nm radius and different mass for these melts were accepted for our calculations. Then, the maximal force FC acting on the particle in convective flow was 3∙10-21 N and 1.2∙10-21 N for GaxIn1-xSb and Si melts, respectively.

It is established that introduction of ultrasound into the liquid forms a standing wave channel under the S/L interface. The force FUS acting on the melt particle in ultrasonic standing waves was calculated according to the equation in [4]. The ultrasonic field parameters for a given system were selected at frequencies of 0.1 – 4.0 MHz. The force FUS increases with increase of the frequency in GaxIn1-xSb and Si melts from 1∙10-20 N to 7∙10-19 N and 6.3∙10-23 N to 4.0×10-18 N, respectively.

The results of these calculations have shown that the ultrasound at frequencies over 0.1 MHz can be used for reduction of striations in pulled GaxIn1-xSb single crystals, since the force FUS is larger than the convective force FC by a factor of 3 – 600 in this melt. However, for striation-free Si single crystals, the ultrasound frequencies must be larger than 0.5 MHz, according to our calculations.

Therefore, we believe that ultrasound at a high frequency can be used for reduction of convection in Czochralski growth of GaxIn1-xSb solid solutions and Si single crystals. The strong reduction of flow motion is associated with the formation of a standing wave channel between the S/L interface and the crucible bottom. This modeling has shown the potential of using ultrasound to damp convection in the melt and decrease striations in semiconductor single crystals. 

References

[1] G.N. Kozhemyakin,V.G. Kosushkin, S.Y. Kurochkin, J. Crystal Growth. 121, 240 (1992).

[2] G.N. Kozhemyakin,L.G. Kolodyazhnaya, J. Crystal Growth. 147, 200 (1995).

[3] G.N. Kozhemyakin, J. Crystal Growth. 149, 266 (1995).

[4] G.N. Kozhemyakin, J. Crystal Growth. 257, 237 (2003).  

 

Legal notice
  • Legal notice:
 

Related papers

Presentation: Poster at 17th International Conference on Crystal Growth and Epitaxy - ICCGE-17, Topical Session 8, by Gennadiy N. Kozhemyakin
See On-line Journal of 17th International Conference on Crystal Growth and Epitaxy - ICCGE-17

Submitted: 2013-03-19 21:48
Revised:   2013-07-18 21:29