Welcome to Visit BULLETIN OF THE CHINESE CERAMIC SOCIETY! Today is

BULLETIN OF THE CHINESE CERAMIC SOCIETY ›› 2025, Vol. 44 ›› Issue (9): 3451-3461.DOI: 10.16552/j.cnki.issn1001-1625.2025.0172

• Functional Materials • Previous Articles     Next Articles

Numerical Simulation of Interface Shape and Constitutional Supercooling on Heavily Phosphorus-Doped Czochralski Silicon with Ultra-Low Resistivity

MA Wuxiang1, MEI Haotian1, LI Xiaochuan2, FAN Jixiang1, WU Yanguo1   

  1. 1. MCL Electronic Materials Co., Ltd., Luoyang 471000, China;
    2. Longmen Laboratory, Luoyang 471000, China
  • Received:2025-02-19 Revised:2025-03-20 Online:2025-09-15 Published:2025-09-19

Abstract: The production of heavily phosphorus-doped silicon single crystals with ultra-low resistivity is primarily constrained by constitutional supercooling induced by high-concentration doping and crystal-melt interface deformation during Czochralski growth. In this study, to investigate the persistent lattice mismatch observed at a crystal length of approximately 1 000 mm during production, simulations were performed using CGSIM finite element software. It was demonstrated that the lattice mismatch originated from the elevation of the crystal-melt interface and phosphorus impurity segregation. And an improvement was proposed. The results indicate that interface height can be reduced by lowering crystal rotation, growth rate, and melt gap, while increasing crucible rotation. Based on the constitutional supercooling criterion V/G (crystal growth rate V divided by melt temperature gradient G), crystal pulling is conducted with a reduced V value to mitigate supercooling effects. By applying optimized parameters, dislocation-free 8-inch (1-inch=2.54 cm) heavily phosphorus-doped silicon single crystals are successfully grown over 5 consecutive cycles. The measured resistivity ranges from 0.000 93 Ω·cm to 0.001 25 Ω·cm, further validating the simulation results.

Key words: ultra-low resistivity, heavily phosphorus doping, silicon single crystal, numerical simulation, interface shape, constitutional supercooling

CLC Number: