ø300 mm直拉硅单晶生长过程中的变晶现象及其影响因素

摘要/Abstract

摘要： 直拉法生长直径300 mm硅单晶过程中,直径均匀是获得高品质硅单晶的关键。在生产实践中发现,当硅晶体进入等径生长阶段,过高的提拉速度会引起晶体发生扭晶现象,导致晶线断裂随即变晶,对等径生长不利。本文采用数值模拟和理论相结合的方法分析了ø300 mm硅单晶生长过程中扭晶现象的成因,建立了不同提拉速度下晶体直径与熔体温度分布的关系,分析了晶体发生扭晶的影响因素。结果表明,随着提拉速度的增加,熔体自由表面产生过冷区且该过冷区随提拉速度的增加不断扩大,过冷区的产生是导致晶体发生扭晶的主要原因。提出了一种基于有限元热场数值模拟的最大稳定提拉速度的判别方法,并给出了通过改变晶体旋转速度来改善熔体自由表面温度分布的工艺措施建议,从而避免晶体扭晶现象的发生。研究结果对设计大尺寸硅单晶生长热场具有一定的指导作用。

关键词: 直拉硅单晶, 扭晶, 数值模拟, 提拉速度, 晶体旋转速度, 过冷区

Abstract: In the process of growing 300 mm diameter silicon single crystal by Czochralski method, uniform diameter is the key to obtain high quality silicon single crystal. It is found in production practice that too high pulling speed causes the phenomenon of twisting crystal at uniform diameter growth stage, which leads to the crystal line fracture and then crystalloblast, and it is unfavorable to uniform diameter growth. In this paper, the origin of the twisting crystal in the growth of ø300 mm silicon single crystal was analyzed by the combination of numerical simulation and theory analysis. For varying pulling speed, the correlation between crystal diameter and melt temperature distribution was established, the affecting factors on twisting crystal were obtained. The results show that with the increase of pulling speed, the supercooling zone is generated on the free surface of the melt, and continue to expand. The generation of supercooling zone is the main reason leading to twisting crystal. A method for identifying the maximum stable pulling speed is presented based on the numerical simulation of the finite element thermal field, and the method of varying the rotation speed of the crystal is proposed to improve the free surface temperature distribution of the melt, so as to avoid the crystal twisting. The results provide a guidance to the thermal field design of large diameter silicon monocrystalline growth.

Key words: Czochralski silicon monocrystal, twisting crystal, numerical simulation, pulling speed, crystal rotation speed, supercooling zone

中图分类号:

O782

张晶, 刘丁. ø300 mm直拉硅单晶生长过程中的变晶现象及其影响因素[J]. 人工晶体学报, 2022, 51(7): 1185-1193.

ZHANG Jing, LIU Ding. Crystalloblast and Its Affecting Factors in ø300 mm Czochralski Silicon Single Crystal Growth[J]. JOURNAL OF SYNTHETIC CRYSTALS, 2022, 51(7): 1185-1193.

参考文献

[1] KALAEV V V. Liquid flow in a cubic cavity generated by gas motion along the free surface[J]. International Journal of Heat and Mass Transfer, 2012, 55(19/20): 5214-5221.
[2] SCHWABE D, UECKER R, BERNHAGEN M, et al. An analysis of and a model for spiral growth of Czochralski-grown oxide crystals with high melting point[J]. Journal of Crystal Growth, 2011, 335(1): 138-147.
[3] JUNG Y J, KIM W K, JUNG J H. Process design and simulation for optimizing the oxygen concentration in Czochralski-grown single-crystal silicon[J]. Journal of the Korean Physical Society, 2014, 65(3): 362-367.
[4] NAM W, HAHN Y, BAIK S. Local optimization of graphite heater to save a power consumption of Czochralski Si ingot grower for PV application[J]. Energy Procedia, 2017, 124: 767-776.
[5] 高宇,朱亮,张俊,等.定拉速生长对ø300 mm直拉硅单晶生长影响分析[J].人工晶体学报,2020,49(5):811-814+823.
GAO Y, ZHU L, ZHANG J, et al. Analysis of ø300 mm Cz silicon single crystal growth with constant growth rate[J]. Journal of Synthetic Crystals, 2020, 49(5): 811-814+823(in Chinese).
[6] KALAEV V, SATTLER A, KADINSKI L. Crystal twisting in Cz Si growth[J]. Journal of Crystal Growth, 2015, 413: 12-16.
[7] TAKAGI K, FUKAZAWA T. Effect of growth conditions on the shape of Bi₄Ge₃O₁₂ single crystals and on melt flow patterns[J]. Journal of Crystal Growth, 1986, 76(2): 328-338.
[8] OKANO Y, IKEYA T, HOSHIKAWA K, et al. Behavior of RAlO₃ (R=Ho, Er and Dy) during crystal growth by the Czochralski method[J]. Journal of Crystal Growth, 1993, 131(3/4): 616-619.
[9] FEI Y T, CHOU M M C, CHAI B H T. Crystal growth and morphology of substituted gadolinium gallium garnet[J]. Journal of Crystal Growth, 2002, 240(1/2): 185-189.
[10] KAMADA K, YANAGIDA T, PEJCHAL J, et al. Scintillation properties of Ce doped Gd₂Lu₁(Ga, Al)₅O₁₂ single crystal grown by the micro-pulling-down method[J]. Journal of Crystal Growth, 2012, 352(1): 35-38.
[11] 李景湧.有限元法[M].北京:北京邮电大学出版社,1999.
LI J Y. Finite element method[M]. Beijing: Beijing University of Posts and Telecommunications Press, 1999(in Chinese).
[12] TU H L, XIAO Q H, GAO Y, et al. Numerical analysis and simulation of Czochralski growth processes for large diameter silicon crystals[J]. Rare Metals, 2007, 26(6): 521-527.
[13] VIRZI A. Finite element analysis of the thermal history for Czochralski growth of large diameter silicon single crystals[J]. Journal of Crystal Growth, 1989, 97(1): 152-161.
[14] BORNSIDE D E, KINNEY T A, BROWN R A. Minimization of thermoelastic stresses in Czochralski grown silicon: application of the integrated system model[J]. Journal of Crystal Growth, 1991, 108(3/4): 779-805.
[15] KIRPO M. Global simulation of the Czochralski silicon crystal growth in ANSYS FLUENT[J]. Journal of Crystal Growth, 2013, 371: 60-69.
[16] 张向宇.ø200 mm和ø450 mm直拉单晶硅生长过程有限元模拟研究[D].济南:山东大学,2015.
ZHANG X Y. Finite Element Simulation of ø200 mm and ø450 mm Czochoraski single silicon crystal growth[D]. Jinan: Shandong University, 2015(in Chinese).
[17] ZHANG J, LIU D, PAN Y N. Suppression of oxygen and carbon impurity deposition in the thermal system of Czochralski monocrystalline silicon[J]. Journal of Semiconductors, 2020, 41(10): 102702.
[18] 刘丁.直拉硅单晶生长过程数值模拟与工艺优化[M].北京:科学出版社,2020.
LIU D. Numerical simulation and process optimization of Czochralski silicon single crystal growth[M]. Beijing: Science Press, 2020(in Chinese).
[19] 滕冉,戴小林,肖清华,等.大直径硅单晶生长过程中固/液界面形状及熔体流动的数值分析[J].人工晶体学报,2013,42(4):611-615.
TENG R, DAI X L, XIAO Q H, et al. Numerical simulation of melt/crystal interface and melt flow during large diameter single silicon crystal growth[J]. Journal of Synthetic Crystals, 2013, 42(4): 611-615(in Chinese).
[20] ZHANG J, LIU D, ZHAO Y, et al. Impact of heat shield structure in the growth process of Czochralski silicon derived from numerical simulation[J]. Chinese Journal of Mechanical Engineering, 2014, 27(3): 504-510.
[21] LIU B T, YU Y, TANG X, et al. Influence of silicon melt convection on interface instability in large-size silicon carbide solution growth[J]. Journal of Crystal Growth, 2019, 527: 125248.
[22] 闵乃本.晶体生长的物理基础[M].上海:上海科学技术出版社,1982.
MIN N B. Physical basis of crystal growth[M]. Shanghai: Shanghai Science and Technology Press, 1982(in Chinese).
[23] TAO G, JIN Z. A daptive cascade generalized predictive control[C]. International Journal of Intelligence Science, 2014, 4(7): 70-79.
[24] 金超花,朱彤.旋转对Cz法硅晶生长过程中热输运和熔体流动影响的数值模拟[J].人工晶体学报,2015,44(1):31-37.
JIN C H, ZHU T. Numerical simulation on the effect of the rotating on the heat transfer and melt flow during the silicon crystal growth process by Cz method[J]. Journal of Synthetic Crystals, 2015, 44(1): 31-37(in Chinese).